AD  A 0 3 2 1 2 8 


AFAL-TR-76-U9 


IDENTIFICATION  OF  DEFECTS  IN  COMPOUND 
SEMICONDUCTORS 


UNIVERSITY  OF  DAYTON 
PHYSICS  DEPARTMENT 
WO  COLLEGE  PARK 
DAYTON,  OHIO  454 <99 


OCTOBER  1976 


TECHNICAL  REPORT  AFAL-TR-76-189 

FINAL  REPORT  FOR  PERIOD  1 JULY  1971  - 81  MARCH  1976 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  AVIONICS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  46488 


NOTICE 


. 

Mmb  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  In  connection  with  a definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  Incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  In 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  Implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture, 
use,  or  sell  any  patented  Invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS),  At  NTIS,  It  will  be  available 
to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


PROJECT  ENGINEER 


For  the  Coamander 


T-  . 


* ^ /> 

DONALD  C.  REYNOLDS 
Senior  Scientist 
Electronics  Research  Branch 
Electronic  Technology  Division 
Air  Force  Avionics  Laboratory 


ROBERT  D.  LARSON 

Chief,  Electronics  Research  Branch 
Electronic  Technology  Division 
Air  Force  Avionics  Laboratory 


Copies  of  this  report  should  not  be  returned  unless  return  Is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a specific  document. 

AW  FORCfc  - 8 NOVO«R  78  - 125 


AFAIHTR-76-139  j- 


T a PERIOD  COVERED 


Final 
1 Jul. 


RACT  OR  GRANT  NUMBER!*) 


A ABSTRACT  ( Continue  on  reverse  side  If  necessary  end  Identify  by  block  number) 

Kadiat ion-damaged, ion- implanted,  doped,  and  as-grown  semiconducting  materials, 
including  CdS,  CdSe,  CdTe,  CdO,  ZnS,  ZnSe,  ZnTe,  ZnO,  HgS,  CuInS^,  In„Te!T, 

GaAs,  GaAlAs,  Si,  and  Ge,  have  been  examined  by  photoluminescence,  cathodolumin- 
escence,  optical  absorption,  thermal  conductivity,  magnetic  susceptibility. 

Hall  effect,  photoconductivity,  photomagnetoelectric  effect,  and  nuclear 


DO  , 1473  COITION  OF  I NOV  65  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


3ArF  READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

2 GOVT  ACCESSION  NO.  3 RECIPIENT’S  CATALOG  NUMBER 


4.  TITLE  (and  Subtitle)  

IDENTIFICATION  OF  DEFECTS  IN  COMPOUND 
SEMICONDUCTORS,  * * 


17.  AUTHORS) 


David  C./Look,  Jon  M./Meese,  Marshall  M./Kreitman 


I Jose  C./Manthuruthil 


».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  10  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREA  A WORK  UNIT  NUMBERS 

Physics  Dept. 

University  of  Dayton  78850613 

Dayton,  OH  45469 ^ 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Air  Force  Avionics  Laboratory/DHR  Oct^B(!flB76 

Wright-Patterson  AFB,  OH  45433  u number  of  pages — - 

Contract  Monitor:  Donald  C.  Reynolds 214  _/  / / 

IA.  MONITORING  AGENCY  name  » ADDRESS!/!  different  Irani  Controllint  Oltice)  Ts  SECURITY  CLTSV%!  fWi  repot 4> 


Unclassified 

I5».  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


fl?^  DISTRIBUTION  STATEMENT  (o I Ihie  Report) 


A — Aporoved  for  public  release;  distribution  unlimited. 


I 17  DISTRIBUTION  STATEMENT  (ol  the  ebetrect  entered  in  Block  20,  it  dllterenl  Irom  Report) 


fie  supplementary  notes 


1 19.  KEY  WORDS  ( Continue  on  reverse  side  il  necessary  and  Identify  by  block  number) 


III-V  Compounds 
II-VI  Compounds 
I-III-V1  Compounds 


Luminescence 
Optical  Absorption 
Thermal  Conductivity 


ISemi-insulating  GaAs  Hall-effect 


■Radiation  Damage 


Nuclear  Magnetic  Resonance 


Solar  Cells 
Photoconductivity 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Entered) 


FOREWORD 


' 


The  contract  effort  described  in  this  report  represents  approximately 
17  man-years  of  effort,  from  1 July  1971,  to  31  March  1976.  At  the  begin- 
ning of  this  period  the  work  was  primarily  concerned  with  electron-irradiation 
damage  in  II-VI  compounds.  Toward  the  end  of  the  period  the  chief  interest 
was  in  the  characterization  of  III-V  compounds,  especially  GaAs.  Thus,  much 
of  our  time  has  been  spent  in  establishing  new  techniques  and  becoming 
familiar  with  different  materials.  In  spite  of  this,  the  effort  has  been 
productive,  resulting  in  30  published  papers  in  refereed  journals. 

The  personnel  who  have  served  in  a full-time  capacity,  for  at  least 
part  of  the  contract  period,  include:  Dr.  David  C.  Look,  Mr.  Marshall  M. 

Kreitman,  Dr.  Jon  M.  Meese,  Dr.  Jose  C.  Manthuruthil,  and  Miss  Mary  Turner. 
Readers  who  desire  more  information  about  a particular  section  of  this  report 
should  contact  Dr.  Meese  (Section  II),  Dr.  Look  (Sections  III,  IV,  and  V), 
or  Mr.  Kreitman  (Sections  VI  and  VII). 

We  wish  to  thank  our  many  colleagues  at  the  Aerospace  Research  Laboratories 
(now  disestablished) , and  the  Avionics  Laboratory,  who  have  unselfishly  helped 
to  make  our  experience  here  pleasant  and  productive.  We  hope  that  our  research 
has  aided,  and  will  continue  to  aid.  Air  Force  efforts  to  solve  relevant  semi- 
conductor materials  problems. 


This  work  was  performed  by  the  University  of  Dayton  Physics  Department, 
300  College  Park,  Dayton,  Ohio,  45469,  under  Contract  F33615-71-C-1877.  The 
project,  task,  and  work  unit  number  was  78850613,  and  the  project  engineer  was 
Donald  C.  Reynolds,  AFAL/DHR,  Wright-Patterson  AFH,  Ohio,  45433.  The  report 


was  submitted  on  1 May  1976. 
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Semiconductor  devices  play  a major  and  increasingly  important  role  in 
many  catagories  of  Air  Force  systems.  Thus,  it  is  important  to  understand  the 
materials  out  of  which  these  devices  are  made,  and  our  contract  effort  has 
primarily  been  concerned  with  the  investigations  of  impurities  and  defects  in 
as-grown,  doped,  radiation-damaged,  and  ion-implanted  crystals.  The  semi- 
conductor materials  we  have  investigated  include  CdS,  CdSe,  CdTe,  CdO,  ZnS, 
ZnSe,  ZnTe,  ZnO,  HgS,  CuInS2,  In2Te^,  GaAs,GaAlAs,  Si,  and  Ge,  while  the 
techniques  we  have  employed  include  photoluminescence,  cathodoluminescence , 
optical  absorption,  thermal  conductivity,  magnetic  susceptibility.  Hall  effect, 
photoconductivity,  photomagnetoelectric  effect,  and  nuclear  magnetic  resonance. 

During  the  first  2-3  years  of  this  contract  effort  a primary  goal  was  the 
solution  of  radiation-damage  problems  in  II-VI  compounds,  and  this  work  is 
described  in  Section  II.  During  this  period,  several  important  atomic- 
displacement  threshold  energies  were  measured.  In  Section  III  we  describe 
some  nuclear-magnetic-resonance  investigations  of  the  properties  of  free  and 
bound  electrons  in  CdS,  CdTe,  and  GaAs.  Among  the  accomplishments  here  was  the 
measurement  of  the  effective  g-factor  of  the  conduction  electrons  in  CdTe, 
an  important  parameter  for  tunable,  spin-flip  Raman  lasers. 

One  of  the  promising  new  solar-cell  materials  is  CuInS2  and  we  were  able 
to  grow  perhaps  the  largest  single  crystals  presently  available.  Furthermore, 
we  carried  out  tha  first  extensive  investigation  of  the  electrical  properties 
of  this  material,  both  n-  and  p-type.  This  work  is  discussed  in  Section  IV. 

In  Section  V we  consider  the  problems  associated  with  GaAs  substrate  materials. 
We  have  devised  a scheme  for  determining  the  donor  and  deep  acceptor  (Cr) 
concentrations  from  the  mixed-conductivity  electrical  data. 
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The  thermal  conductivity  of  a semiconductor  material  is  an  important  parameter, 
especially  for  cases  in  which  the  device  will  be  used  at  very  high  or  very 
low  temperatures.  A relatively  fast  and  convenient  instrument  which  has  recently 
been  developed  is  the  thermal  comparator,  and  we  have  studied  its  applicability 
to  the  II-V1  compounds  and  other  related  materials.  It  has  proved  to  be  a 
useful  and  effective  tool,  as  shown  in  Section  VI.  In  Section  VII  we  look  at 
the  optical  properties  of  several  materials,  including  HgS  and  high-resistivity 
GaAs.  HgS  has  great  potential  in  non-linear  optics  applications,  and  we  have 
grown  and  studied  some  large,  single  crystals  of  a-HgS.  In  the  GaAs  effort 
we  have  been  able  to  clear  up  some  recent  problems  in  the  interpretation  of 
photoconductivity  data. 

Since  most  of  the  work  described  in  this  report  is  already  in'  the  published 
literature  we  have  not  gone  into  great  detail  in  many  areas.  However,  we  have 
attempted  to  carefully  document  the  research  so  that  the  interested  reader 
can  either  refer  to  the  literature  or  contact  the  authors  for  more  information. 
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2.0  ELECTRON  RADIATION  DAMAGE  IN  SEMICONDUCTING  COMPOUNDS 

2.1  Introduction 

At  the  beginning  of  this  contract,  the  accelerator  facilities  at 
AFARL  (now  AFAL/DHR)  were  to  a large  degree,  non-  functional . This  facility 
consisted  of  a 150  keV  Cockroft-Walton  accelerator  which  had  not  been  run 
for  at  least  five  years  and  from  which  most  of  the  beamline  vacuum  equipment 
and  electronics  had  been  stripped,  a 1 MeV  electron-beam  Van  de  Graaff  which 
would  produce  a maximum  beam  energy  of  only  about  400  keV  after  considerable 
arcing  and  vacuum  excursions,  and  a 400  keV  electron- beam  Van  de  Graaff 
which  was  largely  disassembled.  Furthermore,  a small  machine  shop  was  located 
in  the  accelerator  area  causing  considerable  difficulty  from  its  associated 
dust,  metal  chips  and  debris.  Shielding  for  the  accelerators  was  inadequate 
in  that  10-20  mR/hr  radiation  leakage  could  be  detected  near  the  vault  doors 
(maximum  permissible  human  exposure  is  300  mR/week).  If  the  machines  had 
been  operating  on  a continuous  40  hr  week,  a considerable  radiation  hazard  v;ould 
have  existed  for  the  operating  personnel  near  the  vault  doors.  Furthermore, 
the  vault  walls  did  not  extend  up  to  the  ceiling  and  approximately  10  mR/hr 
radiation  bounce  over  the  walls  made  it  impossible  to  work  in  the  Cockroft- 
Walton  vault  if  either  of  the  Van  de  Graaffs  was  running. 

The  accelerators  had  only  a single  beam  port  each,  since  there  were  no 
beam-switching  magnets  installed,  making  it  impossible  to  set  up  multiple 
experiments  on  the  beam  lines.  Accelerator-beam-energy  calibration  had  been 
accomplished  only  approximately  by  measuring  the  range  of  electrons  in  A1 
foil  of  various  thickness.  Although  a sound  technique  in  principle,  the 
theoretical  calculations  of  A.  T.  Nelms  (N.B.S. -Circular  577)  had  been  used  for 
energy  calibration  rather  than  the  Katz-Penfold  experimental  values,  leading  to 
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experimental  errors  in  beam  energy  as  large  as  15%.  Since  none  of  the  accelerators 

had  beam-scanning  circuits,  it  was  nearly  impossible  to  define  fluence  (number 

2 

of  incident  electrons/cm  )due  to  beam-profile  nonuniformities.  Furthermore, 
beam  current  integrators,  a universally  accepted  fluence  measuring  instrument 
for  charged-particle  beams,  were  not  being  used. 

From  these  meager  beginnings  and  in  the  time  span  of  the  present  contract, 
the  accelerator-based  materials-research  program  and  facilities  at  WPAFB  have 
become  internationally  recognized  in  the  areas  of  radiation  damage  in  compound 
semiconductors,  ion  implantation  in  semiconductors,  and  radiation  effects  in 
semiconductor  devices  as  evidenced  by  the  many  acceptances  of  refereed  publica- 
tions at  international  radiation  effects  conferences,  international  ion- 
implantation  conferences, and  the  I.E.F..E.  Annual  Space  and  Nuclear  Radiation 
Effects  conferences.  As  a result  of  the  efforts  of  both  military  and  civilian 
personnel  in  this  group,  contract  personnel  from  Systems  Research  Laboratories, 
and  personnel  from  the  University  of  Dayton  employed  on  the  present  contract, 
a research  program  has  been  established  which  will  be.  a valuable  asset  in  solving 
Air  Force  materials  problems  for  years  to  come.  A broad  base  of  experience 
has  been  developed  in  the  research  group  ranging  from  elemental-semiconductor 
irradiation  problems  through  the  III-V  and  II-VI  semiconductor  materials, 
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areas  where  industrial  and  military  experience  is  still  maturing.  Furthermore, 
the  work  has  extended  from  vary  basic  research  problems,  such  as  annealing- 
kinetics  studies,  to  extremely  pertinent  and  applied  areas  such  as  radiation 
hardening  of  CaAs  and  Si  devices. 

In  addition  to  the  sometimes  necessarily  dull  and  dry  reporting  of  the 
experimental  facts,  we  are  hopeful  that  this  report  will  capture  some  of  the 
excitement  which  we  all  have  experienced  during  the  development  of  this  group 
under  conditions  which  were  frequently  adverse.  It  is  further  hoped  that  the 
results  of  the  present  program  will  stimulate  the  desire  in  others  tc  pursue 


4 


k 


with  renewed  vigor  those  elusive  semiconductor  defect  structures  which  tend 
more  often  than  not  to  play  havoc  with  the  best  of  electronic  systems  exposed 
to  a hostile  radiation  environment. 

2 . 2 Facilities  Modifications  and  Kxperimental  Fquipment  Design 

Considerable  modification  of  the  accelerator  facilities  at  UPAFB  has 
been  performed  under  this  contract  both  directlv  and  by  means  of  suggestions 
and  planning  originated  by  personnel  supported  by  this  contract.  The  purpose 
of  these  modifications  has  been  to  systematically  broaden  the  basis  of  experi- 
mental materials  problems  which  require  accelerator  usage  at  WPAFB.  As  a 
result  of  these  efforts  by  ourselves  and  others  in  the  group,  the  accelerator 
facility  at  AFAL/DHR  represents  perhaps  the  most  versatile  and  complete  facility 
in  the  mid-west  devoted  exclusively  to  materials  problems  of  interest  to  the 
Air  Force.  The  facilities  are  now  suitable  for  a variety  of  tasks  from  basic 
research  on  point-defect  interactions  and  defect-production  mechanisms  to 
the  applied  research  areas  of  device  fabrication  by  ion  implantation,  radiation 
testing  and  hardening  of  semiconductor  devices,  and  accelerator-based  trace- 
impurity  analysis. 

Considerable  effort  has  been  expanded  on  the  maintenance  and  repair  of  the 
accelerators  by  personnel  supported  by  this  contract.  The  initially  poor 
performance  of  the  1 MeV  Van  de  Graaff  in  achieving  a maximum  energy  of  only 
400  keV  before  arcing  was  attributed  to  a lack  of  accelerator-tube  conditioning 
or  a dirty  accelerator  tube.  After  a systematic  conditioning  by  operating  the 
machine  in  the  instability  range  of  energies  just  before  breakdown,  at  high 
beam  currents,  a maximum  energy  of  850  keV  was  achieved.  It  was  determined  that 
a dirty  accelerator  tube  was  the  primary  limitation  in  achieving  a 1 MeV  beam. 

A ohysical  inspection  indicated  that  the  tube  was  covered  with  a layer  of  Hg, 
probably  as  a result  of  contamination  from  the  Hg-dlf fusion  vacuum  pump  which 
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had  at  one  time  been  used  as  the  primary  pumping  system  for  the  beam  line. 

This  contamination  problem  was  due  either  to  improper  maintenance  of  the 
diffusion-pump  cold  trap  or  a vacuum  rupture  of  the  Al-foil  beam  windows, 
although  improper  maintenance  is  the  more  likely  since  a properly  trapped 
diffusion  pump  should  have  stopped  nearly  all  the  Hg  migration  into  the 
accelerator  tube.  The  Hg  was  actually  standing  in  puddles  in  the  acceleration 
tube,  shorting  out  many  of  the  acceleration  planes.  A new  accelerator 
tube  was  procured  and  installed. 

In  addition  to  the  installation  of  the  new  accelerator  tube,  the  belt- 
charging  system  was  changed  from  a charging  comb,  which  depended  upon  estab- 
lishing a corona  between  the  belt  and  the  comb  points,  to  a wiper  screen.  This 

modification  allowed  a stiffening  of  the  insulating-gas-mixture  dielectric  by 
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the  addition  of  about  5 lb/in  of  SF^.  After  all  these  modifications,  the 
maximum  energy  obtained  from  the  1 MeV  machine  was  1.50  MeV,  a million-volt 
gain  from  the  accelerator's  original  condition  at  the  start  of  the  contract, 
and  a 500  keV  gain  over  its  maximum  nominal  voltage  rating.  The  maximum 
voltage  attainable  is  now  limited  by  arcing  between  the  acceleration  planes  inside 
the  accelerator  tube  and  cannot  be  increased  without  lengthing  the  acceleration 
column  and  tank.  It  is  possible  to  make  these  modifications  if  additional 
funds  are  provided,  at  which  point  beam  energies  in  excess  of  2 MeV  should 
be  obtainable. 

Once  the  accelerator  was  functioning  properly,  consideration  was  given 
to  the  other  limitations,  such  as  lack  of  beam  scan,  lack  of  accurate  energy 
calibration,  lack  of  beam  steering,  and  the  lack  of  a sufficient  number  of 
beam  ports  to  accommodate  the  accelerator  usage  envisioned.  A homebuilt  beam- 
switching  magnet  was  installed  on  the  end  of  the  drift  tube  (Fig.  2-1).  The 
circuit  for  the  magnet-current  power  supply  is  shown  in  Fig.  2-2.  The  colls 
were  wound  on  two  6 1/2"  dia.  x 2 1/2"  spools  with  #16  wire.  The  25f2  - 50  W 
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Fig.  2-2.  Switching-magnet  power  supply  and  degaussing  circuit 
for  1 MeV  Van  de  Graaff. 
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resistor  provides  a leakage  path  for  the  back-EMF  when  the  magnetic  field  is 
collapsed.  The  Lambda  power  supply  was  run  as  a voltage  regulator  since  no 
temperature  stability  problems  were  encountered.  An  auto  transformer  (Fig.  2-2) 
can  be  switched  across  the  magnet  coils  for  degaussing  the  residual  magnetism 
in  the  pole  iron  pieces.  It  is  necessary  to  remove  this  field  when  making 
irradiations  using  the  center  beam  port. 

The  beam-energy  calibration  procedure  is  based  upon  measuring  the  range 
of  electrons  in  A1  foils  of  calibrated  thicknesses  as  a function  of  beam- 
energy  dial  setting.  The  dial  setting  is  the  setting  of  the  ten-turn  pot 
shown  in  Fig.  2-3  for  the  energy  stabilization  circuit.  In  this  circuit  the 
signal  from  the  accelerator  generating  voltmeter  is  compared  with  a reference 
voltage  obtained  from  a constant-voltage  supply  and  ten-turn  pot.  The  output 
of  the  voltage  difference  between  these  two  signals  (using  a Brown  recorder 
differential  amplifier  and  power  amplifier)  drives  a motor  which  in  turn  drives 
the  accelerator  belt-charge  control. 


The  Katz-Penfold  range-energy  relation  for  electrons  in  A1  (Table  2-1) 
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is  well  known  experimentally.  A set  of  secondary  standard  foils  with  various 
thicknesses  covering  the  energy  range  from  200  keV  to  1 MeV  can  be  mounted  in  the 
window  of  a Faraday  cup.  A 375-volt  bias  is  applied  to  the  cup  to  suppress 
secondary-electron  currents.  As  the  beam  energy  is  raised,  the  current  from 
the  Faraday  cup  increases  dramatically  once  the  electron  range  exceeds  the  foil 
thickness  (Fig.  2-4).  By  extrapolation  to  the  baseline,  the  dial  setting  for 
the  energy  of  electrons  which  just  penetrate  the  foil  is  determined. 

This  calibration  procedure  was  cross  checked  by  comparing  the  dial  setting 
vs.  beam  energy  curve  to  dial  settings  corresponding  to  various  channel  numbers 
of  a multichannel  analyser  which  was  attached  to  a Si(Li)  detector  and  pre-amp 
sampling  the  accelerator  beam.  The  channel  energies  were  calibrated  using  the 
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Fig.  2-3.  Beam-energy  stabilization  schematic  for  1 MeV  Van  de  Graaff. 
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Table  2-1.  Katz-Penfold  range-energy  values  for  electrons  in 

A1  [Katz  and  Penfold,  Rev.  Mod.  Phys.  24,  28  (1952)]. 


E(MeV) 

2 

R(tng/  cm  ) 

0.01 

0.016078 

0.02 

0.067863 

0.03 

1.5101 

0.04 

2.6135 

0.05 

3.9561 

0.06 

5.5121 

0.07 

0.08 

0.09 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

0.65 

0.70 

0.75 

0.80 

0.85 

0.90 

0.95 
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Fig.  2*4.  Range  measurement  for  a single  A1  foil  thickness. 
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same  detector  and  the  known  emission  energies  of  internal-conversion  electrons 
207  207 

from  Bi  A typical  Bi  spectrum  is  shown  in  Fig.  2-5.  The  agreement 


between  these  two  calibration  procedures  was  within  one  channel  number  (3  keV) 


or  less  than 


Figure  2-6  shows  the  agreement  between  the  beam  energy  as  measured  using 


the  multichannel  analyser  and  solid  state  detector  (solid  line)  and  the  energy 


calibration  obtained  by  using  the  Katz-Penfoid  experimental  range-energy 


relationship  and  the  ranges  as  determined  in  Fig.  2-4  for  various  foil  thick 


nesses.  The  agreement  is  excellent  and  we  believe  essentially  correct  since 


the  energies  are  determined  by  two  separate  and  independent  experiments 


Also  shown  in  Fig.  2-6  is  the  calculated  range-energy  relation  by  A.  T.  Nelms 


(broken  line)  along  with  the.  energy  labels  found  on  the  A1  foils  at  the 


beginning  of  this  contract 


accelerator  prior  to  this  contract  using  the  theoretical  Nelms  curves  rather 


than  the  experimental  Katz-Penfold  curves.  As  a result  of  this  error,  it 


is  likely  that  the  displacement  thresholds  reported  in  the  literature  by 


Kulp,  et.  al,  prior  to  1971  are  in  error.  Figure  2-6  can  be  used  to  correct 


these  values.  As  a result  of  this  error  the  displacement  thresholds  in  CdS 


have  been  remeasured  on  this  contract.  The  results  are  summarized  in  Table  2-2 


Table  2-2  Displacement  thresholds  in  CdS 


(Kulp.  et.al.) 
(Bryant  & Cox) 
(Eisby  (>  Meese) 
(Elsby  & Meese) 
(Kulp) 

(Elsby  & Meese) 


Sulfur 

Sulfur 

Sul  fur 

Sulfur 

Cadmium 

Cadmium 


The  older  data  from  ARL  (Kulp,  et.al.)  appear  to  be  systematically  lower 


than  the  more  recent  results  of  Elsby  and  Meese  and  the  data  of  Cox  which  are  in 
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Fig.  2-6.  Calibration-foils  7 ange  (thickness) -energy  curves.  The  solid 
curve  is  the  correct  one.  The  Katz-Penfold  data  may  be  found 
in  Ref.  1,  and  the  A.  T.  Nelms  data,  Ref.  2. 


agreement.  This  agrees  reasonably  well  with  the  curves  shown  in  Fig.  2-6, 
although  an  even  larger  difference  would  be  expected  for  the  cadmium  thresholds 
than  is  reported.  Since  no  beam  scanning  was  used  on  the  1 MeV  Van  de  Graaff 
prior  to  the  contract,  some  of  the  discrepancy  could  also  be  accounted  for  in 
terms  of  auto-focusing  of  the  beam  at  higher  accelerator  energies  leading  to 
errors  in  fluence  measurements  as  a function  of  energy. 

Figure  2-7  shows  the  dial  settings  vs.  beam  energy  corresponding  to  the 
solid  calibration  curve  in  Fig.  2-6.  In  addition,  the  switching-magnet  current 
required  to  bend  the  beam  cut  a side  port  at  45°  is  shown  (triangles) . The 
residual,  field  of  the  iron  poles  must  be  taken  into  account  for  this  curve 
fo  be  useful.  The. triangle  data  points  are  obtained  after  polarizing  the 
poles  by  running  a maximum  magnet  current  in  an  aiding  direction.  If  the  field 

o 

direction  is  now  reversed,  so  that  the  beam  will  exit  the  other  45  port, 
an  additional  current  (solid  line)  is  requireo  to  overcome  this  residual 
polarization  field.  If  a flip  coil  or  Hall  probe  were  inserted  in  the  beam 
can  so  that  the  magnetic  field  could  be  reproducibly  set,  the  switching  magnet 
could  be  used  as  a beam  energy  analyser.  High  energy  resolution  (accuracy  in 
determining  beam  energy)  would  demand  that  entrance  and  exit  slits  be  installed. 
Such  a procedure  is  highly  recommended  since  a constant  check  would  then  be 
* possible  on  the  generating  voltmeter.  A third  independent  energy  calibration 

check  can  be  obtained  by  building  a simple  Li  6-eutechtic  source  for  the 
entrance  slit.  Since  the  magnetic  analyser  passes  a constant  momentum  beam  for 
a fixed  magnetic  field,  and  since  the  mass  of  a Li  ion  is  very  approximately 
14,000  times  that  of  an  electron,  Li-ion  accelerating  voltages  of  less  than 
500  V will  cover  the  calibration  range  up  to  1 MeV  electrons.  When  comparing 
equal  momenta,  it  must  be  remembered  that  the  electrons  are  relativistic. 
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SWITCHING  MAGNET  CURRENT  (AMPS) 
FOCUS  CURRENT ( AMPS) 


Therefore,  the  magnetic  field  B required  to  pass  a constant  electron  momentum 


P or  Li  ion  momemtum  V. . will  be  given  by 


where  it  is  desired  to  determine  the  proportionality  constant  K.  The  two  momenta 


must  be  equal,  for  a fixed  field  B;  therefore 


where 


In  these  equations  V is  the  voltage  necessary  to  accelerate  a Li  ion  sufficiently 


such  that  it  will  pass  through  the  analyzer  whose  field  is  B,  m is  the  electron 


rest  mass,  is  the  Li  isotope  rest  mass,  c is  the  speed  of  light,  and  E is 

the  electron-beam  kinetic  energy  ( or  beam  energy)  sufficient  to  cause  an 


electron  beam  to  pass  through  the  analyser.  The  separation  of  the  two  Li  isotope 


lines  will  be  an  indication  of  energy  resolution 


In  order  to  be  able  to  make  accurate  fluence  measurements  (total  number  of 


electrons/cm  ),  beam  scanners,  beam  steerers,  and  beam-current  integrators  had 


to  be  constructed.  The  circuit  diagram  for  the  beam-current  integrator  is 


shown  in  Fig.  2-8  while  the  circuit  diagram  for  the  integrator  power  supply  is 


shown  in  Fig.  2-9.  The  integrator  is  basically  a voltage  integrator  which 


Integrates  the  analog  output  voltage  of  a Keithley  410  current  meter  up  to  a 


voltage  determined  by  the  comparator.  Once  this  voltage  is  reached,  the  compar 


ator  furnished  the  proper  pulses  to  discharge  the  integrator  feedback  capacitor 


and  also  pulses  the  electromechanical  counter  to  give  one  count.  The  process 


is  then  repeated.  The  accuracy  of  the  integrator  can  be  adjusted  to  better 
than  1 0.2%  ty  selecting  the  proper  current  at  pt.  B by  adjusting,  R^  so  as  to 
just  counterbalance  the  offset  voltage  of  the  1741G  op  amp.  A complete 
description  of  the  operation  and  adjustment  of  this  integrator  has  been 

g 

published. 

Beam  scanners  for  the  1 MeV  Van  de  Graaff  were  also  designed  and  constructed 
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under  this  project  and  have  been  published.  Figure  2-10  shows  the  basic 
circuit  to  provide  a triangular  wave  form  to  the  driver  power  amplifiers.  The 
basic  circuit  was  constructed  from  modules.  A Wavetek  triangular  wave  generator 
module  for  each  channel  was  used.  selects  the  frequency  of  oscillation 

while  a 5K  gain  control  and  1714G  op  amp  determine  the  amplitude.  Audio  power 
amplifiers  were  used  to  amplify  this  signal  and  power  the  deflection  coils 
(TV  yoke).  The  front  end  of  the  power  amplifier  (Universal  Tiger-175A)  was 
redesigned  so  that  (+)  or  (-)  dc  offset  could  be  provided  (Fig.  2-11)  to 
furnish  beam  steering.  Figure  2-12  demonstrates  the  linearity  of  the  sweep. 

This  figure  was  made  with  a 1 MeV  electron  beam  striking  an  A1  foil  coated 
with  CdS  powder.  Figure  2-13  demonstrates  the  excellent  beam-scan  uniformity. 
This  figure  shows  the  optical  density  produced  in  an  irradiated  microscope 
slide  which  is  proportional  to  the  number  of  electrons  Striking  the  glass. 
Uniformity  of  scan  is  estimated  to  be  better  than  ± 1%. 

In  addition  to  the  above  accelerator  modifications,  several  different 
types  of  experiments  have  been  set  up  at  the  end  of  the  beam  line  at  various 
times  during  this  contract.  Figure  2-14  shows  the  accelerator  vault  set  up 
for  a luminescence-band  defect-production  experiment.  Electrons  from  the 
Van  de  Graaff  are  accelerated  and  shot  down  the  beam  tube  to  irradiate  the 
sample  of  interest  which  is  mounted  at  the  tip  of  the  cold  finger  of  the 

21 


Optical  density  of  irradiated  microscope  slide 


cryostat  (Dewar).  The  cold  tip  makes  an  angle  of  45°  with  both  the  beam 
direction  and  the  optical  axis  of  the  spectrometer  and  lens  system.  Since 
the  electron  beam  penetrates  thin  samples  completely  the  electron-hole  pairs 
are  created  ddep  in  the  sample.  This  cathodoluminescence  therefore  has  the 
advantage  over  photoluminescence  in  not  being  surface  sensitive.  Figure  2-15 
is  a diagram  of  the  vaiable-temperature  irradiation  dewar.  The  A1  window 
on  the  LN2  shield  scatters  the  beam  slightly  for  uniformity.  In  later  exper- 
iments where  fluence  accuracy  was  a requirement,  the  scattering  foil  was 
moved  up  into  the  beam  tube.  An  insulated  beam-defining  aperature  was 
mounted  on  the  LN^  cold  shield  and  grounded.  The  rest  of  the  Dewar,  which 
is  insulated  from  ground,  then  served  as  a Faraday  cup.  Quartz  windows 
allowed  the  luminescence  to  be  collected  by  the  optical  system.  Electrical 
measurements  were  also  made  in  this  Dewar  by  making  use  of  the  feedthroughs 
at  the  top.  Thermometry  was  accomplished  using  carbon-resistor  and  Pt- 
resistance  themometers. 

For  IR  optical  measurements,  a double-pass  monochrometer  was  assembled  from 
several  junked  Perkin-Elmer  IR  systems  to  provide  an  operating  doublepass 
prism  monochrometer  which  could  cover  the  spectral  range  from  3000X  to  2.5um 
using  a RCA  Quanticon  photomultiplier  tube  and  a PbS  photoconductor  as  detectors. 
This  experimental  system  is  shown  in  Fig.  2-16.  The  irradiation  Dewar  is 
shown  at  the  top  of  this  figure.  Photoluminescence  could  be  compared  with 
cathodoluminescence  by  illuminating  the  sample  with  light  from  a high-pressure 
Hg  lamp. 

EPR  samples  were  irradiated  in  an  immersion  Dewar.  The  samples  were  mounted 
on  the  end  of  a stainless  steel  tube  with  clips  which  could  be  opened  from  the 
top.  The  tube  was  keyed  so  that  orientation  could  be  maintained.  A substitute 
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Diagram  of  iri'adiation  dewar. 


flux  collector  could  also  be  inserted  in  the  sample  position  so  that  the  ratio 
of  fluence  to  integrated  current  collected  on  the  Dewar  could  be.  measured. 

This  provided  a measure  of  the  beam  scattering  by  the  L.He.  Once  the  sample 
was  irradiated,  a small  transfer  Dewar  was  lowered  into  the  irradiation  Dewar 
and  a sample  transfer  made  under  L.He.  The  small  Dewar  could  then  be  trans- 
ferred to  the  EPR  Dewar  where  the  sample  was  finally  released  again  under 
L.He.  Although  several  samples  were  lost  in  developing  the  experimental 
technique,  this  system  proved  to  be  quite  reliable  eventually. 

Figure  2-17  shows  the  circuit  diagram  for  a present  fluence  accelerator 
shut-off.  The  homebuilt  integrator  mentioned  previously  was  eventually  replaced 
by  a commercial  unit  (B.I.C.  1000C)  once  the  usefulness  of  integrators  for 
fluence  measurements  had  been  demonstrated  by  use.  This  new  integrator  has  the 
advantage  of  better  accuracy  and  also  a preset  which  produces  a 5V  drop  when 
the  charge  collected  reaches  the  preset  value.  The  circuit  in  Fig.  2-17  used 
this  vdltage  level  change  to  shut  off  the  accelerator.  It  is  now  possible  to 
irradiate  samples  overnight  with  automatic  shut-off  when  the  desired  fluence 
is  reached. 

For  the  last  year,  considerable  effort  on  this  contract  has  been  devoted 
to  an  examination  of  radiation-damage  effects  in  solar  cells.  A reoccurring 
problem  is  to  evaluate  either  the  minority-carrier  diffusion  length  of  life- 
time. These  two  parameters  are  related  by  the  relation  L = >A)t  where  L 
is  thd  diffusion  length,  T is  the  lifetime,  and  D is  the  diffusion  coefficient 
(D  - p(kT/e),  by  the  Einstein  relation).  By  measurement  of  the  solar  cell 
short-circuit  current  denisty,  Jgc»  a determination  of  the  diffusion  length 
can  be  made  from  the  relationship 


where  J is  the  beam  flux(amp/cm  ) and  S is  the  specific  ionization  factor 
sc  o 

which  can  be  determined  either  experimentally  or  theoretically.  For  Si, 

= 225  e-h  pairs/u  at  1 MeV.  Figure  2-18  is  the  cirucit  used  to  make 
these  measurements.  Beam  current,  ID,  is  collected  in  the  Faraday  cup  of 

D 

aperature  A,,.  This  signal,  JL  = ID/A_,  becomes  the  flux  signal  which  is 

O D D D 

converted  to  the  voltage  +V^  by  the  Keithley  110  and  741  op  amp  of  gain  = 5/3. 
The  cell  short-circuit  current  Igc  is  fed  into  the  Keithley  419  and  741  op  amp 
of  gain  = 1.  The  full-scale  voltage  outputs  of  the  Keithleys  are  3V  and  5V. 
This  signal  becomes  -V  . The  ratio  is  formed  by  the  Analog  Devices  voltage 
divider  whose  output  is  fed  to  a recorder.  With  this  system,  continuous 
diffusion-length  degradation  (independent  of  beam  current  f luxuatior.s)  may 
be  observed.  Figure  2-19  shows  a typical  trace.  It  would  be  desirable  to 
use  a digital-to-analog  circuit  to  drive  the  x-axis  ot  the  recorder  using  the 
flux  counting  circuit  of  the  integrator.  This  has  not  been  built  yet. 

Fluence  measurements  are  presently  added  manually  to  the  graph. 

A second  independent  measurement  of  the  minority  carrier  lifetime  has 
also  been  made  or.  Si  cells  by  using  a pulse-bias  open-circuit  technique.  In 
this  measurement,  a solar  cell,  which  always  meets  the  requirement  that  the 


junction  is  abrupt  and  one  sided,  is  forward  biased  sufficiently  that  the 
depletion  width  is  removed.  On  removal  of  this  bias,  which  can  also  be  obtained 


by  a pulse  of  light,  the  depletion  width  formed  produces  a change  in  the  diode 
open-circuit  voltage  which  then  yields  information  about  the  minority-carrier 
lifetime  in  the  low-doping  side  of  the  junction  (the  side  of  interest  here). 

From  simple  p-n  junction  theory,  the  minority  carrier  hole  density  in  the 


n-region  is  given  by 


P « Pne 


qV/kT 
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Fig*  2-18.  Circuit  to  measure  solar-cell  carrier-diffusion  lengths 
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Since 


then 


p = + Ap, 

V = — ln(i  + Ap/p  ). 
q n 


Also,  the  excess  carrier  density,  Ap,  decays  on  the  removal  of  the  forward  bias 
by  the  equation 

Ap  =■  Apoe_t/T 

where  T is  the  minority  carrier  lifetime  and  Ap^  is  the  excess  minority  carrier 
density  at  t = 0 when  the  bias  is  removed.  If  is  the  value  of  the  open 
circuit  voltage  at  t = 0,  then 


and 


kT 

V = — ln(l  + Ap  /p  ) 
o q o n 


V(t)  = — ln[l  + feqV°/kT-l)e  t/T], 

q 


For  t«T  and  VQ>>kT/q,  the  last  equation  yields 


dV  = kT  1 
dt  q T 

or 

q 1 

T * " kT  (dV/dt)  . 

Therefore,  a measurement  of  the  slope  of  the  decay  of  the  open  -circuit  voltage 
vs.  time  yields  the  lifetime  directly. 

Ke  have  used  both  of  the  above  techniques  on  a variety  of  Si  solar  cells 
and  have  obtained  good  agreement  between  the  two  techniques. 

In  addition  to  the  Van  de  Graaff  facility  work  and  the  design  and  construction 
of  experimental  equipment  already  described,  considerable  effort  has  been 
spent  on  the  rest  of  the  accelerator  facility.  New  and  safer  accelerator  vaults 
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for  all  the  machines  were  planned  under  this  contract  and  have  now  been  installed. 
The  400  keV  Van  de  Graaff  was  moved  from  an  unshielded  open  position  to  a new 
vault.  Help  was  furnished  from  this  contract  in  converting  this  accelerator 
to  protons.  This  machine  now  resides  in  approximately  the  same  place  as  the 
old  Cockroft-Walton  which  we  also  helped  remove. 

Considerable  effort  was  also  expended  in  setting  up  the  ion-implantation 
accelerator,  aligning  the  beam,  repairing  the  analyzing-magnet-winding  shorts, 
etc.,  as  well  as  training  personnel  to  operate  the  machine. 

This  concludes  the  experimental  equipment  section  of  this  report  relating 
to  radiation  effects  experiments.  The  following  sections  describe  results  of 
experimental  studies  of  radiation  damage  in  compound  semiconductors. 

2.3  Radiation  Effects  in  Compound  Semiconductors 

In  this  section,  brief  summaries  of  experimental  work  concerned 
with  defect  identification,  defect  production  mechanisms*  and  irradiation 
effects  in  semiconductor  devices  will  be  described.  The  work  will  be  broken 
down  into  the  following  areas: 

2.3.1  Radiation  Effects  in  II-VI  Compounds 

2.3.2  Defect  Semiconductor  Radiation  Damage  - In2Te^ 

2.3.3  Annealing  Kinetics  and  Radiation  Annealing  in  Si,  Ge,  and  GaAs 

2.3.4  Solar  Cell  Radiation  Effects 

2.3.5  Proton-Induced  Characteristic  X-ray  Analysis. 

Since  most  of  this  work  has  been  published,  only  brief  summaries  of 
the  major  conclusions  will  be  presented  here.  Since  the  development  of 
scientific  models  contains  an  evolutionary  element,  attempts  will  be  made  to 
summarize  our  present  understanding  of  these  problems  in  terms  of  the  prospective 
we  have  been  able  to  develop  during  the  period  of  this  contract.  Attempts 
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will  wise  be  made  Co  suggest  the  areas  in  which  understanding  is  still  lacking 


1 


*r 


‘ 
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and  to  indicate  fruitful  avenues  of  future  research. 

2.3.1  Radiation  Effects  In  II-VI  Compounds 

In  the  early  70' s,  at  the  beginning  of  this  contract,  very 
iittie  was  understood  concerning  the  nature  of  the  defects  produced  by 
irradiacion  in  II-VI  compounds  and  how  to  identify  them  by  their  electrical 
and  optical  properties,  etc.  The  displacement  thresholds  had  been  measured 
for  a number  of  the  II-VI' s and  some  luminescence  bands  produced;  however, 
the  correlation  of  these  bands  with  particular  defects  was  sketchy  in  most 
cases,  the  results  being  inferred  only  from  one  kind  of  measurement. 

Radiation  damage  had  been  studied  extensively  in  Si  and  Oe  and  it  had 
been  discovered  that  vacancies  in  these  materials  migrate  at  temperatures 
near  L.i^.  The  interstitials  were  thought  to  migrate  near  or  below  L.He 
temperature.  Since  most  of  the  semiconductor-radiation-effects  scientists 
had  cut  their  teeth  on  these  materials,  there  was  a general  belief  in  this 
community  that  the  rest  of  the  semiconductors  should  behave  like  Si  and  Ge. 

The  rather  obvious  point  that  defects  in  a compound  semiconductor  should  move 
on  their  own  sublattices  (requiring  considerably  more  thermal  energy  than 
in  the  elemental  semiconductors)  was  missed  bv  most.  The  literature  of  the 
60' s Is  full  of  speculation  suggesting  that  new  low- temperature  annealing  stages 
remained  to  be  found  in  the  II-VI  semiconductors.  Very  few  took  the  trouble 
to  look,  however. 

Some  controversy  existed  concerning  the  displacement  thresholds  also. 
Although  many  of  the  II-VI  and  III-V  displacement  energies  had  been  found  to 
be  in  the  5-10  eV  range,  it  was  felt  that  these  numbers  were  too  low  compared 
with  the  values  of  15-30  eV  from  C,  Si,  and  Ge.  To  complicate  matters,  estimates 
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of  the  displacement  thresholds  in  ZnO  had  placed  the  value  for  the  displacement 
energy  of  Zn  in  ZnO  above  600  keV.^  When  later  EPR  experiments^  showed  that 
the  defect  produced  by  these  very  high  energy  electron  Irradiations  was  the 
oxygen  vacancy,  and  not  a Zn  defect,  the  conclusion  was  drawn  that  should 
be  larger  for  ZnO  than  diamond,  a difficult  point  to  defend  for  long.  Further- 
more, in  CdS,  Kulp^  had  identified  the  Cd  displacement  threshold  ar.  290  keV  by 
the  production  of  a luminescence  band  at  1.03um.  Bryant  and  Cox  had  reported 

12 

forming  this  band  at  125  keV  suggesting  that  sulfur  displacements  were  involved. 
The  situation  was  far  from  clear  concerning  the  displacement  energies.  A 
knowledge  of  these  energies  is  essential,  however,  to  an  understanding  of  the 
association  of  various  annealing  stages  with  defects  on  particular  sublattices 
and  to  physical  property  changes  associated  with  defects  on  the  two  sub- 
lattices. 

One  other  point  should  be  mentioned.  Most  of  the  II-VI  compounds  are 
usually  n-type  and  strongly  resist  attempts  to  dope  them  p-type.  It  has  been 
observed  in  Ge  and  Si  that  defect  stability  is  affected  by  the  defect  charge 
state,  thereby  producing  radiation  effects  which  are  quite  different  in  n-  and 
p type  samples  of  the  same  materials.  Exploring  these  differences  in  the 
II-VI  system  is  impossible.  It  is  possible,  however,  to  gain  some  understanding 


of  the  charge-state  effects  by  studying  ZnTe,  the  only  II-VI  which  is  strongly 
p-type.  Furthermore,  at  the  beginning  of  this  work  effort,  no  radiation  effects 
had  been  performed  on  this  substance.  We  therefore  will  begin  this  section  with 
a discussion  of  the  work  done  under  this  contract  on  ZnTe.  This  will  be 
followed  by  a discussion  of  radiation  effects  studies  on  CdS,  ZnO  and  ZnS. 

2 . 3 . 1 . 1 Radiation  Damage  in  ZnTe 

Since  ZnTe  had  not  been  previously  investigated,  it  was 
important  to  determine  the  displacement  thresholds.  This  was  done  by  irradiating 
single-crystal  samples,  which  were  cooled  to  near  L.Pe  temperature,  with  electrons 
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of  various  energies  and  measuring  the  effect  of  this  irradiation  by  monitoring 

13 

changes  in  the  ZnTe  cathodoluminescence . The  near-band-edge  cathodoluminescent 
intensity  was  observed  to  increase  with  increasing  beam  energy  as  shown  in 
Fig.  2-20.  Figure  2-21  shows  a typical  low-resolution  cathodoluminescent 

1 8 

spectrum  before  and  after  249  keV  electron  irradiation  to  a fluence  of  1.2  x 10 

, 2 
e/cm  . 

A displacement  threshold  can  be  derived  from  the  data  in  Fig.  2-20  by 
using  the  relation 

5?“  «<E>  « <E  - Ed)2 

where  dl/d0  is  the  initial  slope  in  Fig.  2-20  at  each  irradiation-electron-beam 

energy  F,,  n(E)  is  the  average  number  of  displaced  atoms  per  electron,  and 

is  the  electron-beam  energy  corresponding  to  the  threshold.  From  the  above 

1/2 

relation,  a plot  of  (dl/d0)  vs  E yields  a straight  line  as  shown  in  Fig.  2-22. 

We  take  the  intercept,  E = 185  i 5 keV,  as  the  displacement  threshold.  This 

d 

electron-beam  energy  corresponds  to  an  atomic  lattice-atom  recoil  energy  of 
7.35  eV  for  a Zn  displacement,  or  3.75  eV  for  Te  displacements.  We  can  identify 
the  atomic  species  displaced  as  Zn  for  the  reasons  outlined  below. 

(1)  The  mass  of  the  Te  atom  it  approximately  twice  that  of  the  Zn  atom. 

We  would  expect,  therefore,  that  the  lowest  threshold  observed  should  belong 
to  the  Zn  atom. 

(2)  Recoil  energies  for  anions  in  the  other  II-VI  compounds  fall  between 

14 

7.8  and  8.7  eV.  If  the  observed  threshold  corresponded  to  Te  displacements, 

then  the  Te  recoil  energy  would  be  approximately  a factor  of  2 too  low  as 
compared  to  other  II-VI  compounds. 

(3)  The  energy  required  to  displace  an  atom  in  a lattice  should  be  related 
to  the  lattice  bond  energy  since  this  energy  is  the  dominant  term  in  any 
calculation  of  displacement  thresholds  from  first  principles.  Figure  2-23  shows 
a graph  of  estimated  bond  energy  vs.  cation  recoil  energy,  T^,  for  a number  of 
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Fig.  2-22.  Square  root  of  production  rate  vs.  electron  energy 
of  2.320  eV  transition. 
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II-VI  compounds.  Our  value  of  7.35  eV  for  Zn  displacements  fits  this  trend 
well . 

lie  therefore  conclude  that  the  threshold  observed  corresponds  to  the 

displacement  threshold  for  Zn  in  ZnTe. 

Further  evidence  is  obtained  for  this  model  by  an  interpretation  of  the 

optical  transitions  involved  in  the  production  of  luminescence  near  the  band 

edge,  figure  2-24  shows  high-resolution  photoluminescence  spectra  of  irradiated 

ZnTe  taken  at  4.2°K  and  77°K.  This  edge  emission  can  be  understood  on  the 

basis  of  the  energy  level  model  for  ZnTe  which  we  were  the  first  to  propose 

in  1971. ^ The  energy  level  diagram  is  shown  in  Fig.  2-25. 

As-grown  ZnTe  contains  acceptor  levels  at  0.5  eV  and  0.15  eV  above  the 
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valence  band,  as  determined  by  electrical  measurements.  Two  edge-emission 
series  with  zero  phonon  lines  at  2.320  and  2.234  eV  are  observed  in  Fig.  2-24. 
Both  series  increase  in  intensity  for  electron  irradiations  above  the  Zn 
displacement  threshold,  as  shown  in  Fig.  2-21.  These  two  transitions  are 
probably  the  free -to-bound  transitions  E.  and  shown  in  Fig.  2-25  to  the 
two  shallow  acceptor  levels  seen  in  as-grown  ZnTe  by  electrical  measurements. 
(These  transitions  could  also  be  from  very  shallow  donor  levels  near  the  conduc- 
tion band).  The  corresp ondiag  optical  absorption  transitions  and  A^  have 
been  seen  in  as -grown  samples  by  Nahory  and  Fan.^  They  have  also  observed  the 

A^  and  optical  transitions  from,  the  F_,  split-off  valence  band  while  we 
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have  observed  the  corresponding  cachodoluminescence  transitions.  These 
transitions  demonstrate  that  we  are  dealing  with  acceptor  levels,  not  donors. 

From  the  radiation-damage  evidence,  coupled  with  the  fact  that  these  defects 
are  produced  at  4.2°K  but  are  stable  at  300°K,  the  only  possible  defects  which 
could  have  these  properties  are  isolated  Zn  vacancies.  (Interstitials  have  now 
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Fig.  2-24.  High-re solution  photoluainescence  oi  Zrtle  edge  emission  at  4.2°K  and  77°K. 


been  shown  to  migrate  above  77  K but  below  room  temperature  in  most  of  the  II-VI 
compounds) . 

We  can  now  make  the  following  important  conclusions  about  ZnTe: 

(1)  ZnTe  is  always  p-type  due  to  an  excess  of  Zn-vacancy  shallow  acceptor 
levels  (double  acceptor  with  both  levels  near  the  valance  band)  which  exist 'in 
as-grown  ZnTe. 

(2)  Electron  irradiation  above  the  Zn  displacement  threshold  of  185  keV 

(T,  = 7.35  eV)  increases  the  concentration  of  these  Zn  vacancies  and  hence  the 
a 

room-temperature  hole  conductivity  (as  we  have  observed) . 

(3)  The  green  edge-emission,  which  dominates  the  luminescence  spectra  0*" 
ZnTe  at  4.2°K,  is  due  to  the  presence  of  Zn  vacancies  in  the  lattice  after 
crystal  growth. 

(4)  The  zinc  vacancy  in  ZnTe  is  stable  at  room  temperature  in  contrast 
to  the  behavior  of  vacancies  in  the  elemental  semiconductors  which  generally 
migrate  below  room  temperature. 

Subsequent  studies  on  ion-implanted  ZnTe  have  modified  the  above  model  in 
detail  but  not  in  substance. ^ ^ The  4.2°K  photoluminescence  of  ZnTe  was 
monitored  as  a function  of  isochronal  annealing  in  an  argon  atmosphere.  Typical 
spectra  before  and  after  300°C  ischronal  annealing  is  shown  ir.  Fig.  2-26.  The 
relative  intensities  of  the  luminescence  transitions  as  a function  of  anneal 
temperature  is  shown  in  Fig.  2-27.  The  temperature  dependence  of  these  trans- 
itions is  shown  in  Fig.  2-28.  All  as-grown  samples  exhibit  a strong  exciton 
line  at  532oX(I^)  and  two  edge  emission  series,  the  first  at  5318&(A^)  and 
the  second  at  5543X(B).  The  A^  and  B transitions  are  related  to  the  two 
transitions  described  above. 

After  300°C  annealing,  a marked  increase  is  noted  in  the  intensity  of 
1^  and  the  B series  transitions.  In  addition,  all  samples  exhibit  a new 
A-series  line  (Aq)  at  5300X.  The  temperature  dependence  of  the  A-series 
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Fig.  2-26.  Edge  emission  showing  effects  of  300°C  argon  anneal. 
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Fig.  2-28.  Temperature  dependence  of  A transitions. 
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transitions,  and  a possible  exciton  line  at  5218$,  indicate  that  two  shallow 

acceptors  are  involved  (corresponding  to  the  0.05  eV  level  mentioned  previously) 

with  ionization  energies  of  0.052  eV  and  0.060  eV.  These  values  correspond 

to  those  found  by  electrical  measurements  which  have  been  attribut'd  to  the 

first  ionization  level  of  the  zinc  vacancy  and  to  a Group  I substitutional 

acceptor.  This  close  proximity  of  vacancy  shallow  levels  and  substitutional 

impurity  levels  of  the  same  electrical  type  (donor  and  acceptor),  will  be 

shown  to  be  a reoccurring  theme  in  the  nature  of  Isolated  lattice  defects  in 

II-VI  compounds.  This  proximity  has  led  to  considerable  confusion  in  the 

II-VI  literature  when  dealing  with  the  separate  identification  of  isolated 

lattice  defects  and  substitional  impurities. 

ZnTe  samples,  ion  implanted  with  90  keV  Cd,  Ar,  and  Xe  ions  at  room 
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temperature  with  fluences  from  10  to  10  ions/cm  , have  also  been  investi- 
gated. Photoluminescence  was  measured  at  4.2°K  as  a function  of  10  min. 

isochronal  anneals  in  flowing  Ar.  None  of  the  samples  show  luminescence  after 
12  2 

a dose  of  3 x 10  ion s/cm  . Complex  recovery  of  the  edge  emission  is  observed 

between  300-500°C.  A broad  band  at  5900X  was  produced  at  the  low  end  of  this 

temperature  range  but  disappeared  completely  after  annealing  to  500°C.  Samples 
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implanted  with  3 x 10  to  4 x 10  Cd  /cm  exhibit  exciton  emission  at 

5293X  after  300-400°C  anneals;  this  was  not  found  in  any  of  the  cleaved  samples 
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or  in  those  samples  implanted  with  Ar  or  Xe.  Samples  implanted  with  Cd 
suggest  another  exciton  line  at  5287$  which  appears  to  be  an  isotope  shift  of 
unexplained  nature. 

The  appearance  of  the  Cd  exciton  lines,  the  broad  radiation-damage  band 
at  5900X,  and  the  two  shallow  acceptor  levels  in  the  0.05  eV  range,  which  are 
enhanced  as  a result  of  annealing  in  an  inert  atmosphere,  all  suggest  that 


considerably  more  information  concerning  the  nature  of  acceptor  levels  in 
II-VI  compounds  can  be  obtained  by  further  study  of  ZnTe. 

2 . 3 . 1 • 2 Radiation  Damage  in  CdS 

A considerable  effort  on  this  contract  has  been  devoted 
to  understanding  the  nature  of  radiation  produced  defects  in  CdS.  The  reasons 
for  this  effort  are  twofold.  First,  CdS  is  a useful  device  material  both 
as  a thin-film  photodetector  in  the  visible  and  also  as  a thin-film  solar-cell 
material.  Second,  the  understanding  of  the  materials  properties  of  CdS  is 
far  advanced  over  that  of  most  of  the  II-VI' s.  It  therefore  holds  an  important 
place  in  obtaining  a general  understanding  of  II-VI  semiconductors. 


As  mentioned  in  a previous  section,  a discrepancy  existed  between  the 

displacement  threshold  measurements  of  Kulp  et.  al.^  and  those  of  Bryant  and 
12 

Cox.  Both  groups  had  observed  the  production  of  a luminescent  band  at 
1.03pm  by  electron  irradiation  but  Kulp  claimed  that  the  threshold  for  this 
production  was  at  290  keV  and  represented  Cd  displacements  while  Bryant  and 
Cox  supported  the  view  that  the  threshold  for  the  production  of  the  1.03pm  band 
was  125  keV  and  represented  S displacements. 

These  experiments  cast  considerable  doubt  concerning  the  measurement  of 
the  energy  required  for  cadmium  displacements  in  CdS.  For  these  reasons 
additional  experiments  were  undertaken. ^ These  measurements  were  made  near 
L.Ke  temperature  in  contrast  to  Kulps  work,  which  was  near  L.^  temperature. 

Typical  before-  and  after-irradiation  cathodoluminescence  spectra,  taken  at 
5°K,  are  shown  in  Fig.  2-29.  Figure  2-29(A)  is  the  spectrum  before  irradiation. 

Fig.  2-29(B)  the  spectrum  after  400  keV  irradiation  to  a fluence  of  1.89  x 10^e/cm", 
and  Fig.  2-29  (C),  after  a room-temperature  anneal.  The  1.03pm  band  is  apparent 
after  irradiation  but  shows  considerable  annealing  at  room  temperature.  The 
series  of  bands  in  the  1.6-2. 2pm  region  have  been  attributed  to  the  internal 


transitions  between  levels  of  substitutional  Cu 
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function  of 


Figure  2-30  shows  the  growth  of  the  1.03pm  band  during  irradiation  for 

various  beam  energies.  The  displacement  thresholds  can  be  determined  as  in 

1/2 

the  previous  section  by  a plot  of  (dl/d0)  vs.  beam  energy,  where  dl/d0 

is  the  initial  slope  of  each  curve  in  Fig.  2-30.  This  plot  is  shown  in 

Fig.  2-31.  It  is  apparent  that  two  displacement  thresholds  are  evident  at 

125  keV  (S-displacements)  and  at  314  ± 5 keV  (Cd-displacements) . The  lower 

12 

value  agrees  with  the  result  of  Bryant  and  Cox  while  the  upper  value  is 

close  to  the  Cd  threshold  of  Kulp.'*  Although  the  reasons  for  the  difference 

7 12 

between  the  results  of  Kulp  and  Bryant  and  Cox  are  not  clear  from  the  data 
presented  here,  it  is  evident  from  Fig.  2-31  that  both  these  experiments  are 
substantially  correct.  Our  results  obviously  represent  an  intermediate 
situation  in  that  both  thresholds  can  be  seen  in  our  experiment.^ 

The  differences  in  these  various  experiments  are  the  result  of  stoichiometry 
differences  among  the  various  starting  crystals.  If  the  1.03um  luminescence 
is  the  result  of  a transition  between  a sulfur-vacancy  level  and  a cddmium- 
vacancy  level,  then  a large  concentration  of  cadmium  vacancies  prior  to  irrad- 
iation would  emphasize  the  S threshold.  Only  the  Cd  threshold  would  be  seen, 
however,  in  a crystal  containing  a large  concentration  of  sulfur  vacancies  but 
no  cadmium  vacancies.  This  modal  is  substantiated  by  the  fact  that  our  crystal 

had  a slight  excess  of  S while  those  of  Bryant  and  Cox  contained  a lirge  excess 
12 

of  S.  Kulp's  crystals  apparently  contained  an  excess  of  Cd,  as  evidenced 
from  the  fact  that  he  made  electrical  measurements  on  some  of  these  samples 
easily  at  low  temperaure.  We  emphasize  again  that  this  model  suggests  the  basic 
stability  of  vacancies  in  CdS  at  room  temperature,  a result  similar  to  the 
conclusions  for  ZnTe. 

The  displacement  threshold  observed  here  for  Cd  at  L.He  temperature  is 


FLUENCE  (e/cm  ) 

Production  of  the  1.03pm  band  in  CdS  as  a function  of  fluence 
for  various  electron-irradiation  energies. 
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Fig.  2-31.  Square  root  of  the  production  rate  vs.  electron 
energy  for  the  1.03Mm  band  in  CdS. 


be  noted  that  both  Kulp's  old  value  at  77LK  and  our  newer  one  fit  the  lattice 
bond  model  well.  The  77°K  threshold  is  lower,  however,  than  that  at  10°K 
suggesting  a temperature-dependent  displacement  threshold  as  has  been  observed 
in  other  semiconductors. 

The  above  model  for  the  production  of  1.03pm  luminescence  presupposes  that 

cadmium  vacancies  are  deep  acceptors  (this  has  been  demonstrated  for  Zn 
23 

vacancies  in  ZnSe  ) and  that  the  sulfur  vacancies  are  shallow  donors.  Further- 

evidence  for  the  donor  nature  of  the  sulfur  vacancy,  its  migration  energy, 

thermal  ionization  energy,  etc.,  will  be  presented  from  low-temperature 

luminescence  experiments  and  electrical  measurements.  A defect  misidentif ication 

by  Kulp  and  Kelley'*  for  the  7200$  band  will  also  be  demonstrated  as  well  as  a 

24 

second  misldentification  by  Kikuchi. 

Figure  2-32  shows  typical  emission  spectra,  before  irradiation,  near  the 
CdS  band  edge  for  undoped  (IT),  Cd-doped  (Cd),  and  S-doped  (S)  CdS  platelets, 
taken  at  7°K.  The  1^,  and  1^  lines  are  excitons  bound  to  neutral  acceptors, 
neutral  donors,  and  ionized  donors  respectively.  The  broad  emission  series  at 
5140&  is  known  as  the  edge  emission  (sie  ZnTe  section)  and  is  due  to  free-to- 
bound  or  bound-to-bound  transitions  in  general.  Eoth  photo-  and  cathodo- 
luminescence  have  been  used  to  monitor  radiation  damage  in  CdS  at  10°K  for 
electron-beam  energies  capable  of  producing  sulfur  displacements  only 
(125  $ E H 300  keV),  or  sulfur  plus  cadmium  displacements  (K  - 300  keV) . 

Plots  of  the  10°K  bound-exciton  luminescence  in  general  show  that  the 
intensity  decays  according  to  the  relation  I/I0  or  0 ^ for  large  fluence,  0 
(see  Fig.  2-33).  This  decrease  is  probably  due  either  to  a disturbance  of 
the  surface  states  by  irradiation  or  an  increase  in  optical-absorption  band  tailing. 
Also  seen  in  Fig.  2-33  is  an  increase  in  the  intensity  of  a new  bound-exciton 
line,  1R  at  4867$,  for  250  keV  irradiations  at  < 10°K  (between  the  thresholds), 


57 


°'-C>°~<250K^> 


'K  500 KeV 
^350  KeV 

700  KeV 
800  KeV 


500  KeV 


4869  A (I5  line) 

A B-l  1-5-69 (Cd) 
O B-I0-3I-69(U) 


350  KeV 


700  KeV 

Rline) 


FLUENCE  <j>  (e  4 m2  ^ 


10 

FLUENCE  <f>  (e  cm2) 

Fig.  2-33.  Decay  of  bound-exci ton  line  intensity  with  irradiation  in  CdS. 
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in  the  photoluminescent  spectrum.  It  was  felt  that  the  decrease  in  intensity 
of  all  the  bound  exciton  lines  for  photo  excitation,  which  is  shallow,  was 
probably  masking  the  production  of  the  I line.  This  was  confirmed  by 
monitoring  the  intensity  of  the  adjacent  1^  ionized-douor  line  with  both 
electron-beam  and  photo  excitation  as  shown  in  Fig.  2-34.  It  is  apparent 
that  the  electron  beam  samples  deeper  into  the  crystal  and  produces  less 
pronounced  intensity  decay. 

The  production  of  the  I bound  exciton  was  further  investigated  with 

K 

cathcdoluminsscence,  as  shown  in  Fig.  2-35.  There  is  a remarkable  energy- 

16  2 

dependent  growth  in  intensity  of  i_  up  to  fluences  of  10  e/cm  . Using  the 

K 

same  analysis  as  previously,  a threshold  for  the  production  of  thii  line  is 
observed  in  the  neighborhood  of  115  keV.  It  is  apparent  that  the  IR  line  is 
a result  of  sulfur  displacements;  all  the  other  bound-exciton  lines  adjacent 
to  the  I line  have  decreased  during  the  same  fluence>as  represented  by  the 

iv 

1^  line  intensity  in  Fig.  2-35. 

The  edge  emission  beyond  514oX  is  related  intimately  to  the  bound-exciton 
emission,  as  shown  in  Fig.  2-36.  This  figure  illustrates  the  two  types  of 
edge  emission  series  found  in  CoS  which  are  labeled  free-to-bound  and  bound-to- 
bound.  The  concentration  cf  shallow  donors,  as  well  as  the  intensity  of  bound- 
exciton  donor  lines,  is  related  to  the  intensity  of  the  edge  emission.  Figure  2-37 
shows  the  growth,  of  the  edge  emission  by  electron  irradiation  above  the  sulfur- 
displacement  threshold.  We  have  been  able  to  demonstrate  that  the  I line  is 

related  to  shallow  donors  which  are  increased  by  irradiation  above  the  sulfur- 

2 

displacement  threshold.'  This  production  of  shallow  donors  also  causes  an 
increase  In  the  edge  emission.  Since  these  irradiations  are  performed  near 
L.Ue  temperatuie,  and  for  energies  above  the  displacement  threshold,  the  sulfur 
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Fig.  2-35.  Production  of  IR  bound  exciton  line  in  CdS. 


CONDUCTION  BAND 


Fig.  2-36.  Energy-level  model  for  two  edge-emission  series  in  CdS. 


vacancy  appears  to  be  the  best  candidate  for  the  shallow  donor  level  (the 

S-lnterstitial  should  be  neutral  or  an  acceptor).  An  analysis  of  the  exclton 

spectrum  suggests  that  the  donor  level  Is  about  0.03  eV  below  the  conduction 
25 

band.  It  has  further  been  observed  by  optical  measurements  that  the  production 
of  the  I line  Is  possible  at  room  temperature.  Zeeman -effect  studies  on  this 
line  indicate  that  it  is  indeed  a donor  level  line. 

25 

Several  annealing  stages  have  been  studied  in  irradiated  CdS.  Figure  2-35 

indicates  a recovery  stage  for  the  intensity  of  the  bound-exciton  line  intensity 

at  180°K  which  can  be  shown  to  be  related  to  the  migration  of  sulfur  sublattice 

defects  (irradiation  energy  was  less  than  Cd-threshold) . In  addition,  a 

luminescence  band  at  7200X  is  observed  to  be  formed  during  this  stage  as  shown 

in  Fig.  2-39(a).  This  7200X  band  is  again  related  to  sulfur  sublattice  defects 

since  it  can  be  produced  below  the  Cd  threshold.  The  defect  responsible  for 

this  emission  appears  to  be  a sulfur  sublattice  defect  complexed  in  some  way 

25 

with  impurities  or  other  sites  near  the  surface  of  the  sample.  Its  intensity 
is  high  for  photoluminescence  but  very  low  for  cathodolvminescence.  Further- 
more, it  can  be  broken  up  by  irradiation  at  10°K  and  20C  keV  as  shown  in 
Fig.  2-39(b).  Kulp  has  attributed  the  7200X  band  to  isolated  sulfur  vacancies.^ 
This  identification  clearly  is  in  error  since  the  band  cannot  be  produced  at 
4.2°K  and  since  it  is  found  as  a result  of  annealing  (defect  migration). 

Kikuchi  has  suggested,  on  the  basis  of  thermal  neutron  irradiation  at  room 

temperature,  that  the  7200X  band  is  the  result  of  Cd  vacancies,  as  shown  in 
24 

Fig.  2-40.  This  clearly  is  also  incorrect  since  this  band  is  produced  at 

energies  below  the  Cd-displacement  threshold,  as  seen  in  Fig.  2-39(b).  The 

best  evidence  to  date  suggests  the  sulfur  interstitial  in  a complex  is  the 

25 

cause  of  this  defect  band. 
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(C)  FLUENCE  * (e/cm2)x I0"16 

Production  of  the  7200X  luminescence  band  (a)  by  isochronal  annealing 
and  (b)  as  a function  of  fluenr.e  before  a 200°K  anneal.  Also,  the 
radiation  annealing  of  this  band  is  displayed  in  (c). 


INITIAL 


Thermal  neutron  production  of  720oX  band  in  CdS  at  room  temperature 


Another  annealing  stage,  at  230  K,  has  been  observed  for  platelets 

irradiated  above  the  Cd  threshold,  as  shown  in  Fig.  2-41.  This  stage  and 

the  180°K  stage  are  seen  by  observing  the  recovery  of  the  1.03bm  bands. 

Minor  annealing  stages  are  also  observed  below  150°K;  these  have  also  been 
26 

reported  by  Vook  for  the  recovery  of  radiation-induced  thermal  resistance. 

This  thermal  resistance  data  is  seen  in  Fig.  2-42 (c).  Vook  has  observed  an 

increase  in  anisotropic  thermal  conductivity  produced  by  annealing  from  60  to 

150°K,  followed  by  a loss  of  this  anisotropy  in  a stage  at  180°K,  followed 

by  a nearly  complete  recovery  in  a broad  stage  extending  from  210-270°K  after 

2 MeV  irradiation.  He  attributes  the  formation  of  this  anisotropic  thermal 

conductivity  to  a re-orientation  of  defects  along  the  c-axis.  If  we  assume 

that  these  defects  are  Frenkel  pairs,  then  the  annihilation  of  these  pairs 

at  180°K  is  consistent  with  the  loss  of  the  anisotropic  thermal  conductivity, 

27 

with  the  larger  percentage  recovery  observed  with  electrical  measurements, 

shown  in  Figs.  2-42(a)  and  2-42(b),  and  with  the  large  thermal  conductivity 

28—30 

recovery  in  this  same  annealing  stage.  Also,  since  from  F,PR  results  it 

is  unlikely  that  vacancies  migrate  in  CdS  below  room  temperature,  we  attribute 

the  180°K  annealing  stage  to  the  migration  of  sulfur  interstitials,  as  is 

consistent  with  the  displacement  energies.  The  annealing  stage  at  230°K  is 

then  most  likely  to  be  the  onset  of  Cd-interstitial  migration  since  the  cadmium 

28 

vaca-'cy  is  known  to  be  stable  at  room  temperature  from  EPR. 

An  additional  annealing  stage  has  been  observed  by  us  in  room-temperature 
31  32 

tron  irradiated  CdS.  ’ Information  concerning  the  position  of  the  defect 
-n»r*'  levels  from  electrical  measurements  is  also  in  agreement  with  the 
.urements  above.  In  the  experiment  to  be  described,  250  keV 
i’  i • 'between  the  Cd  and  S displacement  thresholds)  and  1 MeV  irradiations 


Fig.  2-41.  Isochronal  annealing  of  the  1.03ym  and  1.65ym  bands  for 
samples  irradiated  above  the  Cd-threshold. 
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(a)  Isochronal  annealing  in  high-conductivity  CdS  crystals  alter 
irradiation  by  10  MeV  electrons  near  LN2  temperature;  (b)  isochronal 
annealing  in  sample  106-CG  ((a)  and  (b)  from  Ref.  27);  (c.)  fraction 
of  the  additive  thermal  resistivity  at  12°K  that  remains  after 
successive  15  min.  anneals,  (from  Ref.  26). 
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(above  both  cnresholds)  are  compared.  Figure  2-43  shows  that  both  energy 
electrons  produce  an  increase  in  carrier  concentration,  suggesting  the  production 
of  a shallow  donor  level.  The  ratio  of  defect  production  rates  at  the  two 
energies  is  very  close  to  what  one  would  expect  for  sulfur  displacements, 
suggesting  that  a sulfur-sublattice  defect  is  dominating  the  carrier  concentration. 

Figures  2-44  and  2-45  show  the  carrier  concentration  and  mobility  of 
these  CdS  samples  after  1 MeV  electron  irradiation  as  a function  of  temperature. 

It  appears  from  Fig.  2-44  that  two  donor  levels  exist  in  CdS  before  irradiation 
and  that  the  concentration  of  these  donors  increases  after  1 MeV  irradiation. 

The  room-temperature  mobility  is  dominated  by  lattice  scattering  and  does  not 
change  as  a result  of  1 MeV  irradiation  as  shown  in  Fig.  2-45.  Analysis  of 
these  electrical  data  leads  to  the  energy  level  scheme  shown  in  Fig.  2-46 (D) 
where  the  defect  donor  levels  at  0.022  eV  and  0.13  eV  are  shown  to  dramatically 
increase  after  irradiation  (see  data  in  Fig.  2-46 (E)).  The  concentration  of 
acceptors  does  not  change  appreciably. 

Isochronal  anneals  are  shown  in  Fig.  2-46(A),  (B),  and  (C) . Figure  2-46(A) 
is  for  1 MeV  irradiation,  (B)  for  250  keV  irradiation  (below  the  Cd  threshold 
but  above  the  S threshold),  and  (C)  for  an  unirradiated  sample.  In  the  unirradiated 
sample  a reverse  annealing  SLage  from  270-350°C  is  obtained  which  produces 
additional  donors.  We  attribute  this  to  the  evaporation  of  sulfur  from  the 
crystal.  The  dominant  annealing  stage  at  250°C  is  the  loss  of  irradiation- 
induced  donor  levels.  The  I donor-bound-exciton  line  also  disappears  in  this 
same  temperature  range.  It  follows  from  the  previous  low-temperature  annealing 
stage  analysis,  the  two-donor  character,  the  fact  that  this  defect  is  stable 
at  room  temperature,  and  the  fact  that  this  defect  can  be  produced  for  electron 
irradiation  energies  just  above  the  S displacement  threshold,  that  the  only 
possiblity  for  the  identification  of  this  defect  is  the  sulfur  vacancy. 
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Carrier  concentration  vs.  temperature  for  samples  which  are 
unirradiated,  irradiated  at  1 MeV,  and  annealed  at  260°C. 
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Temperature  dependence  of  mobility  for  same  conditions  as 
specified  in  Fig.  2-44. 
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Isothermal  anneals  of  the  250  C annealing  stage  add  support  to  this 
3 2 

argument.  Figure  2-47  shows  these  isothermal  anneals  at  250°C,  275  C and 

300°C.  An  analysis  of  these  data  indicates  that  the  annealing  kinetics 
32 

are  second  order.  In  Fig.  2-48  the  data  are  fitted  according  to  the 
equation 

ni  = ^i^1  + niVoZVt  exP(-E/kT) ] 1 

where  is  the  concentration  of  defects  at  time  t,  r|°  is  the  initial  concentration 
is  the  iump  frequency  of  the  defect,  z is  the  number  of  sites  from  which 
a jump  can  take  place,  v is  the  volume  per  lattice  site,  and  E is  the  defect 
migration  energy. 

—7  — R 

The  best  fit  to  the  data  yields  E = 1.4  ± 0.1  eV,  and  vqzv  =10  - 10 

13  -1 

For  an  order-of-magnitude  estimate  of  Vqzv  we  may  take  v = 10  sec  , z = 10, 

—22  —3  —8 

and  v - 10  cm  to  get  Vozv  = 10  . The  same  annealing  kinetics  are  observed 

for  the  extra  donors  created  by  heating  the  sample  above  the  reverse  annealing 
stage  and  quenching.  The  activation  energy  of  1.4  eV  is  typical  of  vacancy 
motion.  We  therefore  feel  that  this  annealing  stage  represents  sulfur-vacancy 
motion  and  annihilation  while  the  reverse  stage  represents  sulfur  evaporation. 

We  can  summarize  the  results  of  radiation  damage  effects  in  CdS  as  follows: 

1.  The  displacement  threshold  for  S in  CdS  is  about  120  ± 5 keV. 

2.  The  displacement  threshold  for  Cd  is  about  314  i 5 keV. 

3.  Irradiation  at  10°K  produces  a shallow  donor  level  near  0.03  eV 
(I  bound-exciton  line)  and  a deep  level,  probably  an  acceptor 
(1.03ym  luminescence  band). 

4.  This  shallow  donor  level,  as  well  as  a deeper  donor  level  at  0.13  eV, 
are  associated  with  the  S vacancy.  The  1.03ym  band  is  associated  with 
a defect  on  the  Cd  sublattice,  probably  the  Cd  vacancy. 

5.  In  the  temperature  range  of  60-150°K,  Frenkel  pairs  tend  to  align  along 
a preferential  axis. 


Fig.  2-47.  Isothermal  annealing  of  additional  carrier  concentration 
for  three  annealing  temperatures. 


6.  At  about  180  K,  S Interstitials  migrate  and  remove  perhaps  90%  of 
the  sulfur  vacancies.  In  addition,  a S-interstitial  complex  is 
formed  (7200X  band). 

7.  The  annealing  stage  at  230°K  is  probably  related  to  Cd-interstitial 
migration;  certainly  a Cd  sublattice  defect  is  involved. 

8.  The  n-type  conductivity  of  CdS  at  room  temperature  is  caused  by 
shallow  S-vacancies. 

9.  These  vacancies  migrate  at  250°C  with  an  activation  energy  of 
1.3-1. 5 eV. 

10.  Above  300-350°C,  additional S— vacancies  are  added  to  the  lattice  by 

a mechanism  of  S evaporation  plus  S vacancy  migration  from  the  surface. 

11.  The  dominance  of  conductivity  type  by  a vacancy  on  one  sublattice  is 
similar  to  the  ZnTe  case.  In  fact,  a survey  of  the  literature  suggests 
that  most  irradiated  II-VI  compounds  have  a below-L. ^-temperature 
Frenkel  orientation  stage,  interstitial  migration  between  77°K  and 
300°K,  vacancy  migration  stages  from  100°C-300°C,  and  evaporation 
stages  above  300°C. 

12.  Vacancies  appear  to  be  stable  in  all  the  II-VI' s at  room  temperature. 

In  the  next  sections,  some  of  the  above  ideas  are  confirmed  in  ZnO  and 

ZnS.  Watkins  has  found  similar  results  in  ZnSe.  Although  many  details  are 
not  fully  understood  the  general  picture  of  defect  behavior  in  the  II-VI 's 
appears  to  be  fairly  well  established. 

2. 3. 1.3  Radiation  Damage  in  ZnO 

Radiation  effects  in  ZnO  are  of  interest  because  of 

33 

the  technological  importance  of  this  material.  Laser  transitions  in  the 
UV  have  been  stimulated  in  this  material.  Tts  photoconductivity  properties 
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form  the  basis  of  a copy  machine  process.  ZnO  has  been  used  in  a variety  of 
light-emission  devices,  from  digital  displays  to  TV  picthre  tubes.  ZnO  has 
also  been  proposed  as  a useful  material  for  the  fabrication  of  high-frequency 
delay-line  transducers  and  acoustic  amplifiers  because  of  its  high  electro- 
mechanical coupling  constant  and  low  dielectric  constant. 

One  of  the  oldest  uses  of  ZnO  is  as  a white  paint  pigment.  Radiation 
effects  on  ZnO  used  in  this  application  are  of  particular  interest  since  this 

pigment,  in  combination  with  others,  has  been  applied  as  a thermal  paint  coating 
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to  more  spacecraft  than  any  other  paint. 

Oxides  have  been  found  to  be  some  of  the  most  radiation-resistant  materials 
34 

known.  Chen  et  al.  have  been  able  to  observe,  by  means  of  optical  absorption, 

the  oxygen  displacement  threshold  in  MgO  at  330  keV.  This  value  corresponds 

to  a recoil  energy  T,  = 60  eV.  This  value  is  approximately  six  times  as 
d 

large  as  the  displacement  energies  observed  in  the  IX-VI  semiconductors. 

In  the  case  of  ZnO,  early  estimates  of  T^  were  much  to  high  when  compared 

35 

with  displacement  energies  in  other  materials.  Vehse,  Sibley,  Keller  and 
Chen  irradiated  ZnO  with  0.6  MeV  and  1.7  MeV  electrons.  They  observed  that  the 
optical-absorption  coefficient  at  4100X  increased  for  the  higher  energy  irradia- 
tion. They  were  able  to  correlate  this  absorption  band  with  the  yellow  color 

35 

centers  produced  by  heat  treatment  of  ZnO  in  Zn-vapor. 

They  have  also  shown  that  elastic  collisions  are  the  major  source  of 
radiation  damage  in  ZnO.  Since  they  observed  no  coloration  in  ZnO  for  600  keV 
irradiation,  they  assumed  that  the  displacement  threshold  electron  energy, 
could  be  no  lower  than  600  keV.  This  corresponds  to  a recoil  energy,  T^ , of  at 
least  30  eV  for  Zn  displacements,  or  120  eV  for  oxygen  displacements.  They  argued 
that  120  eV  was  too  high  a value  for  T^  and  therefore  assumed  that  they  were 
observing  Zn  displacements  (a  reasonable  assumption) . They  suggested  that  the 
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4100X  absorption  band  was  associated  with  Zn  interstitials  which  were  produced 
either  by  heat  treatment  in  Zn  vapor  or  by  electron  irradiation  above  600  keV. 

It  has  been  difficult  to  interpret  subsequent  experiments  in  terms  of 
the  abov6  model.  Careful  channeling  experiments  by  Appleton  and  Feldmen  on 
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ZnO,  additively  colored  in  Zn-vapor,  have  failed  to  detect  Zn  interstitials. 

More  surprisingly,  it  was  found  that  ionized  oxygen  vacancies  (F+  centers) 

37 

could  not  be  observed  by  EPR  for  600  keV  electron  irradiations.  If  the  Zn 

threshold  is  approximately  600  keV,  as  orginally  suggested,  then  the  oxygen 

threshold  should  have  been  much  lower  because  of  its  lower  mass.  Therefore, 

600  keV  irradiations  should  have  produced  many  oxygen  vacancies.  Smith  and 
37  + 

Vehse  have  also  found  that  the  F resonance  thermally  anneals  in  the  same 
temperature  range  as  the  color  center  identified  in  the  previous  work.~^ 

The  conclusions  concerning  the  displacement  energies  in  ZnO  to  be  drawn 
from  these  early  studies  are  not  satisfying.  The  work  of  Smith  and  Vehse, 
and  also  that  of  Appleton  and  Feldman,  suggest  that  the  4100&  absorption  band 
should  be  associated  with  oxygen  vacancies,  not  zinc  interstitials.  This 
implies,  however,  that  T^  for  oxygen  in  ZnO  must  be  larger  than  120  eV  since 
no  EPR  resonance  for  the  F+  center  is  observed  after  600  keV  electron 
irradiation.  This  conclusion  is  unsatisfactory  since  T^  is  of  the  order  of 
60-80  eV  in  diamond,  a material  with  a formation  energy  nearly  twice  that  of 
ZnO.  The  work  performed  under  this  contract  has  resolved  many  of  these 
difficulties.  In  addition,  valuable  information  has  been  obtained  from  which 
it  should  be  possible  to  prevent  much  of  the  yellowing  associated  with  ZnO 
paint  pigment  exposed  to  a space  radiation  environment  whose  degradation,  is 
manifested  in  higher  spacecraft  cooling  system  loads. 


in  air  for  16  hrs.  at  5-10)jA/cm  in  steps  of  50  keV  from  0.6  to  1.0  MeV.  The 
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The  investigation  of  irradiation  effects  began  with  a study  of  degradation 


of  electrical  properties  of  ZnO  during  room-temperature  irradiations.  The 


conductivity  for  fixed  beam  energies  was  found  to  decrease  linearly  with  fluencn 
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for  fluences  up  to  3 x 10  e/cm  so  long  as  the  beam  power  density  was  kept 


constant.  The  carrier  removal  rate,  dn/d0,  was  found  to  vary  quadrat ically  with 


bombardment  energy  above  a displacement  threshold  of  310  t 20  keV.  This  Is 


shown  in  Fig.  2-49  where  the  square  root  of  carrier  removal  rate  (dn/d0) 


vs.  beam  energy  is  shown 


A second  displacement  threshold  is  observed  at  higher  energies  by  the 


onset  ot  three  phenomena.  The  sampj.es  used  in  this  experiment  were  irradiated 


EPR  spectrum  was  taken  after  each  irradiation.  For  energies  above  900  keV, 


a sharp  reduction  in  the  free-carrier  microwave  absorption  is  observed. 


3+  2+ 

Accompanying  this,  the  EPR  resonances  for  Fe  and  Mn  are  observed  to  increase 


drastically  due  to  the  shift  in  Fermi  level  detected  by  the  reduction  in 


microwave  absorption.  The  F center  is  also  observed  for  the  first  time  and 


is  accompanied  by  the  sample  taking  on  a yellow  color  due  to  an  increase  in  the 


4100X  absorption  band.  We  emphasize  that  only  F centers  (ionized  oxygen 


vacancies)  can  be  seen  by  EPR.  The  neutral  F center  would  not  be  observed  if 


present.  It  is  also  not  possible  to  make  optical  absorption  transitions  to  energy 


levels  which  are  occupied. 


To  understand  the  relationships  between  these  experiments,  the  crystal  was 


thermally  annealed  in  air  at  300  C for  30  min.  This  procedure  is  known  to  remove 


most  of  the  color  centers  produced  by  neutron  or  electron  irradiation.  After 


this  anneal,  the  crystal  became  clear,  indicating  the  removal  of  the  color 


centers.  The  microwave  free  carrier  absorption,  did  not  recover,  however, 
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Fij>.  2-49.  Square  root  of  the  carrier  removal  rate  vs.  beam  energy 
for  ZnO  after  a room-temperature  irradiation. 
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spectrum  of  ZnO  irradiated  at  309  K by  600  keV  electrons 
r a previous  irradiation  at  1 MeV  and  300°C  thermal  anneal 


(PRODUCTION  RATE  OF  F*  CENTERS) 


200  300  400  500  600  700 

BOMBARDMENT  ENERGY  (keV) 

Fig.  2-51.  Square  root  of  the  production  rate  of  ionized  oxygen 
vacancies  (F+  Centers)  versus  beam  energy. 
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being  displaced.  Using  F^  (oxygen)  = 310  keV  and  E^(zinc)  - 900  keV,  we 

find  that  1,(0)  =»  T , (Zn)  = 57  eV.  Since  T,  represents  the  maximum  kinetic 
da  d 

energy  transferred  to  a lattice  ion  by  the  minimum  electron  energy  necessary 
to  produced  displacements,  we  would  expect  the  values  of  to  be  nearly 
equal  for  both  sublattices, as  observed.  This  is  to  be  expected  since  T^  is 
the  energy  required  to  break  the  ionic  bonding  and  this  energy  must  be  nearly 
the  same  for  either  ion. 

The  time  constant  for  the  decay  of  the  F+  resonance,  after  removal  of 
white  excitation  light,  has  also  been  measured.  The  results  are  shown  in 
Fig.  2-52.  The  decrease  in  intensity  is  composed  of  a rapid  decay  followed 
by  a slow  decay  which  is  exponential  with  a time  constant  of  about  1 sec.  We 
believe  that  this  double  exponential  decay  indicates  that  the  zinc  vacancy  is 

a double  acceptor.  A model  for  this  double  acceptor  has  been  presented  else- 
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where.  A further  confirmation  of  this  modei  is  found  by  studying  the  optical 
excitation  of  the  F+  resonance.  Two  excitation  peaks  are  observed  at  5300X 
and  ~4000X.  These  may  indicate  the  positions  of  the  zinc-vacancy  levels  or 
even  the  oxygen-vacancv  levels.  These  transitions  are  clearly  irradiation-defect 
related. 

It  should  be  noted  again  that  the  oxygen  vacancy  is  stable  in  ZnO  at  room 
temnerature,  and  anneals  in  the  300°C  range  as  was  the  case  for  ZnTe,  CdS, 
and  Uatkin's  work  on  ZnSe.  Chen  et  al.  have  also  observed  an  annealing  stage 
between  L.N^  ar.d  room  temperature  in  which  some  90%  of  the  F-center  vacancies 
anneal.  This  stage  is  probably  the  migration  stage  for  the  oxygen  interstitial, 
a situation  similar  to  CdS.  Thus  a pattern  of  defect  migration  energies  is 
beginning  to  emerge  for  all  the  II-VI  compounds  investigated. 


i 
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2. 3. 1.4  Radiation  Effects  in  ZnS 


These  experiments  represent  an  attempt  to  exploit  the 
suspected  similarity  between  What  was  firmly  exstablished  in  the  previous 
section  for  ZnO  and  some  previous  experiments  which  have  been  unexplained 
ir  ZnS.  FPR  had  been  previously  used  to  conclusively  identify  the  sulfur 

vacancy  (P+  centers)  in  ZnS  produced  by  heat  treatment  in  a Zn  atmosphere, 

39 

and  by  neutron  irradiation  at  room  temperature.  ' Note  again  that  the  zinc 

vacancy  is  stable  at  room  temperature  in  ZnS.  Furthermore,  an  annealing 

stage  at  about  240°X  was  observed  previously  from  luminescence  data  but  was 
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not  conclusively  identified.  It  had  also  been  reported  that  1 Mev  electron 

+ 39 

irradiation  at  room  temperature  produced  no  F centers  in  ZnS  and  it 

was  suggested  that  the  displacement  thresholds  were  very  high  (a  conclusion 

which  seems  faulty  in  terms  of  the  difficulty  of  observing  F+  centers  in  ZnO) . 

Single  crystals  of  ZnS  were  irradiated  at  77°K  with  1 MeV  electrons  to 
18  2 

a tluence  of  about  10  e/'cm  . Several  impurity  and  defect  spectra  were 
observed  and  identified  as  being  modified  or  produced  bv  irradiation.  Addition- 
ally, a center  though  to  be  the  F+  center  was  observed.  Nearly  100%  recovery 
of  this  resonance  was  observed  to  anneal  at  240°K  in  agreement  with  previous 
optical  measurements.  We  believe  that  this  experiment  explains  the  difficulty 

in  obtaining  F+  centers  in  ZnS  by  1 MeV  irradiation  at  room  temperature.  Since 

•f  o *1* 

most  F centers  are  lest  in  this  annealing  stage  at  240  K,  but  since  the  F 

can  also  be  observed  after  neutron  irradiation  at  room  temperature,  we  tentatively 

conclude  that  the  S interstitial  migrates  near  240°K  but  that  the  vacancy  does 

not  rove  until  above  room  temperature. 

It  must  be  added,  however,  that  in  the  case  of  ZnS,  these  conclusions  are 

very  tentative  since  not  a large  enough  co-  centration  of  F+  centers  could  be 


produced  to  see  the  hyperfine  structure  necessary  for  positive  defect  identi- 
fication. The  reasons  for  this  an  not  understood  at  this  time. 

2. 3. 1.5  Summary  of  Radiation  Damage  Experiments  in  II-VI  Compounds 
Although  much  work  remains  to  be  done  to  completely 
identify  all  the  radiation-produced  effects  in  Il-VI's,  a pattern  of  damage 
effects  has  begun  to  emerge.  The  low-temperature  annealing  which  has  been 
observed  below  L.^  temperatures  has  been  found  to  be  associated  with  Frenkel- 
pair  reorientation  whenever  these  stages  have  been  observed.  It  is  expected 
that  large  recovery  stages  in  which  interstitials  migrate  are  found  in  most  of 
the  II-VI 's  in  the  temperature  range  near  200°K.  Vacancies  have  always  been 
found  to  be  stable  in  the  II-VI' s at  room  temperature  and  begin  to  migrate  from 
100-300°C.  Beyond  this  range,  evaporation  of  one  of  the  cor.stitutents  of  the 
crystal  usually  is  found.  The  displacement  thresholds  in  the  II-VI  compounds 
are  generally  in  the  8 eV  range  for  T^.  ZnO  is  the  one  exception . The  very 
high  value  of  T^  for  ZnO  can  be  understood  in  terms  of  the  unique  strength 
of  the  oxygen  bond  which  places  the  thresholds  in  ZnO  near  those  found  for  other 
oxide?,  regardless  of  structure.  The  high  radiation  hardness  of  II-VI's  can 
be  inferred  from  the  above  experiments.  Although  T^  is  generally  lower  in  the 
compounds  than  the  elemental  semiconductors,  it  is  now  apparent  that  Frenkel 
pairs  do  exist  in  the  II-VI's  and  that  defect  migration  is  much  more  difficult. 

The  annealing  stages  found  around  200  K in  the  compounds  tend  to  have  large - 
percentage  recoveries.  This  is  interpreted  by  us  as  an  indication  that  most 
interstitials  find  their  parent  vacancies  in  this  stage,  through  correlated 
annihilation.  This  one  fact  accounts  for  the  high  radiation  hardness  found 
for  compound  semiconductor  devices.  As  a comparison,  consider  the  situation 
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in  Ge  and  Si.  In  these  materials,  the  interstitials  are  thought  to  migrate 

freely  at  temperatures  as  low  as  10°K.  The  vacancies  and  interstitial  separate 

immediately  upon  irradiation  at  room  temperature.  Furthermore,  the  vacancies 

in  the  elemental  semiconductors  also  migrate  below  room  temperature,  as 

opposed  to  the  compounds.  is  a result  of  this  migration,  Frenkel  pairs  are 

not  found,  vacancies  and  interstitials  have  difficulty  finding  each  other, 

and  a whole  spectrum  of  defect  impurity  complexes  form  which  become  difficult 

to  anneal.  Additional  evidence  for  this  model  is  found  in  the  experimental 

fact  that  SI  and  Ge  can  be  driven  amorphorous  by  ion  implantation  more  easily 

than  the  II-VI’s  can.  This  implies  that  a considerable  amount  of  self-healing 

occurs  in  IT-VI's  which  are  irradiated  at  room  temperature.  In  fact,  the 
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defect  concentrations  found  in  most  II-VI's  after  1 MeV  fluences  of  10  e/cm 
are  comparable  to  the  concentrations  produced  in  Si  after  only  10^  to  lO^e/cm^. 

Although  most  of  the  radiation  damage  community  is  under  the  impression 
that  the  III-V's  in  general,  and  GaAs  in  particular,  should  be  similar  to  Ge 
and  ci,  strong  evidence  exists  to  suggest  that  the  III-VI's  behave  more  like 
II-VI’s  It  will  be  shown  in  a later  section  that  GaAs  solar  cells  are 
considerably  harder  than  Si  solar  cells  of  similar  efficiency.  It  has  also 

9 

been  observed  that  GaP  starlight  detectors  show  no  degradation  after  10 
rad  °^Co  gamma  irtadiation.  A study  of  the  1TI-V  radiation  damage  literature 
suggests  that  the  annealing  stages  are  found  in  temperature  ranges  similar 
to  the  TI-VI's,  not  31  and  Ge!  It  Is  apparent  that  advantages  are  to  be  gained 
by  studying  radiation  effects  in  a consistent  system  of  semiconductors,  such 
as  the  II-V1 ' s, since  the  intercomparisons  of  defect  behaviors  are  lost  when 
attention  is  confined  to  a single  element  such  as  SI. 

2.1.2  Defect  Semiconductor  Radiation  Damage  - Ii^Tej 

The  III, -VI  binary  compounds  such  as  In-Te„,  Ga0S_,  Ga»Se,,  and 
Z j Z 5 Z _5  5 
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C.a0Te.  form  a class  of  semiconducting  compounds  with  unusual  physical  properties. 

These  crystals  grow  in  a tetrahedral  phase,  based  either  on  the  zinc  blend 

3 

or  wurtize  structures,  because  of  the  tendency  toward  sp  hydrldization  of 

the  crystal  bends.  Valency  considerations  require,  however,  that  there  be 

only  two  cations  for  eacn  group  of  three  anions.  Because  of  the  tetrahedral 

crystal  habit,  this  implies  tnat  one-third  of  all  available  cation  sublattice 

Sites  are  vacant,  thereby  producing  a stiochiomecric  cation  vacant-site 

21  -3 

concentration  of  5.5  x 10  effl 

Since  these  "stoichiometric  vacancies"  are  a structural  component  of  the 
crystal,  they  do  not  carry  a net  charge.  Nevertheless,  this  unusual  structural 
element  predominates  over  many  or  the  physical  properties  of  these  crystals. 

These  "stoichiometric  vacancies"  produce  an  extremely  low  thermal  conductivity 
which  is  associated  with  the  high  anharmonicity  of  the  asymmetric  vibrations 
of  atoms  wl'.ich  surround  the  vacant  sites.  Even  heavily  doped  InjTe  exhibits 
intrinsic  conductivity  to  very  low  temperatures.  The  stoichiometric  vacancies 
and  cations  cn  the  In  sublattice  can  be  either  ordered  or  disordered  but  the 
conductivity  of  both  is  intrinsic. 


Attempts  to  alter  the.  conductivity  of  I^Te^  by  irradiation  near  room 
temperature  with  large  fluences  of  100  MeV  electrons,  1 MeV  neutrons  and 


Co  gammas  have  been  unsuccessful. 


c 41,42 

RlU'.rPQR,  11,  . 


This  is  an  interesting  result 


suggesting  that  further  study  of  the  radiation-damage  properties  of  these  defect 
compounds  could  eventually  lead  to  a new  approach  to  the  radiation  hardening 
of  semiconducting  devices.  It  has  been  suggested  that  non-equilibrium  point- 
defect  concentrations  cannot  be  preserved  in  the  defect  tetrahedral  semiconductors 


by  irradiation  of  quenching.  ’ In  the  present  experiments,  we  have  attempted 
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to  test  this  hypothesis  under  the  conditions  of  low-temperature  irradiation 


In»Te_  grown  from  stoichiometric  amounts  of  very  high  purity  In  and  Te.  The 


starting  material  was  first  reacted  in  vacuum  in  a sealed  quartz  ampoule  at 


700  C.  The  temperature  of  the  reacted  material  was  then  raised  to  800  C 
for  24  hrs  followed  by  a slow  cool  at  a rate  of  2-5°C/hr  to  550°C.  The  ingot 


was  then  rapidly  cooled  to  room  temperature.  The  ingot  formed  was  3 x 1 x 0.5  cm 


and  consisted  of  dark  crystallites  about  1 nm  . Near  the  center  of  the  ingot 


larger  single  crystals  were  found  with  the  largest  being  Icmx3mmx3mm 


It  was  estimated  from  inspection  of  the  quartz  ampoule  that  the  crystals  were 


nearly  stoichiometric  since  the  amounts  of  starting  materials  were  in 


stiochiometric  proportions  to  better  than  t 0.02  atomic  % and  since  no  evidence 


of  an  unreacted  residue  was  found  in  the  ampoule 


An  extensive  investigation  of  the  physical  proerties  was  performed  to 


insure  that  we  had  indeed  grown  ln,Te.  rather  than  InTe,  In.Te, , or  one  of  the 


Table  2-3  lists  physical 


other  stoichiometric  combinations  of  In  and  Te 


properties  of  the  a(ordered)  and  B(disordered)  phases  of  In.Te.,  compared  to 


4 lists  the  properties  of  other  possible  compounds  from 


our  samples.  Table 


In  and  Te  which  are  likely  to  form  during  growth.  It  is  clear  from  these  two 


tables  that  our  samples  were  1;:  , 1 e and  that  they  possessed  some  degree  of 


ordering  of  the  vacant  sites.  Figure  2-53  shows  the  intrinsic  electrical 


properties  of  our  sample  before  irradiation 


Figure  2-54  shows  the  77  K KPR  spectrum  of  In0Tea  before  and  after  77  K 


elect  ron  irradiation  of  1 MeV  to  a fluence  of  10 e m . In  addition  to  the 


very  broad  resonance  which  is  probably  the  result  of  a small  percentage  of 


In  on  the  vacant  sites,  four  new  resonances  labeled  A,  B1 , B,  and  C appear 


Density 


5.79  g/cm 


5.73  g/cm 


5.772  g/cm" 


Constant 

18.50  A 

6.158  A 

6.158  + .010  A with 
superstructure  lines 
corresponding  to 
18.45  A 

Cleavage 

<111> 

<111> 

<111> 

Electrical 

Conductivity 


“10"V1cm"1 


~10  3cm  3 


1.5  X 10-7.Q-1cm  1 


at  300  K 

(intrinsic) 

(intrinsic) 

(intrinsic) 

Thermal 

Conductivity 

- | 

2.8  X 10-3 

1 

1 X 10-3 

~2  X 10"3 

at  300  K 

cal/cm-sec-deg 

cal/ cm* sec-de g 

cal/cm-sec-deg 

Band  Gap 

1.11  eV 

-1.00  eV 

1.11  eV 

Thermo- 

0 mV/deg  at  59.75 

at.  Z Te 

Electric 

-650  mV/deg 

-•-600  mV/deg 

<|t  3 mV/deg  | 

Power 

(59.3  at  % Te) 

61  at  % Te) 

Mobility 

■~50  i 10  cm^V  ^sec 

-1 

2 -1  -1 
41.2  cm  V sec 

Mobility 


Table  2-4.  Physical  properties  of  other 
possible  In^Ta^  ^ compounds. 


MATERIAL 

DENSITY 

(gtn/cm^) 

LATTICE 

CONSTANT 

(A) 

CONDUCTIVITY 

(ft-lcm“l) 

OUR  SAMPLE 

5.772 

6.158  and 
additional 
18.50  lines 

1.5  X 10-7 

ot-In^Te, 

A.  J 

5.79 

18.50 

lO"7 

B-InjTe^ 

5.73 

6.158 

vO 

1 

O 

rH 

InTe 

6.29 

8.437  a 
7. 139  b 
High  pressure 
phase  exists 
which  is  cubic 
with  a = 6.16, 
not  stable  at 
STP 

14  X 10~3 

In3Te4 

5.88 

4.26 

? 

— ■ - 


The  defect-production  rate  at  77  K,  as  estimated  from  the  A resonance,  is  cf 
the  order  of  0.1  defects/1  MeV  electron.  This  production  rate  is  typical 
of  1 MeV  defect  production  rates  in  other  semiconductors. 

No  positive  defect  identifications  can  be  made  from  these  EPR  spectra 
because  the  polycrystalline  nature  of  the  sample  removes  most  of  the  structural 
information  needed  for  such  an  analysis.  It  should  be  noted,  however,  that 
singly  charged  vacancies  on  the  anion  sublattice  might  be  expected  to  be  axiai 
along  a [111]  direction  while  the  cation  vacancies  would  not.  It  is  therefore 
possible  that  the  A resonance  is  a Te  vacancy  while  the  other  resonances  ara 
possibly  associated  with  In  sublattice  defects.  Further  investigation  of 
single  crystal  In2Te^  is  needed  to  determine  this  with  certainty. 

Isochronal-annealing  experiments  were  conducted  from  100  to  300°K.  Figure 


2-55  shows  the  recoveries  of  the  resonances  labeled  in  the  previous  figure. 

Some  annealing  of  all  the  resonances  possibly  occurs  at  lower  temperatures  since 


the  intensities  all  decrease  from  the  100  K reference.  If  the  stoichiometric 
vacancies  on  the  cation  sublattice  tend  to  preferentially  lower  the  radiation- 
induced  defect  stability  on  this  sublattice,  then  annealing  of  these  defects 
at  a lower  temperature  is  in  agreement  with  the  sublattice  assignments  suggested 
above.  The  annealing  of  the  A resonance  is  extremely  broad  suggesting  that 
a fundamentally  different  annealing  mechanism  may  be  occurring  ir.  this  material. 

It  should  be  noted  that  all  the  resonances  have  recovered  by  room  temperature, 
and  that  at  200°K,  only  about  20%  of  the  A resonance  remains.  It  should  also 
be  noted  that  no  new  resonances  are  formed  as  a result  of  annealing.  This 
suggests  that  long-range  defect  migration  and  the  formation  of  new  defect 
complexes  is  not  occurring.  These  observations  also  explain  why  radiation 
damage  was  not  detected  in  the  previous  expel iments  in  Ref.  42  and  Ref.  43. 
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When  lightly  doped  (n  = 101  / cm  ) n-type  Ge  is  electron  irradiated  at 
4.2°K,  acceptor  levels  are  produced.  Those  acceptors  account  for  approximately 
90%  of  all  the  radiation  damage  in  this  material.  These  acceptors  can  be 
removed  following  irradiation  by  the  following  mechanisms: 

1.  Thermal  annealing  at  65°K. 

2.  Radiation  annealing  at  7-10°K  where  the  electron-beam  energy  stimulating 
the  annealing  is  below  the  displacement  energy. 

3.  Intrinsic  light  (hw>E^)  optical  annealing  at  about  27-30°K. 

4.  Extrinsic  light  (hw<E  ) optical  annealing  at  4.5°K  or  lower. 

g 

Although  it  has  been  known  for  many  years  that  vacancy-interstitial 

recombination  is  occurring  during  these  various  annealing  mechanisms,  the 

particular  nature  of  the  defect  migration  and  annealing  mechanisms  has  not  been 

well  understood.  However,  the  results  outlined  above  made  it  clear  that  the 

charge  states  of  the  defects  were  important  in  determining  their  stability. 

Although  it  was  a common  belief  that  these  defects  were  closely  spaced 

Frenkel  pairs,  which  could  be  annealed  by  the  above  mechanisms,  in  1970  it 

was  discovered  that  for  all  of  the  above  annealing  stages,  the  annealing  rates 
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were  dependent  on  the  type  of  n-dopant  in  the  crystal.  This  result  suggested 

that  long-range  interstitial  migration  might  be  occurring  at  low  temperatures 

in  Ge  as  had  been  suggested  for  Si  at  a much  earlier  date. 

A model  for  the  annealing  kinetics  for  all  the  above  annealing  mechanisms, 

and  based  upon  defect  charge  states  and  a freely  migrating  interstitial,  has 

44 

been  published  during  this  contract.  Although  several  close-pair  models 
exist,  the  present  model  is  tbe  only  one  from  which  all  the  annealing-rate 
constants  can  be  calculated  from  first  principles.  Only  a brief  sketch  will 


be  presented  here,  since  many  of  the  calculations  are  quite  lengthy.  Further 
details  may  be  obtained  from  Ref.  44. 
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Figure  2-56  lists  the  defects  and  charge  states  assumed  to  play  a role 


in  the  various  annealing  mechanisms.  In  the  following  treatment,  V represents 


the  vacancy  concentration,  I the  free  interstitial  concentration,  and  C 


the  concentration  of  interstitials  trapped  at  impurity  sites  of  concentration  i 


Figure  2-57  shows  the  assumed  kinetics  equations  for  these  various 


concentrations  as  a function  of  time.  The  rate  constant  K, , and  K~  determine 


been  to  find  approximate  solutions  to  the  system  of  differential  equations 


(exact  solutions  in  closed  form  do  not  exist)  and  to  calculate  the  rate 


constants  K. , K„  and  K_  from  first  principles  for  all  the  annealing  situations 


mentioned  above.  That  is  is  possible  to  accomplish  this  for  the  multitude 


of  experimental  data  available  in  the  literature  tends  to  lend  strong  support 


to  the  essential  features  of  the  model 


Two  approximate  solutions  are  possible  for  the  fraction  of  vacancies 


f - V/V  remaining  as  a function  of  time.  The  approximations  are  given  as 


I and  II  in  Fig.  2-57.  A comparison  of  the  rate  constants  calculated  from 


the  carrier  trapping  model  is  comnared  to  the  experimental  results  in 


Fig.  2-58.  It  is  clear  that  the  model  is  capable  of  predicting  the  rate  constants 


which  can  be  measured  to  within  a factor  of  2 or  so.  This  is  remarkable  since 


the  rate  constants  can  conceivably  vary  over  many  orders  of  magnitude 


Figure  2-59  shows  a fit  of  the  approximate  solution  I in  Fig.  2-57  where 


the  rate  constant  is  calculated  from  first  principles  for  a typical  high 


temperature  thermal  annealing  case 


The  optical  annealing  case  at  4.5  K for  light  of  lower  energy  than  the 


band  gap  is  shown  in  Fig.  2-60.  Again,  the  fit  is  calculated.  Figure  2-61 


shows  the  fit  of  calculated  and  experimental,  annealing  data  for  30  K intrinsic 


Fi-8*  2-56.  Processes  and  defect  charge  states  for  low-temperature 
annealing  in  n-type  Ge. 


KINETICS  EQUATIONS  FOR  MODEL 


dC/dt  = - K.C  + K0  (i  -C)  I 
1 2 o 

dV/dt  - - K3IV 


1 V = I + C 

where V =V°+  v"  + v“",  I = 1°  + l + 


= PROBABILITY  FOR  COMPLEX  DISSOCIATION/UNIT  TIME 
K2  = P|(o)  crV) 

K3'{Pl(+D  PV  Hl  %(-!) +FIQI  PVMoV(-2)}  Vl 


TWO  APPROXIMATE  SOLUTIONS 


STEADY  STATE  SOLUTION  (LARGE  RETRAPPING  PROBABILITY) 

KJK<<\  , V «i  (at  t=0,  0.15  i ) 

1 2 o o o 

then  dC/dt  0 

— dV/dt w-(M  ) V3  (THIRD  ORDER) 

i_Qj \ ,y 

f = VA/  = [l  + (2K-V  /i  )t] 

o * 3 o o J 

TRANSIENT  SOLUTION  (NEGLECT  RETRAPPING) 

near  t=o,  K9(i  -C)  I «K,(i  -C)  (V-C)~0 
2 o 2 o 


-(u-  u) 


f =V/V 


-J. o , E (x)  = T -1  -p 

euo  [E.(u)-E  (u  )]  1 JO** P 

0 L 1 l 0 1 X 


u = LI  e"Klt  U « K_V  IK  . 
o o 3 o 1 

Fig.  2-57.  Kinetics  equations  for  annealing  models. 
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COMPARISON  WITH  EXPERIMENT 


RADIATION  ANNEALING  AT  7-8°  K 


K . = 2K,V2/i  = 8.5xl0'4  sec  _1 
CdlC  3 0 0 

K = 11.2xl0_4sec  _1 
exp 


MONOCHROMATIC  LIGHT  AT  4.2°K 

Kca|c  — 10  6sec  1 or  1000  times  slower  than  radiation  annealing 


MONOCHROMATIC  LIGHT  AT  30°K 

P-Ap  = aIQT  = 5xl09  cm  a = 5 cm  1 IQ  = 104  photons/cm2  - sec 

ip  ~10  '’sec  at  30°  K 

1/  -1/  - r-..  . / L \ _ O r..in  ^ ” 1 


l/t  = pvhoh(+)  = 2.5x10  sec 


1.67x10'^  sec’1 


14  -3 

n=2xl0  cm  due  to  thermal  ionization 
P (-1)^1  Pj  (+1)  • 2.5x10  ® 


“ ICV  /K.  ‘1.7 
°calc.  J 0 1 


uo  -1.0 

exp 

2-58.  Comparison  of  calculated  and  experimental  rate,  constants. 
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FRACTIONS  OF  DEFECTS  UNANNEALED 


.’/(FRACTION  REMAINING) 


Fig.  2-60.  Calculated  fit  to  extrinsic-light  optical  annealing  at  4.5°K. 
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light  stimulation.  The  solid  annealing  curve  is  derived  from  approximate 
solution  II  of  Fig.  2-57  for  this  case. 

Although  the  kinetics  problem  outlined  above  cannot  be  solved  in  closed 

form,  a special  case  of  the  above  model  has  been  found  to  have  a definite 
45 

solution.  This  solution  will  be  presented  here  since  it  has  a wider 
application  to  a number  of  possible  situations  in  a number  of  materials. 

In  a number  of  irradiated  semiconductor  systems,  interstitials  migrate 
at  a lower  temperature  than  vacancies.  Although  a large  number  of  interstitials 
annihilate  vacancies  in  the  temperature  range  in  which  the  interstitial  begins 
to  migrate  freely,  inevitably,  some  interstitials  will  escape  this  annihilation 
to  become  trapped  at  various  sinks  which  remain  stable  up  to  higher  temperatures. 
If  these  interstitials  are  bound  in  defect  conm  1 raxes  by  Coulomb  attraction 
then  an  externally  induced  change  in  the  complexes'  charge  states,  by  photon 
or  fast  particle  irradiation,  can  release  these  interstitials  with  very  little 
probabilitv  for  retrapping.  The  kinetics  can  be  described  by  a first-order 
decrease  in  interstitial  complexes  accompanied  by  a second  order  decrease  in 
vacancy  interstitial  pairs.  A solution  in  closed  form  is  possible  even  when 
the  concentration  of  interstitials  is  less  than  the  concentration  of  vacancies. 

Let  V be  the  concentration  of  vacancies,  I the  concentration  of  free 
interstitials,  C the  concentration  of  interstitials  trapped  in  complexes,  and 
N the  concentration  of  interstitial  trapping  sites  unoccupied  by  interstitials. 
We  assume  that  V(t)  - I(t)  + C(t),  or  that  V = I + C + V where  V is  the  final 


concentration  of  vacancies  of  the  end  of  the  anneal.  V will  be  the  initial 

o 

concentration.  In  general 


dC/dt  - -KjC  + KIN 
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(2-2) 


t i 

! 


dl/dt  = KjC  - KINt  - K.2IV 

dV/dt  = -K2IV  (2—3  J 

where  att  = 0,V*V,C  = C,  and  1=0. 

o o 

If  we  assume  the  retrapping  term  +KINt  in  Eq.  2-1  can  be  neglected  due 
to  a radiation-induced  change  in  the  complex  charge  state,  which  drives 
interstitials  away  from  the  complexes,  then  K = 0. 

We  can  now  solve  Eq.  2-2  at  once  with  the  solution  given  by 

C = (2-4) 

By  substitution  of  V » I+C+V  into  Eq . 2-3  we  have 

dV/dt  = -K2IV 

= -K 2V(V-Vm-C).  (2-5) 

This  is  the  differential  equation  which  we  must  solve  in  order  to  obtain 
a complete  solution.  We  make  the  following  substitutions.  Let  y = i , 
x = K2t,  a = Kj/K2. 

The  Eq.  2-5  becomes 

dy/dx  + g (x)y  = 1 (2-6) 


where 


g(x)  = Vot  + C(x)  = Vco+  CQe 


-ax 


Equation  2-6  can  be  separated  if  we  multiply  by  the  integration  factor 

G(x)  /«Wdx 


where  g(x)  = — ^x^  . 

Thus  d[eG(x)y3  = eG(x)dx 


I 
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G(x)  G(o) 
or  e y - e y 


G(x), 
c dx. 


Since  G(0)  = 0 and 


G(x)  = j o g(x)  dx 


= j (V  + C e-aX)  dx, 
/ o <*  o 


we  can  solve  Eq.  2-7  for  y ^ = V;  i.e., 


V = y 


■1  e 


I + (*  G(x) 

V \ e dx 
o ' o 


We  summarize  the  following  solutions  for  the  various  concentrations : 

-L  «G(>) 

VQ  1+R(x) 


where 


G(x)  - V^x  + — (1  - e_0lX) , R(x)  = vl  aG(,0dx,  and  = e"01*  (2-9) 


I = (V  - V„)  - C . 


Therefore,  all  the  concentrations  are  known  as  a function  of  time  once 


R(x)  is  evaluated. 

Consider  now  the  explicit  form  of  R(x): 


R(x)  - j eG^  dx 


( efvmax]  e-c0[e  aX--l]{ 


where  v = V /a,  c * C /a, and  e = c exp (-ax).  Then 
00  °°  o o o 
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x = “ — ln(c/c  ) 
Ol  o 


dx  • - — dc 
ac 


R(x)  = R (c)  - v [c  ]Vne  expfc  I ) z_(v°°  + 1)e"Zdz 

o o o )q 


v0[co]V“  expf^Jfrf-v^,  c(x)]  - rt-v,  CQ]>.  (2-11) 


r(a,x)  = Jx  xa  e X dx 


Is  the  Incomplete  Gamma  function  which  is  tabulated  in  math  tables. 

The  solution  takes  on  a particularly  simple  form  if  V - 0.  Then  we 
consider  integrals  of  the  form 


/* 00 

),x)  * ) x x 1 e x dx  = F^  ( 


which  is  the  exponential  integral,  also  found  in  math  tables.  Closed  solution 
approximations  exist  for  E^(x)  and  are  given  by 

5 

-lnx  + ^ aixi;  0 1 x 5.  1 


/ . x +b1x+b. 

V x lnx  -y- i ; 1 < x 1 » 

X d-C^X+C^ 

where  the  accuracy  is  better  than  about  1 part  in  10^. 

This  model  i9  reapplied  to  the  optical  annealing  at  30CK  to  take  into 
account  the  fact  that  only  an  88%  recovery  is  observed  fer  this  stage.  This 
is  represented  by  V^/V^  ■ 0.12.  A plot  of  this  model  is  shown  in  Fig.  2-62 
which  should  be  compared  to  the  results  of  Fig.  2-61  which  assumed  a 100%  recovery. 

This  model  has  also  been  applied  to  a radiation-annealing  case  in  Si 

47 

which  has  never  before  been  fit  by  typical  annealing  kinetics.  The  data 


l 
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in  Fig.  2-63  are  for  float-zone  Si  where  the  trapping  site  is  presumably 

carbon.  Here,  the  parameters  and  were  adjusted  to  give  the  best  fit 

A7 

and  the  other  listed  parameters  were  from  experimental  data.  Figure  2-64 

shows  the  fit  to  radiation  annealing  in  a pulled  sample  of  high  oxygen  content. 
The  percentage  recoveries  are  quite  different  for  these  two  samples,  presumably 


due  to  different  interstitial  trapping-site  concentrations  and  types.  The 
curve  in  Fig.  2-64  has  been  fit  with  no  adjustable  parameters  and  although 

not  a perfect  fit,  it  is  approaching  the  accuracy  of  the  experiment. 
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It  has  also  recently  been  demonstrated  by  Kimerling  and  Lang  that  GaAs 
will  optically  anneal.  We  believe  this  to  be  an  important  reason  for  the 
forward-bias  injection  annealing  demonstrated  in  GaAs  lasers  and  LED's. 

2.3.4  Irradiation  Effects  in  Solar  Cells 

The  objective  of  this  work  has  been  to  determine  the  behavior 
of  various  solar  cell  designs  subject  to  fast  (1  MeV)  electron  bombardment  in 
order  to  determine  which  physical  parameter  changes  induced  by  irradiation 
will  be  most  detrimental  to  the  useful  lifetime  of  solar  cell  devices  fabricated 


from  these  materials  and  subsequently  subjected  to  a hostile  radiation  environment. 

Silicon  solar  cells  are  presently  being  used  as  the  primary  power  source 
for  Air  Force  satellites.  There  currently  exists  within  the  Air  Force  a 
need  to  improve  both  the  efficiency  of  power  conversion  of  these  cells  and  the 
device  lifetime  in  orbit  (which  degrades  primarily  through  Van  Allen  belt 
radiation  damage).  This  Air  Force  need  for  improvement  is  based  upon  anticipated 
increases  in  satellite  power  requirements  and  orbit  lifetimes  of  future  missions, 
since  present  mission  limits  are  already  dictated  primarily  by  the  power 
output  versus  orbit  time  of  the  solar  cell  arrays.  Because  a redesign  of  the 
satellite  solar  cell  arrays  for  higher  output,  achieved  by  increases  in  the  number 


115 


>l> 


RADIATION  ANNEALING  OF 
I MeV  (138°  K ) DAMAGE 
IN  DEGENERATE  n-TYPE  Si 


NF  9 ( 35°K ) 


K,  = 0.08 
K2V0=  0.12 
V0=  1.0 
Co  = 0.28 


2.0 


<{>(  e/cm2  ) - 350  KeV 


Fig.  2-63.  Fit  to  radiation  annealing  data  in  n-type  SJ 
for  a float-zone  sample. 


RADIATION  ANNEALING  OF 
I MeV  (138°  K ) DAMAGE 
IN  DEGENERATE  n-TYPE  Si 


NC  8 ( 35  ° K ) 


(X>(e/cm2)  - 350  KeV 


Fig.  2-64.  Fit  to  radiation  annealing  data  for  a pulled 
Si  sample  (no  adjustable  parameters). 
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of  cells,  would  require  a major  redesign  of  the  total  satellite  package,  at 
a great  expense,  economics  dictates  that  higher  conversion  efficiencies  and 
radiation  hardness  (both  of  which  determine  mission  lifetime)  be  obtained 
whenever  possible  from  solar  cell  arrays  compatible  with  present  satellite 
designs.  Unfortunately,  incremental  improvement  of  Si  solar  cell  array 
efficiencies  and  radiation  hardness  are  very  difficult  to  achieve  since  many 
of  the  theoretical  cell  design  limits  have  been  met. 

A recent  experimental  cell  sponsored  by  AF  Aeropropulsion  Laboratory, 

the  multi-vertical-junction  Si  cell,  is  designed  to  optimize  minority-carrier 

collection  efficiency  as  well  as  to  reduce  photon  surface  reflection.  Since 

the  generation  of  minority  carriers  in  the  vertical-junction  cell  occurs  in 

closer  proximity  to  a junction  element  than  in  conventional  solar  cell  designs, 

an  added  benefit  of  the  geometrical  advantage  of  this  design  should  be  less 

sensitivity  to  the  radiation  induced  drop  in  minority  carrier  lifetime  and 

diffusion  length.  We  have  recently  irradiated  some  of  these  experimental 

16  2 

cells  with  1 MeV  electrons  to  a total  dose  of  1 x 10  - 1 MeV  e/cm  and  find 

that  the  radiation  hardness  of  these  cells  is  considerably  improved  over  more 
conventional  Si  solar  cell  designs.  (Data  on  the  performance  of  these  cells 
and  radiation  degradation  may  be  obtained  by  contacting  Dr.  Patrick  Rahilly 
of  AFPL) . Since  the  efficiencies  of  these  prototype  cells  are  not  as  high 
as  eventually  anticipated,  it  is  not  possible  to  precisely  assess  the  improve- 
ment in  radiation  hardness  of  this  design  at  this  time.  Further  experiments 
on  these  cells,  particularly,  a measurement  of  the  radiation-induced  degradation 
of  the  minority  carrier  lifetime  and  a measure  of  the  unirradiated  effective- 
diffusion  length,  are  planned  for  the  near  future  as  soon  as  additional  cells 
become  available  from  the  vendor. 
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In  addition  to  the  above,  approximately  40  Si  solar  cells  of  various 
designs  have  now  been  irradiated  to  fluences  of  lO"^,  10^,  lO^  and  10^e/cm2 


for  AFPL.  Also,  over  30  Li-doped  Si  Samples  have  been  irradiated  for 

Dr.  Brown  of  Ohio  State,  at  the  request  of  AFPL,  to  10^,  10^,  2 x 10^,  and 
16  2 

3 x 10  e/cm  . These  samples  will  be  studied  by  NMR  and  EPR  at  Ohio  State  to 
determine  specific  properties  of  the  radiation-induced  defects  and  their 
annealing.  The  object  of  this  program  is  to  investigate  the  self  annealing 
induced  by  Li  doping  of  Si  in  an  attempt  to  enhance  the  radiation  hardness 
of  Si  cells  by  this  technique.  We  have  therefore  completed  approximately 
310  irradiations  for  this  study. 

The  recent  production  of  many  new  semiconductor  materials  of  reasonable 

size  and  quality  has  prompted  a search  for  alternative  solar  cell  materials 

for  both  space  and  terrestrial  applications.  Based  on  an  ideal  diode  model, 

CuInS2  appears  to  be  a material  deserving  further  experimental  study  for 

49 

solar  energy  conversion.  The  maximum  theoretical  efficiency  for  energy 

conversion  has  been  calculated  for  CuInS2 , using  the  model  of  Wysocki  and 

Rappaport,^  to  be  27.4  to  32.1%,  corresponding  to  minority  carfier  lifetimes 
—8  —6 

of  10  to  10  sec.  The  efficiencies  for  Si,  GaAs,  and  CdS,  based  on  this 
model,  are  24%,  28%, and  19%  respectively.  Since  both  shallow  donor  and 
acceptor  states  appear  to  be  related  to  lattice  defects  in  CuInS2,  solar 
cells  from  this  material  may  prove  to  be  more  radiation-damage  resistant 
in  a space  environment  than  Si  cells. 

Recent  interest  has  developed  in  GaAs  as  a highly  efficient  solar  cell 
material.  Figure  2-65  shows  the  behavior  of  a typical  Varian  GaAlAs-GaAs 
solar  cell  as  a function  of  concentration  ratio,  C'.  In  this  figure, "C'  = 312" 
means  a concentration  ratio  of  312  suns.  The  efficiences  of  the  GaAlAs  window 
cells  have  reached  15%  unconcentrated  and  17%  in  a concentrator  configuration. 

2 

It  bhould  be  noted  that  the  current  output  for  C'  » 312  is  approximately  lOA/cm 
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CURRENT  (AMPERES) 


From  Che  radiation-annealing  model  of  the  previous  section,  it  might  be  expected 
that  such  large  minority-carrier  currents  might  provide  an  annealing  mechanism 

for  radiation  damage  exposure  from  the  Van  Allen  Belts. 

52 

C.  E.  Barnes  has  observed  that  GaAs  laser  diodes  will  anneal  when 

forward  biased.  This  is  shown  in  Fig.  2-6G.  Note  that  the  forward-bias  current 

is  of  the  same  order  of  magnitude  as  the  output  current  for  concentration 

GaAs  cells  shown  in  rig.  2-65.  The  annealing  of  Barnes  was  after  a ^Co 
7 14  2 

irradiation  of  1 x 10  rad  or  3bout  5 x 10  1 MeV  e/cm  . Appreciable,  annealing 

is  seen  in  1 hour.  Figure  2-67  shows  additional  data  of  Barnes  for  the  anneal- 
ing as  a function  of  diode  temperature.^  Barnes  reports  that  the  230°C 
annealing  stage  in  GaAs  is  shifted  to  room  temperature  by  forward  biasing  the 
diodes.  Hum  and  Barry  have  also  assumed  this  effect  on  LED's  (GaAs)  after 
^Co,  2 MeV  neutrons,  and  100  MeV  protons.  ^ 

Since  a solar  cell  is  automatically  forward  biased  when  illuminated  (the 
injection  current  of  interest  is  approximately  the  solar  cell  short-circuit 
current),  the  possibility  exists  that  significant  annealing  of  the  space 
radiation  damage  is  to  be  expected  and  that  the  rate  of  annealing  should 
increase  as  the  output  current  of  the  cell  is  increased.  It  should  be  noted 
that  the  damage  rate  to  solar  cells  in  a space  environment  is  many  orders  of 
magnitude  lower  than  in  Barnes  experiment.  Experiments  have  been  planned  to 
test  this  effect  in  solar  cells  once  some  concentrator  cells  can  be  obtained. 
Preliminary  results  on  several  1-sun  cells  have  not  been  conclusive  to  date; 
however,  there  seems  to  be  a slight  annealing  effect  for  lightly  damaged  cells. 

Although  GaAs  solar  cells  and  their  irradiation  effects  have  been  studied 
for  over  a decade,  only  recently  have  large-area  cell  growth  and  the  very  high 
surface-recombination  problem  been  solved.  In  particular,  the  surface  recombination 
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problem  has  been  overcome  by  the  growth  of  a layer  of  greater  bandgap  GaAlAs 

over  the  cell  which  tends  to  reflect  carriers  toward  the  junction  and  away 

from  the  surface.  These  improvements  have  lead  to  the  fabrication  of  high 

efficiency  cells  which  are  serious  contenders  for  space  applications.  We 

report  here  the  first  experiments  on  the  effects  of  radiation  damage  on 

solar  cell  performance  for  GaAs  cells  with  GaAlAs  windows. 

Performance  evaluations  were  made  of  AlGaAs/GaAs  heteroface  P on  N solar 

cells,  exposed  to  0.5  MeV  and  1 MeV  electron  irradiations  with  fluences  up 
16  2 

to  6 x 10  e/cm  . Current-voltage  data  were  obtained  under  air  mass  zero 
(AMO)  sunlight  conditions.  Also  measured  were  changes  in  electron-beam- 
generated  cell  current  during  irradiation  which  provided  a measure  of  the 
relative  changes  in  the  minority-carrier  diffusion  lengths  as  a function  of 
fluence. 

The  cells  tested  were  supplied  by  Varian  Associates,  Palo  Alto,  Calif. 

These  cells  were  formed  by  liquid-phase  epitaxial  deposition  6f  n-type  GaAs 

on  n+  GaAs  substrates,  followed  by  Al^a,  ^As  deposition  and  p-dopant  diffusion 

from  the  AlGaAs  layer.  This  reduces  the  effects  of  surface  recombination 

and  enhances  minority-carrier  collection  frotr  the  p-region.  These  cells  were 

not  optimized  for  AMO  sunlight  but  did  display  AMO  conversion  efficiencies  of 

2 

12  to  13%.  Cell  sizes  were  approximately  1.5  cm  in  area  and  0.038  cm  in 
thickness.  The  contacts  used  on  these  cells  were  Ge-doped  Au.  All  measurements 

Q 

were  taken  at  25  C. 

The  electron  irradiations  were  performed  at  the  Van  de  Graaff  facility 
at  AF  Avionics  Laboratory,  while  the  AMO  and  spectral  response  data  were  obtained 
using  the  AF  Aero  Propulsion  Laboratory  Spectralab  Mark  I and  Mark  II  solar 
simulators. 


Results  for  the  1 MeV  electron-damage  experiments  are  shown  in  Fig.  2-68. 
Relative  data  ara  presented  with  initial  data  given  in  the  insert.  The  initial 
diffusion  lengths  for  the  minority  carriers  in  the  p and  n materia’  were  not 
measured.  However,  the  relative  changes  in  the  sum  of  the  diffusion  lengths 
(L^  + L + Lpo)  do  follow  a simple  relationship  given  by 


/I  + 84> 


(2-21) 


where  8 = Lq  K,  K is  defined  as  the  electron  damage  coefficient,  $ is  the 

electron  fluence,  and  L is  the  value  of  L = L -PL  for  <b  = 0 . The  solid  line 

o a p 

in  Fig.  2-68  for  (L/Lo)  was  based  on  calculations  using  Eq.  2-12  with 
-15  2 

B “ 3.43  x 10  cm  /e;  the  ex4  rimentally  determined  points  were  used  to 
obtain  8 from  a least  squares  fit. 

The  relative  changes  in  the  short-circuit  current,  I , ouen-circuit 

sc 

voltage,  V , and  maximum  power,  P , follow  the  characterisitc  form 
oc  max 


Relative  Value  = 1 - aln(l  + £$) 


(2-13) 


15  2 

until  > 10  e/cm  , where  deviations  become  quite  drastic.  The  data  suggest 
15  2 

that  beyond  <f>  = 10  e/cm  , the  V , I , and  P become  linearly  dependent 

oc  sc  max 

on  the  minority-carrier  diffusion  length. 

Figure  2-69  illustrates  the  voltage  and  current  at  maximum  power 


(V  , I ) as  a function  of  1 MeV  electron  fluence.  The  V drops  by  nearly 
mp’  mp  mp  * 7 3 

13  2 

4%  after  10  e/cm  and  then  rises.  Simultaneously  I degrades  in  a manner 

mp 

described  by  Eq.  2-13.  The  observed  behavior  of  V can  be  explained  by  series- 

resistance  effects  wherein  the  rapid  drop  in  1 reduces  the  voltage  loss 

across  the  series  resistance,  and  thus  produces  the  observed  increase  in 

V at  higher  fluences. 
mp 
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2-69.  V and  I versus  1 MeV  fluence 


w ' > w 


Figure  2-70  shows  I (4>)/I  (0)  and  L/L  versus  electron  fluence  for 

sc  sc  o 

0.5  MeV  electron  irradiations.  It  was  noted  that  the  damage  by  the  lower 
energy  electrons  becomes  significant  only  at  higher  fluence,  compared  to 
1 MeV  irradiation.  This  is  expected  since  the  displacement  cross  section 
is  lower  for  0.5  MeV  than  for  1.0  MeV  electrons  in  GaAs. 

Overall,  the  results  suggest  that  AlGaAs/GaAs  solar  cells  are  more 
radiation  resistant  to  electron  damage  than  conventional  GaAs  cells.  This 
effect  may  be  due  to  the  presence  of  A1  and  is  in  agreement  with  radiation 
effects  observed  in  both  Ga^l^^As  and  GaAs  ligl t-emitting  diodes. 

Compared  to  conventional  Si  cells,  the  GaAs  cells  are  harder  to  1 MeV 
electron  irradiation,  as  shown  in  Fig.  2-71.  When  compared  to  the  Si 
vertical-lunction  cells,  the  situation  is  not  yet  clear.  More  cells  of  both 
types  must  be  tested.  It  appears  that  the  GaAs  cell  is  definitely  ahead  of 
the  Si  vertical-junction  cell  in  efficiency,  however. 

2.3.5  Proton  Induced  Characterisic  X-Ray  Analysis 

A minor  effort  has  been  spent  on  this  contract  in  aiding 
Dr.  Ron  3erliner  of  Systems  Research  Laboratories  in  his  project  to  develop 
impurity  analysis  and  implanted  profile  techniques  based  upon  characteristic 
x-ray  emission  generated  by  proton  irradiation.  The  effort  from  this  contract 
has  been  in  converting  the  400  keV  Van  de  Graaff  from  electrons  to  protons. 


and  in  doing  some  of  the  mathematical  analysis  of  the  implanted-ion  profiles. 
We  will  not  report  extensively  on  these  subjects,  however,  since  Dr.  Berliner 
has  done  most  of  the  work.  We  reproduce  below  a copy  of  an  abstract  submitted 
to  the  March  1976  meeting  of  the  American  Physical  Society  in  Atlanta. 

Protons  with  energies  in  the  range  50-500  keV  have  been  used  to  excite 
characteristic  X-  tty  omission  from  bulk  specimens  of  Si,  diffusion  doped  with 


J 
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As.  The  X-ray  emission  from  the  As-L  and  Si-K  lines,  as  analyzed  by  a Si (Li) 
X-ray  spectrometer,  was  measured  a?  a function  of  incident  proton  energy. 


Decomposition  of  the  strongly  overlapping  spectra  was  accomplished  by  means 

54 

of  the  computer  code  "SAMPO".  The  X-ray  emission  from  bulk  specimens  under 
proton  bombardment  has  been  calculated.  These  calculations  include  the 
effects  of  proton-energy  loss  and  X-ray  absorption  and  can  be  performed  for 
any  impurity  profile.  The  specimen-impurity  profile  was  obtained  by  fitting, 
in  a least-squares  sense,  the  measured  X-ray  intensities  to  the  results  of 
the  calculations.  An  X-ray  filter  was  used  in  certain  cases  to  enhance  the 
detection  of  the  As-L  line  relative  to  the  Si-K  X-ray  line.  Estimates  of 
the  accuracy  of  these  profile  determinations  will  be  presented. 


r 
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3 . 0 NUCLEAR  MAGNETIC  RESONANCE  IN  II-VI  COMPOUNDS  AND  GaAs 
3.1  Introduction 

Nuclear  magnetic  resonance  (NMR)  techniques  have  been  used  quite 
successfully  in  elucidating,  some  of  the  conduction-electron,  paramagnetic-ion, 
anil  chemical-bonding  properties  of  the  II-VI  compounds,  CdS,  CdSe,  CdTe,  and 
CdO.  At  lease  half  of  the  present  published  literature  on  this  subject  has 
come  from  this  contract  effort.  We  have  been  able,  among  other  things,  to 
study  the  semiconductor-to-metal  transition  as  a function  of  donor  doping  in 


CdS. 


, 55 


This  type  of  study  also  leads  to  values  of  g* , the  effective  g-factor 


of  conduction  electrons.  For  example,  we  were  the  first  to  measure  g*  in 
CdTe,"’*’  following  unsuccessful  attempts  by  others  using  EPR  techniques. 

Another  area  of  applicability  of  NMR  techniques  involves  the  study  of  para- 
magnetic impurites  in  the  ppm  range.  For  example,  10  ppm  of  Co^+  in  CdS  will 
dominate  the  Cd^^^  nuclear  spin-lattice  relaxation  from  1-300°K.^  We  have 
also  investigated  the  effects  of  Cr,  Zn,  and  Te  impurites  in  GaAs.  Below  we 
describe  the  results  of  several  NMR  studies.  The  details  of  most  of  these 

may  bs  found  in  the  published  literaruje. 

2+ 

3 . 2 Study  of  Co  in  CdS 

The  following  abstr  ct  is  reprinted  from  Ref.  57. 

Nuclear  spin-lattice  relaxation  times,  T^,  have  been  measured  over 

a temperature  range  T - 1.1  - 300  K and  frequency  range  V = 2 - 15  MHz  in 

a single  crystal  of  CdS  doped  with  13 -ppm  cobalt.  Minima  in  T^  V3  T are 

observed,  and  absolute  values  of  the  effective  electron  relaxation  time,  t . may 

n —ft 

be  calculated  at  the  temperatures  of  the  minima;  e.g.  at  14  K,  Tg  = 8.0  x 10- 

sec. . At  low  temperatures,  T - 5°K,  is  dominated  by  a resonant  Orbach  process 

3 1 

involving  the  two  ground-state  Kramers  doublets  (S  = + — and  S « + y)  which, 

according  to  our  measurements,  are  split  by  A = (4  1'  1)°K  at  zero  magnetic 

field.  At  higher  temperatures,  T > 20°K,  T is  dominated  by  a nonresonant 
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process,  and  these  data  are  well  fitted  from  40  to  300°K  by  « T^J^ (210/T) , 
where  is  a transport  Integral  and  210°K  is  the  CdS  Debye  temperature. 

3.3  A Measurement  of  the  Conduction-Electron  g-Factor 
in  CdTe 

The  following  abstract  is  reprinted  from  Ref.  56. 

The  effective  g-factorf  g*  of  conduction  electrons  in  degenerate 
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CdTe  has  been  determined  by  using  measurements  of  the  Cd  nuclear  spin- 

lattice  relaxation  time  T^  and  the  Knight  shift  K.  It  is  shown  that  the 

2 2 

magnitude  of  g*  is  given  by  the  Korringa  product  T^TK  = C(g*)  , where  C 

is  a known  constant  and  T is  the  absolute  temperature,  and  that  the  sign 

of  g*  is  given  by  the  sign  of  K for  a spherically  symmetric  conduction  band. 

The  measured  value,  g*  = - 1.1  t 0.1,  is  within  the  range  allowed  by  effective- 

mass  theory.  Also,  the  electronic  probability  density  at  the  nucleus, 

2 25 

normalized  to  unity  in  an  atomic  volume,  is  calculated  to  be  1^(0))  - 6.5x10 

-3  2 

cm  , about  70%  of  that  found  for  the  free  Cd  ion  in  a 5s  state. 

3-4  Studies  of  Paramagnetic  Dopants  in  GaAs 
69 

The  Oa  nuclear  spin-lattice  relaxation  times  (T^)  have  been 
measured  in  GaAs:Te  over  a temperature  range  T = 6-300°K  and  a frequency  range 
v = 4-16  MHz.  Above  60°K  T^  is  dominated  by  quadrupolar  interactions  while 
below  this  temperature  the  paramagnetic  Te  are  the  dominant  mechanism.  There 
appear  to  be  minima  in  T^  vs  T at  8MHz  but  more  data  are  needed  to  confirm  this. 
The  concentration  of  Te  in  this  sample  is  not  known,  as  yet.  Relaxation  effects 
have  also  been  seen  in  Cr-  and  Zn-doped  GaAs  ingots.  The  results  of  these 
studies  have  not  been  fully  analyzed  but  it  is  clear  that  paramagnetic  impurities 


in  the  ppm  range  in  GaAs  can  be  studied  investigated  by  NMR  techniques. 
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4.0  ELECTRICAL  PROPERTIES  OF  CuInS? 


4 . 1 Introduction 

The  ternary  I-III-VI  compounds  have  recently  received  much  attention 
because,  unlike  their  binary  II-VI  analogues,  they  can  often  be  made  usefully 
both  n-  and  p-type."^  ^ The  largest-bandgap  I-III-VI  compound  having  this 
property  is  CuInS2  (Eq  = 1,5  eV  at  room  temperature)  which  has  evoked  technol- 
ogical interest  both  as  an  electroluminescent  device^  (note  that  the  bandgap 

is  close  to  that  of  GaAs),  and  as  a photovoltaic  device  with  a high  theoretical 

49 

solar-cell  efficiency.  In  spite  of  the  importance  of  its  electrical  properties 
there  have  been  no  studies,  to  our  knowledge,  of  their  temperature  dependences, 
with  subsequent  identification  of  donor  and  acceptor  energies,  and  hole  and 
electron  scattering  mechanisms.  We  have  attempted  to  fill  this  gap  by  producing 
both  n-  and  p-type  CuInS^ , including  some  oriented  single  crystals,  and  then 
measuring  the  hole  and  electron  concentrations  and  mobilities  as  functions  of 
temperature . 

4.2  Experimental  Results 

The  CuInS^  crystals  were  grown  from  the  melt  by  sealing  stoichiometric 

quantities  of  Cu,  In,  and  S in  an  evacuated  quartz  tube  and  heating  to  1050, 

1125,  or  1175°C  for  about  1 day,  and  then  cooling  at  a rate  of  about  2-5°C/hr. 

Single  crystals  of  approximate  dimensions  5mm  x 5mm  x 1mm  could  easily  be 

found,  especially  in  the  1175°C  growth.  Good  p-type  behavior  (p  ~ lC!-cm, 

Op  ~ 20  cm^/V  sec)  could  then  be  obtained  by  annealing  at  550°C  in  a sulfur 

2 

overpressure,  and  good  n-type  characteristics  (p  ~ lQ-cm,  ~ 100-200  cm  /V  sec), 
by  annealing  in  indium  and  excess  CuInS2  powder  at  725-800°C. 

In  Fig.  4-1  we  present  Op  vs.  T for  a single  crystal  grown  at  1175°C 
and  then  annealed  at  550°C  in  a sulfur  overpressure.  The  circles  and  squares 
represent  two  pieces  of  the  same  crystal,  and  the  solid  lines  represent  theoretical 
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Fig.  4-1.  The  hole  mobility  vs.  temperature  Cor  a 
single  crystal  of  CuInS^. 
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expressions  for  ionized-impurity  scattering  (tij.)  , polar-optical  scattering 

(y  )^,  and  acoustical  scattering  combined  with  non-polar-optical  scattering 
62 

(u  ).  These  three  contributions  to  the  scattering  are  combined  by 
anpo 

Matniessen's  approximation  to  obtain  the  sblid  line  through  the  experimental 

points.  The  only  undetermined  parameter  is  the  hole  effective  mass,  m*^, 

which  was  taken  to  be  1.3  m to  fit  the  data. 

e 

3/2 

An  acceptor  activation  energy  of  0.145  eV  is  obtained  from  a p/T  vs 

T ^ plot  (not  shown)  and  by  using  our  measured  m*  = 1.3  m we  also  obtain 

P e 

17  -3  17  -3 

Np  - 1.6  x 10  cm  , - 3.4  x 10  cm  It  is  interesting  that  the  calculated 

"hydrogenic"  acceptor  level  is  0.146  eV,  but  for  a level  this  deep  the  agree- 
ment is  probably  fortuitous. 

In  Fig.  4-2  we  present  y^  vs  T for  samples  annealed  in  1^  at  various 

temperatures.  Here  the  solid  lines  are  simply  drawn  through  the  experimental 

points.  One  salient  feature  of  these  plots  is  the  dominance  of  ionized- 

3/2 

impurity  scattering,  as  evidenced  by  the  T dependence  (for  two  of  the  samples) 

and  the  relatively  high  temperatures  for  the  y^  vs  T maxima.  At  the  lower 

annealing  temperatures  (725  and  745°C)  y^  is  a very  steep  function  of  T, 

63  3/2  “1 

suggesting  high  compensation.  The  n/T  vs  T data  (<iot  shown)  give  a 
deep  donor  at  0.35  eV,  and  a shallow  donor  whose  energy  we  could  not  measure 
accurately. 

4 . 3 Conclusions 

The  conclusions  from  this  study, as  supported  by  the  above  data  and 

other  data  and  calculations  from  kef.  64,  are  as  follows. 

CuInS?  crystals,  grown  from  the  melt  at  1175°C,  have  very  low  mobility, 

18  _3 

evidently  due  to  high  concentrations  (>  10  cm  ' ) of  native  donors  and  acceptors 
which  are  closely  compensated.  However,  annedls  in  S at  550°C,  or  in  In  at  800°C, 
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produce  good  p-  and  n-type  conductivity,  respectively.  An  acceptor  level  at 

0.15  eV  and  a donor  level  at  0.35  eV  have  been  identified,  but  there  also  is 

a shallow  donor  level  and  probably  a deeper  acceptor  level.  The  hole-mobility 

data  are  well  fitted  below  300°K  with  an  effective  mass  m*  - 1.3  if  all 

P 

undetermined  parameters  in  the  mobility  equations  are  taken  as  averages  of 

those  same  parameters  in  CdS  and  ZnS,  the  II-VI  analogues.  This  high  value 

of  m*  can  be  explained  by  using  simple  two-band  k*p  perturbation  theory 

if  the  valence  bands  have  about  65%  d-character.  However,  interactions  with 

other  bands,  which  will  also  affect  m*  have  not  been  taken  into  account. 

P 

Above  300°K  begins  to  fall  below  the  expected  temperature  dependence, 
evidently  due  to  multiband  hole  conduction  and  interband  scattering  arising 
from  a pair  of  degenerate  bands  about  0.02  eV  below  the  upper  (spin-orbit  split) 


band. 


5.0  ELECTRICAL  PROPERTIES  OF  GaAs  SUBSTRATE  MATERIALS 


5.1  Introduction 

Although  semiconductor-device  technology  is  still  centered  on 

silicon,  the  growing  need  for  higher-frequency  (>  10  GHz)  and  higner-temper- 

ature  devices  has  stimulated  work  on  compound  semiconductor  materials  with 

higher  mobilities  and  larger  bandgaps  than  those  of  Si.  By  far  the  most 

promising  of  the  compound  semiconductors  at  this  stage  is  GaAs,  with  a band- 

2 

gap  of  1.43  eV  and  an  electron  mobility  as  high  as  8,500  cm  /V  sec,  at  room 

temperature.  Among  the  other  advantages  of  GaAs  are  a direct  bandgap,  easy 

formation  of  p-n  homojunctions,  good  solid  solution  formation  with  other  III-V 

compounds  (to  obtain  a wide  range  of  bandgaps) , and  good  solubility  of  deep 

donor  and  acceptor  centers  so  that  very  high-resistivity  bulk  material 
9 

(>  10  fi-cm)  can  be  grown.  This  latter  property  makes  the  future  of  large-scale 
GaAs  integrated  cirucits  very  bright,  since  the  necessary  epitaxial  layers 
do  not  have  the  disadvantages  of  being  grown  on  a foreign  material  or  else 
isolated  by  a p-n  junction,  as  is  the  case  for  Si  integrated  circuits. 

Unfortunately,  materials  problems  are  holding  back  the  development  of 
many  GaAs  devices,  especially  the  microwave  devices.  It  is  perhaps  easy 
to  understand  why  industrial  sources  have  not  seen  fit  to  invest  large  sums 
of  money  to  improve  the  material  purity  and  stoichiometry  since  the  main 
commercial  application  at  present,  and  in  the  forseeable  future,  is  for 
light  emitting  diodes  (LED),  which  do  not  require  especially  good  material. 

Thus,  during  the  last  two  years  of  this  contract  period  we  have  carried  out 
a GaAs  materials  characterization  program,  using  electrical,  optical,  and 
magnetic-resonance  techniques.  The  goal  has  been  to  identify  and  help  solve 
material  problems  in  order  to  produce  better  microwave  devices. 


We  have  concentrated  on  substrates  for  reasons  outlined  below.  It  is 


safe  to  say  that  virtually  every  GaAs  microwave  vendor  believes  that  "poor" 

substrates  are  a major,  if  not  the  major,  cause  of  poor  device  performance. 

Recently,  workers  at  several  major  industrial  and  university  laboratories  have 

experimentally  verified  this  assertion.  For  example,  one  scientist  at  an 

industrial  laboratory  reports  that  of  five  substrate  (GaAsrCr  and  GaAsrO) 

brands  that  he  has  used,  one  of  these  brands  is  the  best  if  a field  effect 

transistor  (FET)  is  fabricated  directly  onto  the  substrate  surface,  and  the 

worst  if  an  epibuffer  layer  is  used  for  isolation.  The  reasons  are  unclear. 

As  another  example,  a quite  definitive  study  has  been  performed  by  J.  Barrera, ^ 

at  Hewlett-Packard,  in  which  he  closely  correlates  poor  FET  performance  with 

the  formation  of  a conduting  layer  on  the  surfaces  of  certain  brands  of 

substrates,  under  liquid-phase-epitaxy  (LPE)  growth  conditions.  Barrera  and 

3 

others  have  also  noted  gross  nonuniformities  in  dislocation  counts  (from  10 
5 2 

to  10  /cm  over  a set  of  ingot  slices),  impurity  clustering  at  dislocations, 

Cr  compound  precipitates,  and  other  defects.  Furthermore,  light  irradiation, 
especially  near  l.ly,  can  adversely  affect  FET  performance.  Barrera  has 
suggested  a method,  of  "eliminating"  poor  substrates  hy  actually  measuring  the 
surface  layer  resistance  after  a simulated  epicycle,  and  rejecting  those 
substrates  that  end  up  with  greater  than  a specified  surface  conductivity. 

We  have  tried  to  develop  some  room-temperature  methods  of  quickly  and  accurately 
determining  the  suitability  of  a particular  substrate  material  for  a 
particular  device  to  be  fabricated  under  particular  growth  conditions.  This 
goal  has  not  yet  been  reached,  of  course,  but  we  believe  that  we  have  made 
substantial  progress. 


5 . 2 Hall-Effect  and  Magnetoresistance  Measurements 

Probably  the  singly  most  applied  characterization  and  research  technique 

for  semiconductors  is  the  Hall  effect.  Briefly  stated,  a magnetic  field 

(B  = 3z)  applied  at  right  angles  to  a current  (I  = lx)  in  a crystal  will 
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produce  a voltage  (V  ) perpendicular  to  I and  B.  Generally,  a voltage  (V„) 

H C 

parallel  to  I is  also  measured  so  that  the  resistivity  (p) , mobility  (y  ), 

P 

and  carrier  concentration  (n)  may  be  deduced  as  follows  (practical  units} : 


P = (^r)  ~y  (Q-cm) 


8 V 2 

= I°_  £ _H  ( cm 

Mn  B y Vc  ''V  sec' 


, ,,  in10  IB  -3 

n = 6.25  x 10 cm 

VHz 


where  x,  y,  z are,  respectively,  the  crystal  length  (or  points  between  which 
is  measured),  width,  and  thickness  in  centimeters,  and  I is  measured  in 
amperes,  Vc  and  V in  volts,  and  B in  gauss.  Here  an  n-type  sample  is  assumed, 
but  the  generalization  is  obvious. 

When  the  conductivity  is  mainly  of  one  type,  i.e.,  either  n or  p,  the 
temperature  dependence  of  the  above  quantities  can  be  used  to  deduce  donor 
and  acceptor  concentrations  and  energy  levels  in  a nondegenerate  crystal, 
according  to  well-known  relationships.^  Unfortunately,  the  carrier  concen- 
trations and  mobilities  in  GaAs:Cr,  the  most  useful  high-resistivity  substrate, 


are  not  given  by  Equations  5-2  and  5-3,  and  herein  lies  a common  error  in  the 
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analyses  of  these  crystals.  The  reason  is  that  the  conductivity  is  mixed; 
i.e.,  both  hSles  and  electrons  are  contributing.  Although  this  fact  has 
long  been  suspected.  It  was  first  proved  by  our  research  during  this  contract. 
The  Fall  coefficient,  extrapolated  to  zero  B,  is  given  by  (if  r^  * r “ 1) , 


a p 

PU  “ - ny 


e (pu  + ny  )" 
P n 


V z 3 

= io8  JL  <£!_> 

IB  coul 


and  the  Hall  "mobility"  is 

2 2 
pU'  - nu 


e(pu  + nu  ) 
P n 


108  x , cm2 

B y \ 'V  sec^ 


where  o = l/o.  Thus,  since  y^/p  > 1,  it  is  evident  that  if  an  acceptor  is 

dominant,  the  apparent  mobility,  ffcj,  could  be  either  positive  (p-type)  or 

negative  (n-type) , depending  upon  the  value  of  n/p.  It  is  also  evident  that 

we  cannot  obtain  four  quantities,  n,  p,  U^,  and  y , from  two  measurements, 

V and  V„.  The  resolution  of  this  problem  lies  in  measuring  R and  o(L.e.,V„ 
H v/  * H 
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and  V^,  with  I constant)  as  a function  of  B.  For  example,  it  is  known  that 
R O + R O + R R 0 * (R  + R)B 

R - -3-s — 2P . p.  — hr1-  (5-6) 

(a  + a ) + a a (R  + R ) B 

n p n p n p 


where  a = ney  , U = neu  , R = r (ne)  \ and  R = r (pe)  ^ . A similar 
r.  n p p nn  ’ ppr 

equation  can  be  written  for  0 (3) . During  this  contract  period  we  have 


formulated  a systematic  method  of  obtaining  u , u , n and  p from  the  measurements. 

n p , 

The  results  are: 

A . 2 fltt.  +-  g-i  T = iw!  (5_7) 

1 - T/0  * S (5  7' 

P 


b 1 h mTa  (A  ^ ~ 4} 

P 
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_ n 1 + ab8 
c ” p “ b(ab  + 6) 


% " 


1 + 8 


-1 


(-R  0 ) 

O O 


1 - (ab) 


-1 


y **  y /b 

P n 


ey  (c+b_1) 
n 


n = cp,  = /np 


(5-9) 

(5-10) 

(5-11) 

(5-12) 


where  8 Is  the  x-lntercept  of  an  R versus  Ap/p  (magnetoresistance)  plot,  S 
2 

is  the  slope  of  1/B  versus  Po/A p plot,  and  r^  are,  respectively,  the  Hall 

factors  for  electrons  and  holes,  a = r /r  , n.  is  the  intrinsic  concentration, 

’ n p i 

and  the  subscript  zero  denotes  a measurement  at  B - 0,  The  quantities  r , 

n 

r^,  and  a are  close  to  unity  and  can  be  taken  as  such  if  they  are.  not  explicitly 
known.  Figures  5-1  and  5-2  illustrate  the  process  for  sample  A. 

More  can  be  learned  from  the  temperature  dependence  of  these  quantities, 
especially  y^  and  p,  shown  in  Fig.  5-3.  For  thermally  activated  holes,  well 
below  the  exhaustion  region  (p  <<  N - N ) , p should  be  given  by 


P - <jfT 


1)  S.  e-VkT 


(5-13) 


D 


where  and  N'  are  the  acceptor  (chromium)  and  net  donor  concentrations, 

3/2  3 

respectively,  g is  the  degeneracy  factor,  and  N = 2 (277m  *kT;  /h  is  the 
V v p 

valence-band  density  of  states.  Here  (=N^S)  is  the  concentration  of 

all  acceptors  lying  below  the  chromium  level.  The  quantities  and  g^  are 

known  reasonably  well,  and  EA  can  be  determined  from  a plot  of  In  p versus 
“1  s 

T , thus  leaving  N^,  N^,  and  N as  unknovrns.  The  various  relationships,  to  be 
discussed  later,  are  illustrated  in  Fig.  5-4.  For  most  samples  we  can  ignore 
the  upper  Cr  level  shown  in  Fig.  5-4  and  we  will  also  assume  that  if  oxygen  is 
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SAMPLE  A 
370°K 


I / B2  x I08  ( gauss2  ) 


( cm 


present  its  levels  are  mostly  empty 


Further  information  may  be  obtained  from  the  electron  mobility,  which 


in  Mathiessen’s  approximation  is  given  by  U 


is  the  contribution  of  lattice  modes  to  the  electron  scattering,  y is  the 


contribution  of  neutral  C r centers  (the  Cr  level  is  only  partially  filled 


he  contribution  of  ionized  centers,  including 


with  electrons') , and  U 


all  donors  lying  above  the  Cr  level,  all  acceptors  lying  below  the  Cr  level 


and  the  ionized  fraction  of  Or  centers.  By  using  the  Conwell-Weisskopf 


relationship  for  y (the  mere  familiar  Brooks-Herring  formula  is  not 


caled"  Erginsoy  relationship  for  y , we 


applicable  here),  and  a 


write,  at  room  temperature, 


are  measured  in  ppm,  and  v;e  have  used  the  commonly 


that  there  are  three  unknowns,  N 


, and  N, 


5-13  and  5-14  reveal 


while,  unfortunately,  there  are  only  two  measured  quantities,  p and  y 


problem  is  not  serious  in  Eq.  5-14,  since  the  second  term  is  about  a factor 


10  smaller  than  the  third,  but  must  be  dealt  with  in  Eq.  5-13.  As  a first 


approximation  we  can  use  the  typical  compensation  ratio  for  undoped  hulk  GaAs 


C.5.  Then,  using  Equations  5--7  through  5-14,  it  is  possible  to 


estimate  U (che  chromium  concentration)  and  N_  (often  silicon  and  oxygen) 


To  our  knowledge  we  are  tha  first  group  to  obtain  such  data,  although  many 


groups  have  previously  investigated  the  electrical  properties  of  GaAs.’Cr.  We 


illustrate  the  method  below 


3/2  -1 

For  one  sample  (A),  an  p/T  versus  T plot  gives  F.  - 0.69  eV  from  the 

valence  band  (This  is  a T = 0 value  if  E (T)  = E (0)  + aT) . Then  p and  u 

A A p 

are  measured  only  at  room  temperature  for  other  samples  (B,  C,  D,  and  E)  and 

and  are  calculated.  The  results  are  shown  in  Table  5-1  along  with 

spark-source  mass  spectrographic  (SSMS)  data  obtained  on  Samples  B,  D,  and 

E.  The  SSMS  data  are  considered  accurate  to  within  only  a factor  three.  It 

is  clear  that  for  even  our  simple  model,  which  can  be  improved,  the  results 

are  consistent  with  the  SSMS  data.  The  data  for  each  sample  may  be  accumulated 

and  analyzed  with  a hand  calculator  in  less  than  1-1/2  hours  from  receipt  of 

the  sample  and  thus  the  cost  is  less  than  one-tenth  that  of  SSMS.  Of  course, 

SSMS  gives  a "total"  spectrum  of  elements,  whereas  our  method  gives  only  Cr 

and  the  total  donor  content.  In  this  illustration,  the  data  are  of  particular 

interest  because  they  show  that  the  bottom  ends  of  the  boules  (D  and  E)  contain 

higher  Cr  concentrations  than  the  top  ends  (B  and  C) . 

Because  both  epigrowth  and  post-ion-implantation  annealing  require  high 

temperatures  (700-900°C),  it  is  important  to  know  whether  or  not  the  GaAs:Cr 

substrates  retain  their  high  resistivity  after  exposure  to  such  temperatures. 

Several  groups,  as  mentioned  earlier,  have  already  investigated  this  problem, 

at  least  from  a practical  point  of  view.^  To  check,  we  heated  Sample  A for 

10  minutes  at  600°C  in  an  Ar  atmosphere  and  the  resistivity  dropped  from 
8 7 

8 x 10  to  2 x 10  fi-cm.  Upon  removal  of  about  20p,  the  resistivity  returned 
to  its  initial  high  value.  In  an  independent  study  we  have  shown  that  heating 
"undoped"  GaAs  ingots  in  vacuum  for  24  hours  at  600  or  700°C  produces  a p-type 
sample  with  an  acceptor  level  at  0.14  eV  from  the  valence  band.  We  believe 
that  this  center  is  associated  with  a Ga  vacancy,  although  further  work  needs 


to  be  done  to  verify  this. 


Table  5-1  The  relsitivity  (p) , electron  mobility  (u  ),  hole  concentration 
(p),  donor  (Np)  , and  acceptor  (%)  concentrations  (calculated 
from  Hall-effect  measurements),  and  Si(Nsi),  0(N~)  and  Cr(Ncr) 
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The  apparatus  necessary  to  perform  the  Hall-effect  and  magnetoresistance 
measurements  on  high-resistivity  samples  has  been  completed  and  working  for 
about  1-1/2  years.  It  makes  use  of  electrometers,  operated  in  a unity  gain 
mode,  guarding  each  voltage  contact.  Triaxial  cable  is  employed  with  the 

inner  shield  of  each  cable  held  at  its  respective  center-lead  potential  to 

1 3 

effectively  reduce  cable  capacitances.  The  system  can  measure  10  ohms  with 
1 percent  accuracy,  limited  by  a time  constant  of  only  10  seconds  at  this 
resistance.  Our  He  exchange-gas  dewar  can  operate  between  5.5  and  600°K 
with  an  efficient  cool-down  period  of  45  minutes. 


5.3  Photomagnetoelectric  Measurements 


The  photomagnetoelectric  effect  is  essentially  the  Hall  effect  of  a 
— ► 

diffusion  current.  If  light  I is  absorbed  near  one  surface,  the  created 

o 

holes  and  electrons  will  diffuse  toward  the  opposite  face.  If  a magnetic  field 
- ► -► 

B is  applied  perpendicular  to  I , the  carriers  will  experience  a force  perpendicular 

-f  + 

to  both  B and  I and  we  can  measure  a short-circuit  current  I„,„,,  or  an  open- 
o PME  K 

oc 

circuit  voltage  in  this  perpendicular  direction.  These  quantities  depend, 

of  course,  upon  the  intensity  of  the  light  and  the  conversion  efficiency, 

both  hard  quantities  to  determine  accurately.  One  way  to  circumvent  this 

problem  is  to  measure  the  photoconductivity  under  the  same  conditions.  Then, 
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it  can  be  shown  that 


xBI  . . 

E£  /L_  _ -i_ 

v 5 H 

PME  V V 

pc  pc 


— ) a — (1  + — i— ) 
d ' L U * 


(5-15) 


1+Sx/L 


(5-16) 


Here  I is  the  photoconductivity  current,  held  constant,  and  are 

pc  y ’ pc  pc 

the  PC  voltages  in  the  light  and  dark,  respectively,  between  two  points  a 
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the  carrier  lifetime,  and  L the  diffusion  length 


distance  x apart,  I 


assumed  equal  for  holes  and  electrons,  a is  the  absorption  constant,  and 


S is  the  front-surface  recombination  velocity.  We  are  assuming  here  that 


od  >>  1,  where  d is  the  sample  thickness  in  the  light  direction 


approaches  a constant,  V. 


it  can 


be  shown  that 


the  ambipolar  diffusion 


The  quantities  L and  T are  related  by  L 


onstant,  is  given  by 


relationships.  Note  that  for  an  n-type  sample,D  = D , and  for  p-type,  D = D ; i.e 


ill:  is  a minority  carrier  effect.  From  plots  of  either  (I  /I 


versus  ct,  we  can  obtain  L (or  x)  and  S 


In  case  trapping  centers  are  present,  x 4 x 4 X and  x 4 x. 


be  shown  that  Eq.  5-15  can  be  revised  by  letting  x -*■  X 


where  x. 


the  measurement  of  all  of  the  quantities  I. 


hy  Equations  5-15  and  5-17,  and  thus  x and  x 


The  apparatus  we  have  been  using  for  the  investigation  of  the  PC  and  PM£ 


effects  includes  the  electrical-measurement  system,  described  in  the  previous 


paragraphs,  along  with  a high-intensity  B&L  monochromator  with  a lens.  The 


monochromator  is  capable  of  providing  a constant  1 x 10  photons/cm  sec  over 


the  wavelength  range  0.7  to  1.6m,  with  a resolution  of  about  0.02  eV 


L 


For  Sample  A (Table  5-1)  we  nave  performed  detailed  measurements  of 

I , and  V as  a function  of  wavelength  (or  a).  The  results  are: 

PME  pc  pc 

-8  4 

T - T - 1.4  x 10  sec  and  S - 5 x 10  cia/soc.  We  have  also  measured  i 
n p 

-g 

for  a total  of  ]0  other  GaAs:Cr  samples  and  have  found  values  from  5 x 10 
to  3 x 10  ^ sec.  Some  of  the  samples  exhibit  ttapping  effects.  Considering 
this  range  of  T values  it  seems  clear  that  the  PME  Technique  is  sensitive  to 
the  differing  properties  of  the  substrai  e materials  3nd  should  be  valuable 
as  a characterization  tool. 

5 . 4 Photo-Hall,  Photomagnetorcsistancc . and  Photoconductivity 
Measurements 

In  the  presence  of  light  irradiation,  new  equilibrium  concentrations 
of  holes  and  electrons  will  be  established,  and  consequently  all  of  the 


quantities  discussed  in  Sec.  5.2  may  fce  affected.  In  chia  section  we  show  how 
the  various  effects  can  be  isolated  and  we  present  a model  for  the  deep 
levels  in  GaAs  which  is  consistent  with  our  present  results.  However,  more 


data  will  be  necessary  to  confirm  the  model. 

Most  of  the  GaAsrCr  wafers  which  wa  have  analyzed  yield  room-temperature 
photoconductivity  (PC)  and  photo-Hall  (PH  or  Ra)  curves  similar  to  those  shown 
in  Fig.  3-5.  For  this  sample  the  PC  increases  from  a threshold  at  about 
0.72  eV,  with  perhaps  anot.her  threshold  existing  between  0.8  and  0.9  e\r.  The:, 
the  intrinsic  (above-band-gap)  PC  takes  over  above  1 t eV.  In  a mixed- 
conductivity  system  like  this,  an  increase  in  PC  could  ba  due  to  increases  ir. 
p,  n,  Up,  or  u^,  or  even  some  complicated  combination  of  increases  and 
decreases  of  these  quantities.  Usually,  it  is  assumed  that  this  threshold  is 
due  to  electrons  being  excited  from  a deep  impurity  level  to  the  conduction  band. 


but  in  most  cases  this  is  wrong.  The  Ro  data  help  to  clarify  the  situation 
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photoconductivity  and  photo-Hall  mobilit' 
light  energy  for  GaAs:Cr  Sample  F at  300' 
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somewhat  by  showing  that  this  sample  is  more  p-type  (i.e.,  R a is  less  negative) 
in  above-threshold  light  than  in  the  dark.  It  is  doubtful  that  this  could 
occur  if  the  threshold  corresponds  to  electrons  being  excited  to  the  conduction 
band. 


Amother  sample,  with  similar  Rct  data  (not  shown),  is  illustrated  in 

Fig.  5-6.  Here  two  thresholds  at  0.72  and  0.81  eV  are  clearly  seen.  Also, 

p and  p/n  (not  shown)  are  rapidly  increasing,  and  is  decreasing  as  the 

light  energy  is  increased.  These  results  are  consistent  with  electrons  being 

excited  from  the  valence  band  to  acceptor  levels,  thus  increasing  p and 

decreasing  Pr,  since  the  acceptors  are  becoming  ionized.  For  this  sample, 

3/2  -1 

p/T  vs  T data  give  an  acceptor  level  at  0.83  eV  above  the  valence  band. 

Even  though  this  is  a T = 0 value  of  E^,  still  it  is  natural  to  identify  it 

with  the  0.81  eV  level  measured  by  the  room  temperature  PC.  On  two  other 
3/2  -1 

GaAs:Cr  samples,  p/T  vs  T measurements  have  yielded  E^  = 0.69-0.72  eV, 

which  probably  corresponds  to  the  PC  value  of  0.72  eV.  If  this  model  is 

correct,  then  the  two  Cr  centers  must  remain  relatively  fixed  to  the  valence 

band  as  a function  of  temperature.  It  is  probable  that  one  of  the  two  levels 
3+  2+ 

is  Cr  on  a Ga  site  (becoming  Cr  after  accepting  on  electron)  while  the 
other  may  be  an  interstitial  or  complex.  Further  speculation  must  await 
more  data. 

A few  samples,  such  as  Sample  B shown  in  Fig.  5-7,  exhibit  differences 
compared  with  those  in  Figures  5-5  and  5-6.  One  difference  is  a lower  threshold 
energy,  about  0.68  eV,  on  the  average.  Another  difference  is  that  Re  is  more 
n-like  at  the  threshold,  than  in  the  dark,  and  n/p  is  increasing.  A third 
difference  is  a relatively  low  cr  at  the  impurity  maximum.  These  features 
are  all  consistent  with  the  presence  of  a significant  amount  of  oxygen,  which 
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E(eV) - 

Fig.  5-6.  The  photoconductivity,  electron  mobility  (un) * anc*  hoL 
concentration  (p)  vs.  light  energy  for  GaAs  Sample  G 
at  300°K. 


I 
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is  thought  to  be  a deep  donor  and  recombination  center.  In  fact,  the  SSMS 
analysis  of  this  sample  showed  about  7 ppm  0 and  only  about  4 ppm  Cr. 

Although  the  SSMi  analysis  is  not  accurate  enough  to  conclude  that  there 
is  more  0 than  Cr,  still  the  data  seem  to  show  the  influence  of  a deep 
donor.  In  another  sample  (not  shown)  doped  with  oxygen  alone  (no  Cr) , 

3/2 

the  PC  data  showed  an  electron  threshold  at  about  0.69  eV,  and  the  n/T 
vs  T 1 data  showed  a donor  level  lying  about  0.75  eV  from  the  conduction 
band  (at  T = 0) . Thus,  our  present  data  are  consistent  with  oxygen  producing 
a deep  donor,  remaining  just  above  mid-gap  as  a function  of  temperature. 

However,  it  must  be  emphasized  that  more  data  are  necessary  to  confirm  the 
location  of  this  level  as  well  as  the  two  Cr  levels.  Our  present  model  for  these 
centers,  at  T * 0,  is  depicted  in  Fig.  5-4. 

Thus,  it  appears  that  the  techniques  described  in  this  section  are  quite 
valuable  for  the  characterization  of  GaAs:Cr.  It  will  be  necessary  to  confirm 
the  explanations  for  the  various  PC  thresholds,  peaks,  and  valleys  hut,  once 
this  is  done,  the  next  step  should  be  to  obtain  correlations  between  these 
results  and  device-reliability  results.  Then  it  should  be  possible  to  suggest 
improvements  in  the  materials-growth  techniques. 


6.0  THERMAL  PROPERTIES  OF  II-VI  COMPOUNDS  AND  RELATED  MATERIALS 


6.1  Introduction 

Studies  of  the  thermal  conductivity  of  II-VI  type  crystals  are 

important  in  their  own  right  since  this  type  of  information  is  most  necessary 

for  the  engineer  concerned  with  device  applications,  especially  in  a variable 

temperature  environment.  Additionally,  the  use  of  low-temperature  thermal 

conductivity  changes  as  a probe  in  studies  of  the  character  of  lattice 

imperfections  induced  by  radiation  has  been  well-known  for  several  decades. 

However,  partly  because  of  the  nonavailability  of  good  crystal  specimens, 

such  research  on  almost  all  of  the  II-VI  compounds  has  been  neglected. 

Another  reason  for  this  neglect  is,  of  course,  the  fact  that  much  of  the 

early  work  was  confined  to  the  elemental  semiconductors  such  as  Ge  and  Si 

because  of  the  less  complicated  associated  problems.  Many  of  the  conclusions 

of  these  preliminary  studies  indicate  the  advantages  of  extending  this  method 

to  other  classes  of  materials.  It  has  been  predicted  with  confidence  that 

many  interesting  and  important  findings  will  be  made  in  compound  semiconductors, 

74 

in  particular,  in  the  II-VI  materials. 

Thermal  conductivity  theory  indicates  that  this  property  is  particularly 
sensitive  to  impurities  and  defects  in  the  temperature  region  where  the 
thermal  conductivity  is  a maximum.  Although  in  some  cases  the  maximum  is 
near  room  temperature,  for  most  dielectric  solids  the  important  temperature 
region  is  below  50°K.  Thus,  low-temperature  techniques  are  desirable. 

Radiation  produces  an  additional  thermal  resistance  and  to  a first 
approximation  this  addition  can  be  regarded  as  a quantity  which  is  additive; 
to  the  normal  lattice  resistance.  In  quartz,  for  example,  the  additional 
thermal  resistance  has  a minimum  at  approximately  15°K.  Above  15°K.  it 
increases  directly  with  the  temperature  and  it  has  been  suggested  that  in 
this  region  the  resistance  is  caused  by  the  presence  of  single  lattice  defects. 
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Below  15°K  the  radiation-induced  increase  in  the  thermal  resistance  increases 
more  rapidly  than  i/T  and  it  is  believed  that  the  cause  is  due  to  the  presence 
of  large  defects  at  Jower  temperatures. 

The  effectiveness  of  the  thermal  conductivity  tool  results  from  the 
ability  of  the  different  types  of  lattice  defects  to  produce  correspond- 
ing temperature-dependent  differences  in  the  phonon-scattering  probabilities. 
For  example,  the  several  well-known  behaviors  and  defects  may  be  tabulated 
as  follows: 


Table  6~1  Thermal  resistance  behaviors  of  various 

defect  types 


Lattice  Defect 

Point  defect 

Dislocations,  elongated  disorders, 
etc.,  linear  regions 

Two-dimensional  boundaries 


Thermal  Resistance  Behavior 

Increase  is  proportional  to 
temperature,  T 

Increase  is  inversely  propor- 
tional to  T 

Increase  is  inversely  propor- 
tional to  T^ 


Since  any  physical  or  chemical  defect  can  be  a source  of  phonon  scattering, 
those  materials  which  have  the  fewest  defects  are  best  suited  for  extensive 
thermal  conductivity  measurements.  In  recent  years  small  single  crystals  of 
some  of  the  II-VI  compounds  have  been  prepared  in  a high  purity  state. 

Feedback  from  thermal  conductivity  measurements  on  these  crystals  has  aided 
in  the  development  of  better  growth  tec’  ’’ques.  At  the  same  time  such 
measurements  are  importan:  because  of  ability  to  form  a basis  for 

comparing  thermal  conductivity  results  on  irradiated  crystals  and,  thereby, 
identify  radiation-induced  defects.  Such  results  give  information  on  the 
nature  and  properties  of  the  interstitial,  the  role  of  impurities,  and  perhaps 
help  to  identify  the  temperature  at  which  long-distance  defect  migration  begins 
to  take  place. 
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Furthermore,  thermal  conductivity  measurements  can  provide  much  insight 


in  differentiating  between  various  possible  processes.  For 


measuring  the  change  in  thermal  resistance  along  different  crystallographic 


directions,  before  and  after  irradiation  and  after  temperature  pulses  of 


an  isochronal  anneal,  the  presence  of  anisotropic  defects  as  well  as  the 


such  defects  are  anisotropic  phonon-scattering  centers  which  produce  an 


anisotropic  increase  in  the  thermal  resistance  of  the  sample.  The  anihilation 


of  interstitials  at  vacancy  sites  can  also  be  detected  as  a decrease  in  the 


thermal  resistance.  Thermally-induced  defect  clustering  would  produce,  on 


the  other  hand,  an  increase  in  thermal  resistance.  Such  defects  might 


have  a preferential  orientation  in  a wurtzite  structure  or  might  be  prefer 


entially  oriented  by  an  externally  applied  strain  in  the  case  of  zinc-blende 


lattices.  Thermal  conductivity  measurements  are,  therefore,  vital  to  a 


complete  understanding  of  the  defect  structures  produced  as  a result  of 


bombardment  and  annealing.  By  combining  this  technique  with  NMB,  luminescence 


optical  absorption,  and  electrical  measurements,  a much  better  understanding 


6 . 2 Thermal  Properties  Equipment 


The  equipment  with  which  the  absolute  thermal  conductivity  measure' 


ments  have  been  made  has  undergone  modif ications 


(some  of 


it  to  be  used  to  make  measurements  of  magnetic  susceptibility 


these  measurements  have  been  carried  out  on  samples  at  temperatures  below 


and  the.  specific  heat  of  various  II-VI 


compounds 


Preliminary  thermal  conductivity  experiments  have  also 


been  made  on  CdS  with  different  concentrations  of  magnetic  impurity.  In 

addition,  the  experimental  apparatus  has  been  used  to  obtain  absolute 
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thermal  conductivity  results  for  thermal  insulators  such  as  nylon  and 
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thermal  bonding  agents,  such  as  Cry-Con  and  the  Apiezon  greases. 

These  results  have  been  especially  useful  in  the  research  since  it  is 

necessary  to  know  the  thermal  transport  properties  of  the  technological 

materials  when  they  are  used  to  support  the  experimental  specimens. 

Furthermore,  anomalies  observed  in  the  heat  conduction  of  thermal  eret.ses 
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have  special  application  when  grease  is  used  as  a thermal  bond.  Likewise, 

information  on  methods  of  establishing  low-temperature  mechanical  contacts 
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with  the  II-VI  compounds  has  also  been  useful.  Other  improved  aspects 
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of  the  equipment,  such  as  the  incorporation  of  a superfluid  value  into 
a variable  temperature  dewarx<^,  have  been  useful  in  helping  to  carry  out 
cryogenic  experiments. 

6.3  Thermal  Conductivity  Measurements  on  II-VI  Compounds 


Thermal  conductivity  results  at  room  temperature  on  undoped 
polished  single  crystals  of  ZnO,  ZnS,  ZnSe,  ZnTe,  CdS,  CdSe,  CdTe,  HgS  and 
crystals  with  Li,  In,  Co,  and  Mr.  impurities  have  been  obtained  with  a thermal 
comparator.  Although  the  duration  of  each  measurement  was  approximately  less 
than  10  to  20  seconds,  numerous  readings  were  taken  in  order  to  obtain 
statistical  information  on  the  reliability  of  the  data.  In  most  cases  the 
percent  probable  error  associated  with  the  reproducibility  of  a measurement 
was  found  to  be  in  the  neighborhood  of  1%,  and  the  results  are  in  reasonable 
agreement,  though  slightly  higher  than  previous  experimental  work.  The 

higher  conductivities  are  interpreted  as  due  to  the  more  recent  availability  of 

v , 101 
better  crystals. 
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Because  it  is  not  always  possible  to  use  large  samples,  we  have 
explored  the  use  of  the  comparative  method  of  measuring  the  thermal  conductivity 
of  small  samples.  Plexiglass  disks  are  used  to  hold  various  II-VI  compounds 
mounted  in  woodsmetal.  The  crystals  are  ordinarily  mounted  this  way  for 
polishing  and  no  further  preparation  is  necessary  for  a conductivity 
measurement.  We  have  kept  the  samples  on  a large  copper  plate  in  a room 
where  the  temperature  drifts  and  air  currents  are  minimized.  A sample  is 
lifted  by  tongs  to  the  sample  stage  of  the  comparator.  The  arm  trip  lever 
is  depressed  and  the  probe  assembly  is  released  and  comes  into  contact  with 
the  sample. 

The  probe  is  kept  at  some  temperature  about  20°C  above  room  temperature. 
When  the  sensing  tip  comes  into  contact  with  the  sample,  heat  is  extrs . ted 
from  the  probe  at  a rate  which  is  dependent  on  the  thermal  conductivity 
of  the  sample.  The  response  of  the  probe  produces  a microvolt  reading 
which  is  recorded.  The  entire  operation  takes  about  15  seconds.  Standards 
used  for  calibration  are  armco  iron,  stainless  steel,  titanium  alloy,  quartz, 
glass* and  ebonite. 

After  taking  a reading  on  a sample,  we  have  found  it  convenient  to 
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move  the  sample  slightly  on  the  stage  so  that  the  probe  may  come  into  contact 
with  another  point  of  the  sample's  surface.  Normally  this  second  reading 
should  be  lower  because  of  the  warming  of  the  sample  during  the  first  reading. 

If  the  second  reading  is  larger,  it  may  imply  that  the  contact  of  the  probe 
with  the  surface  is  not  good  and  in  this  case  both  readings  are  discarded. 

With  a properly  polished  surface,  less  than  10%  of  the  readings  need  to  be 
discarded  when  this  technique  is  used. 


Table  6-2  illustrates  some  results  on  zinc  compounds.  Note  that  the 
ratio  of  successful  readings  to  total  reading,  using  the  criteria  just 
mentioned,  varies  from  1 to  about  0.5.  The  ZnSe  mounted  crystal  had  small 
cracks;  and  these  and  other  surface  defects  could  have  decreased  its  ratio. 

The  reprodueiblity  of  the  remaining  readings  is  often  better  than  1 
percent  and  seldom  worse  than  1.5  percent,  as  can  be  seen  in  the  table. 
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Table  6-3  shows  our  results  on  zinc  oxide  along  with  those  of  Koenig, 
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Slack,  ' and  Martin  and  Wolf.  Martin  and  Wolf  appear  to  have  the 
highest  value  at  300°K.  Table  6-3  also  compares  data  for  additional 
zinc  compounds.  It  is  expected  that  polycrystalline  and  hot-pressed  materials 
should  nave  lower  conductivity  values  and  this  is  borne  out  by  the  comparison. 
Martin  and  Wolf  have  reported  their  results  on  ZnO  recently. 

Table  6-4  gives  results  for  some  cadmium  compounds.  There  is  very 
good  agreement  with  CdS.  Our  CdS  crystal  was  analyzed  by  Bell  and  Howell 
and  it  was  found  to  have  6.3  ppm  carbon,  2.9  ppm  H,  and  2.0  ppm  V.  All 
other  impurities  found  were  less  than  1 ppm. 

6.4  Thermal  Conductivity  Results  from  Comparator  Measurements 
on  CdTe  Crystals 

The  work  described  in  this  section  is  partly  representative  of  a 
coupling  effort  of  the  Solid  State  Physics  Laboratory  of  the  Aerospace  Research 
Laboratories  with  the  Laser  Window  group  of  the  Materials  Laboratory.  The 
latter  supplied  the  CdTe  crystals  upon  which  the  measurements  reported 
below  were  measured.  Comparator  results  on  other  II-VI  compounds  have 
been  reported  above. 

The  thermal  conductivity  of  15  different  CdTe  samples  was  measured  and 
the  highest  value  observed,  0.117  i 0.007  Wcm  ^ @ 298  K was  obtained  on 
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Table  6-3.  Thermal  conductivities  of  2n  compounds  at  300°K 
(L  designates  the  linear-heat-flow  method,  and  C,  the  comparator  method) 


I ' 

Crystal  HWcm^K-1) 

Method,  Crystal  Data,  Source 

Reference 

ZnO 

0.31 (polated) 

, Yellow  color,  Rutgers  U. 

102 

Koeing  (1953) 

0.25 

, Grey  color,  Rutgers  U. 

102 

Koenig  (1953)  U 

0.54 

L,  Single,  colorless,  N.J.  Zinc 

103 

Slack  (1966)  J 

• 

0.58 

L,  Single,  hydrothermal,  Ventron 

Martin  6.  Wolf  (1972)104 

0.62 

Li  doped 

104 

Martin  & Wolf  (1972) 

'O-SC^l^d) 

L,  Single,  vapor,  3M 

104 

Martin  & Wolf  (1972)  ” 

. 

0.36+6% 

C,  Single,  hydrothermal,  slight 
yellow  color,  Litton  Ind. 

Present  Results101 

ZnS 

0.265 

L,  Natural  Cubic 

Eucken  & Kuhn  (1928) 106 

0.27 

L,  Single,  vapor,  colorless,  G.E. 

103 

Slack  (1966) 

0.27+7.0 7. 

C,  Single,  melt,  mounted,  E.P. 

Present  Results101 

0.22+9.3% 

C,  Single,  melt,  1 cm  cube,  E.P. 

Present  Results101 

0.19±7.4% 

C,  Single,  melt,  #833,  E.P. 

Present  Results101 

ZnSe 

0.14 

Polycrystalline,  hot-pressed 

Ioffe  & Ioffe  (I960)106 

0. 13@  327K 

Polycrystalline,  hot-pressed 

Ladd  (1963) 107 

0.19 

L,  Single,  vapor,  lemon-yellow,  G.E. 

103 

Slack  (1966)  J 

0.23±6.8% 

C,  Single,  melt,  321  (M7102  BQ),E.P. 

Present  Results101 

0. 1818. 1Z 

C,  Single,  melt,  P-188  'M7103  ER)E.P. 

Present  Results101 

0.15+8.25: 

C,  Polycrystalline,  Coors 

Present  Results101 

0.14±7.4% 

C,  Polycrystalline,  vapor,  Raytheon 

Present  Results101 

ZnTe 

0.17 

1 Oft 

Kelemen  et.al.  (1965) 

0.18 

L,  Single,  vapor,  red-brown,  G.E. 

103 

Slack  (1966) 

0.22+5.9% 

C,  Single,  melt,  E.P. 

Present  Results101 
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Table  6-4.  Thermal  conductivities  of  Cd  compounds  at  300°K 


Crystal 


Method,  Crystal  Data,  Source 


Polycrystalline,  hot-press 


(Cubic) 


L,  Single 


C.  Single,  ARL 


0.08714* 


C,  Single,  Hughes 


0.10+62 


C,  Single,  100  ppm  In,  ARL 


0.20 

L,  Single,  Vapor,  ARL 

0.2 

L,  Single,  Clevite 

0.20+8.4% 

C,  Single,  melt,  semi-elements 

0.043  (cubic) 

, Polycrystalline,  hot-press 

0.063 

, Polycrystalline,  hot-press 

(at  280°K) 

0.12+6.2% 

C,  Single,  Vapor,  ARL 

0.11+7.32 

C,  Single,  Vapor  10  ppm  Li,  ARL 

Holland 

(1964)115  ! 

Vook  (1968)116 

Present 

D !01 

Results 

Ioffe  & 

Ioffe  (I960)106 

Ioffe  & 

Ioffe  (I960)106 

Present 

Q , _ id 

Results 

Present 

Results 

Ioffe  & 

Ioffe  (I960)106 

Slack  Galginaitis  (1964) 117 

Holland 

(1964)115 

Horch  & 

Nieke  (1965)118 

Present 

D . 101 

Results 

Present 

d w 101 

Results 

Present 

Results101 

an  undoped  single  crystal.  Other  single  crystal  samples  which  also  gave 
high  conductivity  values  contained  small  amounts  of  Ga  (20  ppm)  or  In 
(10  ppm) . The  largest  conductivity  values  for  polycrystalline  samples  were 
observed  to  lie  below  0.100  Wcm  The  data  were  obtained  using  a 

commercial  instrument  with  a capability  of  making  a measurement  on  a polished 
sample  surface  within  a duration  of  approximately  10  seconds.  Out  of  a 
total  of  302  measurements,  267  or  91%  were  found  acceptable  (according  to 
criteria  adapted  to  eliminate  erroneous  readings)  and  the  average  reproduci- 
bility of  all  of  these  measurements  was  calculated  to  be  + 0.92%.  The 

average  value  for  the  precents  of  error  for  the  15  samples  was  determined 
119 

to  be  5.3%. 

Table  6-5  describes  15  CdTe  crystals  which  were  measured.  Eleven  of 

the  CdTe  crystals  are  designated  as  single  crystals  in  the  table.  However, 

2 

this  means  that  t ley  contain  large  single  crystal  parts  about  1 cm  in  area. 
The  results  of  thermal  conductivity  measurements  on  these  CdTe  crystals  by 
the  comparator  method  are  shown  in  Table  6-6.  As  is  expected,  single 
crystals, either  pure,  or  compensated  with  small  amounts  of  gallium  or  indium, 
yield  the  highest  conductivities. 

The  criteria  adopted  to  eliminate  poor  readings  can  be  illustrated 
with  sample  number  18.  (See  the  top  of  the  table.)  Only  27  out  of  31  sets 
of  measurements  were  accepted.  The  remaining  4 sets  were  discarded  because 
the  second  reading  in  each  set  was  larger  than  the  first.  And  fhis  was 
taken  as  an  indication  that  the  probe  had  not  come  to  rest  on  a good  spot 
on  the  crystal. 

In  general,  it  might  be  inferred  that  the  condition  of  the  surface  is 
related  to  the  percentage  of  measurements  accepted.  For  all  of  the  CdTe 
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Table  6-5.  Description  of  CdTe  samples  measured 


by  the  thermal-comparator  method 


Sample 

Number 

Material 

Nominal 

Dopant 

Crystal- 

linity 

Size 

Diam. 

(mm) 

Thick 

17 

CdTe 

Undoped 

- Single 

48 

7 

j 

18 

CdTe 

Undoped 

- Single 

54 

10 

23 

33Ge55Sel2As 

Chalcogenide 

Glass  TI-20 

25 

) 

6 

25 

28Ge60Sel2Sb 

Chalcogenide 

Glass  TI-1173 

25 

* 7 

29 

ZnSe 

Undoped 

Poly 

39 

6 

36 

CdTe 

Undoped 

Poly 

39 

4 

38 

CdTe 

lOppmln 

- Single 

27  1 

2 

i 

39 

CdTe 

lOppmln 

- Single 

27 

2 

40 

CdTe 

lOppmln 

- Single 

28 

. 13a 

42 

CdTe 

1500ppmln 

Poly 

49 

4 

43 

CdTe 

Undoped 

- Single 

38 

4 

44 

CdTe 

2 OppmGa 

- Single 

49 

8 

45 

CdTe 

•SOppmGd 

' Single 

49 

n 

56 

CdTe 

1500ppmln 

Poly 

50 

4 

58 

CdTe 

lOOppmCl 

Poly 

51 

10 

68 

ZnSe 

Undoped 

Polyb 

26 

7 

86 

CdTe 

Undoped 

- Single 

47 

87 

CdTe 

2 OppmGa 

' Single 

48 

6 

88 

CdTe 

2 OppmGa 

- Single0 

48 

7 

Or  j gi  r.a. 
Source 


Gould 


1C9 


_ ij10? 

Gould 

Tex.  Inst. 

Tex.  Inst. 
Ill 


no 


110 


Coors 


Coors 


Hughes 

Hughe 


111 


112 


112 


Hughes 

113 

Tyco 


112 


Gould 


109 


Gould109 

Gould109 


Tyco 

Tyco 


113 


113 


Ratytlieon 
,109 


114 


Gould 


Gould 


109 


Gould 


109 


Maximum  dimension  of  edge.  The  crystal  was  in  the  shape  of  & rustrum  ot  a cone 
Irregular  area. 

Etched  surface.  _ 
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crystals  included  here  there  was  a total  of  302  sets  of  measurements  and 
276  or  91.4%  of  these  were  acceptable.  The  average  reproducibility  was 
within  0.92%  and  the  average  error  was  within  5.3%. 


6.5  Thermal  Comparator  Measurements  on  Ion-Implanted  CdS 

Work  on  the  thermal  conductivity  of  ion- implanted  samples  of  CdS 

was  initiated  during  the  report  period  in  order  to  aid  a probe  of  the 

physical  characteristics  of  the  implanted  layers.  This  work  was  carried 

out  with  the  collaboration  of  Y.  S.  Park  and  W.  R.  Woody.  D.  E.  Johnson 

120 

assisted  by  orientating,  cutting  and  polishing  the  crystals.  A paper 
on  the  subject  of  this  section  was  presented  at  the  XII  International 
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Conference  on  Thermal  Conductivity  in  Birmingham,  Alabama.  A summary 
of  the  work  follows. 

Recently,  a growing  interest  has  been  generated  in  type-converting  II-VI 
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compound  semiconductors.  Ion-implantation  techniques  appear  to  be  particu- 
larly attractive  for  II-VI  compound  semiconductors,  because  of  the  difficulty 
of  achieving  type  conversion  and  p-n  junction  formation  by  thermal  diffusion 
and/or  other  conventional  growth  and  doping  techniques.  Furthermore  the 
accomplishment  of  such  transformations  in  the  electrical  characteristics 
when  combined  with  the  typically  large  band  gaps  associated  with  this  group 
of  semiconductors  offers  the  promise  of  optoelectronic  devices  associated  with 
electroluminescence  in  the  visible  range.  However,  in  order  to  use  the  ion- 
implantation  technique  to  its  best  advantage,  many  questions  such  as  the  extent 
of  radiation  damage,  the  lattice  location,  and  the  penetration  depth  of  the 
implanted  ions  remain  to  be  answered.  In  the  past,  thermal  comparators  have 

been  used  for  comparative  thickness  determinations  of  materials  and  surface 
122-125 

deposits.  We  have  used  the  thermal  comparator  to  observe  changes  in 

the  thermal  conductivity  of  the  implanted  layers  in  CdS  in  the  hope  that 


these  measurements  may  contribute  to  the  understanding  and  to  the 
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characterization  of  the  profile  of  the  implants.  A dose  of  singly-ionized 
atomic  nitrogen  (N^)  ions,  measured  to  be  lO^"*  ions/cm^,  was  implanted  at 
90  KeV  into  a polished  CdS  surface,  the  normal  to  which  was  parallel  to  the 
c-axis  of  the  crystal.  Thermal  conductivity  measurements  on  this  crystal 
and  on  three  other  unimplanted,  but  otherwise  identical,  crystals  were  taken 
and  a statistical  summary  of  the  results  produced  the  conclusion  that  the 
thermal  conductivity  of  the  implanted  sample  decreased  by  approximately  25%. 

A second  experiment  produced  a similar  result.  A third  experiment,  in  which 
the  unimplanted  area  adjacent  to  the  implanted  area  was  tested,  resulted  in 
a thermal  conductivity  similar  to  the  unimplanted  samples.  The  results  of  the 
preliminary  measurements  were  encouraging  because  they  demonstrated  the  potential 
of  this  technique.  They  show  that  the  comparator  technique  can  easily  detect 
the  presence  of  an  implant. 

In  order  to  expedite  the  study  of  the  physical  properties  of  implanted 
layers,  thermal  conductivity  measurements  have  been  carried  out  on  single- 
crystal cadmium  sulfide  samples  implanted  along  the  c-axis  direction  with  doses 

+ 14  15  2 

of  90  KeV  N ranging  from  5 x 10  to  5 x 10  ions  per  cm  . The  electrical 

resistivity,  mobility,  and  carrier  concentration  were  found  to  be  5.32  fi-cm, 

2 15  -3 

292  cm'  /Volt-sec,  and  4.03  x 10  cm  , respectively,  typical  of  the  n-type 
base  material.  The  thermal  conductivity  was  measured  by  the  comparator 
method  which,  prior  to  implantation,  gave  a value  of  0.20  W cm  1 K \ in 
good  agreement  with  published  values.  The  decrease  in  thermal  conductivity, 

occurring  vrith  an  increase  in  dose,  was  observed  to  be  equal  to  12.5%  for  a 

1A  2 15  2 

dose  of  5 x 101  ions  per  cm  and  25%  for  a dose  of  10  ions  per  cm  . A 

15  2 

sample  having  a dose  of  5 x 10  ions  per  cm  was  annealed  in  vacuum  for  15 
minutes  at  500°C.  While  the  annealing  removed  to  darkened  region  left  by 
the  implantation,  subsequent  observations  on  the  annealed  sample  disclosed 


that  most  of  the  implant  damage  (as  measured  by  the  remaining  thermal  resis- 
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tance)  was  still  present.  More  recent  measurements  on  similar  samples  in- 
dicate identical  results  for  an  annealing  temperature  of  250°C.  It  is 
possible  that  the  nitrogen  impurities,  left  within  the  implanted  layers,  may 
be  responsible  for  that  part  of  the  increased  thermal  resistance  left  after 
the  annealing.  Argon-  and  nitrogen- implanted  CdS  crystals  have  been  investi- 
gated to  uncover  differences.  ZnTe  implanted  with  A1  has  also  been  scanned. 

It  was  concluded  that  the  comparator  method  of  measurement  offers  good 
potential  as  a useful  tool  in  the  study  of  implanted  layers.  A detailed 

summary  of  the  use  of  the  comparator  to  obtain  thermal  conductivity  data  on 
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II-VI  compounds  has  been  published  during  the  report  period. 

6.6  Thermal  Conductivity  of  Support  Materials 

The  necessity  of  knowing  the  thermal  properties  of  technological 

materials  such  as  nylon  and  greases  has  already  been  mentioned  in  Section  6.2. 

Preliminary  experiments  on  such  materials  have  already  led  to  thermal 
90-96  97 

conductivity  and  heat  capacity  data.  Some  of  this  research  was 

completed  during  the  period  of  this  report.  Although  these  studies  do  not 

involve  II-VI  semiconductor  compounds  directly  they  may  be  considered  as 

valuable  spin-off  connected  with  the  research  since  the  results  have  not 

only  aided  in  the  research  on  II-VI  compounds  but  they  have  also  helped  all 

similar  research.  They  are  mentioned  briefly  in  the  next  few  paragraphs. 

Measurements  of  the  temperature  dependence  of  the  thermal  conductivity 

of  nylon  have  been  carried  out  over  a temperature  range  from  liquid  helium 

temperatures  up  to  350°K.  A maximum  ^max  = 0.371  Wm  ^ was  observed  at 

the  glass-transition  temperature,  T = 330°K.  Because  the  results  gave 

8 

values  which  were  significantly  larger  than  previous  results  at  room  temper- 
atures, and  because  the  data  were  obtained  by  using  the  linear- heat- flow 


method  instead  of  methods  based  on  a radial  flow  of  heat,  it  was  suspected 


that  orientation  anisotropy  in  the  nylon  (possibly  due  to  extrusion  effects) 

could  account  for  the  difference.  Although  measurements  by  a comparator 

method  on  sample  surfaces  cut  perpendicular  and  parallel  to  the  direction 

of  extrusion  failed  to  show  any  anisotropy,  all  of  the  results  were  in 

reasonable  agreement  with  thermal  conductivity  results  obtained  by  the 

linear-heat-flow  method.  Use  was  made  of  X to  compare  thermal  transport 

max . 

of  nylon  with  that  of  other  polymers.  In  particular,  the  dependence  of 

X/X  on  T/T  for  nylon,  natural  rubber,  silicon  rubber,  poly  (methyl 
nicix  • c 

methacrylate)  and  poly  (vinyl  chloride) , was  determined  and  it  was  observed 

that  the  tendency  for  the  curves  to  peak  at  T = was  common.  During  the 

report  period  a paper  on  this  subject  was  presented  at  the  XII  International 

Conference  on  Thermal  Conductivity  in  Birmingham,  Alabama.  The  work  reported 

was  performed  in  collaboration  with  T.  Ashworth  of  the  Physics  Department 

at  the  School  of  Mines  ai  i Technology  in  South  Dakota.  Summaries  of  this 
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work  have  been  published. 

Greases  are  often  used  to  bring  about  improved  thermal  contact  between 

experimental  components  such  as  thermometers,  samples,  heat  sinks,  and 

shields  over  a temperature  range  extending  from  less  than  1°K  up  to  their 

melting  points.  A study  of  the  effectiveness  of  these  materials  has  also 
129 

been  published. 

The  work  described  in  this  paragraph  involves  a coupling  effort  with 
Mr.  J.  Becsey  of  the  ARL  Chemistry  Laboratory  and  workers  in  the  Avionics 


Laboratory  who  were  searching  for  a suitable  heat-sink  liquid.  The  thermal 
conductivity  of  the  organic  dimethylsulf ite  is  of  interest  because  of  its 
potential  use  as  a nonaqueous  battery  solvent  owing  to  its  stability  with 
high  energy  density  electrodes  and  its  ability  to  dissolve  electrolytes. 
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Mahefkey  and  Kreitman  have  noted  that  the  important  thermal  related 

properties  of  battery  performance,  viz.,  cell  cycle  life,  and  cell  capacity, 

as  well  as  electrical  efficiency,  can  be  strengthened  by  operation  at  lower 

battery  temperatures  which  they  achieve  with  the  introduction  of  an  intercell 

heat  pipe.  The  removal  of  waste  heat  from  a battery  is  an  important  design 

characteristic  and  therefore  the  thermal  conductivity  of  the  electrolyte 

solution  is  a significant  factor  in  the  design  and  operation  of  a battery. 

However,  very  little  information  exists  on  the  thermal  properties  of  dimethyl 

sulfite.  In  order  to  help  fill  this  gap,  a thermal  comparator  was  used  to 

obtain  measurements  on  the  room-temperature  thermal  conductivity  of  a liquid 

sample  of  dimethylsulf ite  using  a liquid-cell  adapter  and  the  basic  thermal 
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comparator  apparatus  described  elsewhere.  Mr.  Becsey  supplied  the 

dimethylsulfite  sample  and  a number  of  the  calibration  standards,  i.e.,  H2O, 

glycerol,  methanol,  and  CCl^.  The  value  obtained  for  the  thermal  conductivity 

of  dimethylsulfite  was  0.00178  Wcm  * . For  more  detailed  information  the 
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reader  is  referred  to  the  published  article. 

6 . 7 Some  Notes  on  the  Comparator  Method 

We  have  used  the  comparator  method  on  many  other  crystals  and,  aside 
from  its  advantages,  there  are  definite  problems  that  the  user  should  be  aware 
of.  One  problem  is  the  necessity  of  sufficient  and  clean  calibrating  standards. 
This  is  particularly  true  when  measuring  CdTp  which  has  a conductivity 
between  stainless  steel  and  armco  iron.  Because  of  a bend  in  the  calibration 
curve  in  this  region  more  calibrating  standards  are  necessary.  Another 
difficulty  with  some  standards, such  as  armco  iror*  is  the  oxidation  which 
affects  its  surface.  It  may  be  possible  to  replace  such  standards  with  II-VI 
compounds.  Our  results  show  that  the  reproducibility  of  readings  on  CdS,  for 
example,  is  much  higher  than  that  for  the  iron  or  the  stainless  steel. 
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The  thermal  comparator  appears  to  yield  its  highest  reproducibility 
when  measuring  plastics. 

6.8  Linear-Heat-Flow  Measurements  on  ZnSe  and  Type  316  Stainless  Steel 


With  a modified  form  of  the  apparatus  used  to  measure  the  thermal 
conductivity  of  nylon  and  greases  (described  in  reference  95),  a number  of 
linear-heat-flow  measurements  have  been  carried  out  on  bar  samples  of  ZnSe 
and  type  316  stainless  steel,  from  room  temperatures  down  to  pumped  nitrogen 
temperatures.  The  former  material  was  supplied  by  the  Laser  Window  Group 
of  the  Materials  Laboratory.  Mr.  D.  Locker  of  the  Solid  State  Physics 
Laboratory  collaborated  in  the  measurements. 

The  results  of  our  linear-heat-flow  measurements  on  the  ZnSe  bar  (dimensions 
0.62  cm  x 5.32  cm  x 0.28  cm)  yield  room-temperature  values  differing  (by  more 
than  our  limits  of  error)  with  the  literature  values  of  ZnSe.  Although  the 
addition  of  insulating  powder  and  increased  vacuum  were  observed  to  decrease 
the  disparity,  a significant  variance  remained.  In  order  to  test  the  accuracy 
of  the  modified  linear-heat-flow  apparatuj  a bar  of  type  316  stainless  steel 
was  machined  and  measured  in  place  of  the  ZnSe.  In  addition,  a disc  of  the 
same  316  material  was  machined  and  polished  to  the  dimensions  of  the  316 
calibration  standard  of  the  thermal  comparator  instrument.  The  results  of 
comparator  measurements  on  the  316  sample  were  found  to  be  identical  to  those 
on  the  comparator  316  calibrating  standard,  yielding  a room  temperature  conductivity 
of  approximately  0.14  Wcm  ^ for  the  316  test  bar  and  thus  validating  its 
use  as  a means  of  checking  the  modified  linear-heat-flow  apparatus.  Subsequent 
linear-heat-flow  measurements  on  the  316  test  bar  yielded  thermal  conductivity 
data  slightly  larger  than  its  accepted  value.  In  other  experiments,  insulating 
powder  (to  reduce  radiation  losses)  and  increased  vacuum  have  also  been  observed 
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to  reduce  the  discrepancy  in  the  stainless  steel  data.  Although  che  experiment 
had  to  be  terminated  in  order  to  attend  to  the  other  investigations  (below) 


it  appeared  that  the  preliminary  results  of  the  measurements  indicated  a ZnSe 

103 


thermal  conductivity  which  was  slightly  higher  than  that  found  by  Slack 
in  Table  6-3. 

6 . 9 Debye  9 and  Heat  Capacity  Studies  in  Il-VI  Compounds 

The  II-VI  compounds  have  aroused  a good  deal  of  interest  over  the  past 
few  years  for  a variety  of  reasons.  In  order  to  analyze  the  vibrations  in  these 
compounds  and  to  analyze  properties  of  the  vibrational  frequency  spectra, 
knowledge  of  the  specific  heat,  thermal  expansion,  and  elastic  constants  are 
required,  ideally  over  a wide  range  of  low  temperatures.  An  examination  of  the 
available  thermodynamic  data  for  the  oxides,  sulfides,  selenides,  and  tel1 ur ides 
of  zinc,  cadmium,  and  mercury  reveals  the  existence  of  a -large  number  of  gaps. 
Only  for  zinc  oxide  and  zinc  sulfide  do  sufficiently  complete  sets  of  data 


exist  to  justify  thermodynamic  analysis,  and  even  in  these  cases  the  analyses 
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can  only  be  embarked  upon  with  reservations. 

Furthermore, theoretical  studies  of  low-temperature  specific  heat  anomalies 
arising  from  simultaneous  effects  of  exchange  and  magnetic  fields  had  been 
■initiated  to  study  factors  which  influence  impurity  peaks  in  the  1 ow -temperature 
specific  heat  of  samples  in  zero  and  non-zero  magnetic  fields.  One  conclusion 
was  that  the  simultaneous  presence  of  exchange  interdictions  and  external  magnetic 
fields  on  doped  crystals  may  lead  to  two  peaks  in  the  low-temperature  dependence 
of  the  heat  capacity.  These  two  peaks  may  best  be  seen  when  tb®  ratio  of  the 
magnetic  energy  to  the  exchange  energy  is  close  to^ but  not»_qual  to, certain 
integers,  the  number  of  which  varies  directly  as  the  spin  of  the  magnetic  atoms 


in  a pair  system.  Further  information  may  be  obtained  from  the  work  whi'ch  was 
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published  during  the  report  period.  The  application  or  test  of  these  ideas 
rests  on  attempts  approximate  a dilute  state  of  magnetic  material  by  adding 
magnetic  atoms  to  a diamagnetic  lattice,  and  descriftTng  the  resulting  magnetic 
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cooperative  phenomena  by  the  cluster  method,  which  considers  the  magnetic 
impurity  atoms  in  various  nearest  neighbor,  next-nearest  neighbor,  and/or 


higher  order  clusters. 

In  order  to  satisfy  a need  for  more  experimental  data  on  II-VI  compounds 

and  on  low-temperature  heat-capacity  anomalies,  work,  had  been  initiated  on 

CdS.  During  the  report  period  work  on  the  specific  heat  and  Debye  9 of  CdS:0.45% 

Mn  was  completed.  Experimental  heat-capacity  measurements  permitted  values 

of  Debye  0 to  be  tabulated  down  to  3.32°K.  The  lo^-temperature  limiting  value 

of  (T  = 0)  = 216  t 4°K  was  found  to  correspond  satisfactorily  with 

an  experimental  value  of  215°K  obtained  by  Gerlich^^^  from  low-temperature 

elastic  constant  data  on  CdS.  This  result  differed  by  7.4%  with  the  theoretical 

value  of  200°K  obtained  by  Goel  and  Singh'*'^^  based  on  a simple  model  involving 

a radial  nearest-neighbor  approximation.  A theoretical  value  of  263°K 
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published  by  Nusimovici  et.al.,  was  found  to  be  19.2%  different.  Additional 
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information  may  be  obtained  from  our  published  work. 


7.0  OPTICAL  INVESTIGATIONS  IN  II-V1  COMPOUNDS  AND  RELATED  MATERIALS 


Optical  absorption  studies  can  be  used  to  give  important  information. 

Although  the  chief  interest  in  the  effects  produced  by  irradiation  of  semi- 
conductors is  generally  agreed  to  lie  in  the  area  of  electrical  properties, 
studies  of  the  optical  absorption  (as  well  as  other  spectroscopic  methods) 
can  give  detailed  information  related  to  the  microscopic  character  of  the 
defects.  Irradiation  normally  produces  a variety  of  new  absorption  bands 
which  can  be  divided  into  "Localized  Vibrational  Bands"  and  "Electronic  Excitations" 
The  identification  of  the  defects  associated  with  the  bands  could  help  in 

solving  many  unanswered  questions.  Such  studies  have  led  to  information  about 
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the  local  symmetry  of  particular  defects.  Once  identified,  the  band 
can  serve  as  a tool  for  monitoring  the  damage.  Optical  spectra  can  also  give 
information  about  the  excited  states  to  which  a transition  is  made. 

7. 1 Raman  Scattering  in  ZnTe 

Work  in  this  area  of  research  was  initiated  in  order  to  carry  out 
electron-spin-flip  scattering  measurements  and  study  the  associated  physical 
parameters  of  II-VI  semiconductors.  In  one  experiment,  measurements  were  carried 
out  on  a ZnTe  crystal  at  4.2°K  in  various  magnetic  fields  up  to  89  kG,  approximately 
At  least  3 clear  spin-flip  peaks  were  observed  to  occur  at  short  distances  from 
the  laser  line  in  the  spectra, for  slit  widths  less  than  50  p.  Other  experi- 
mental work  performed  within  the  report  period  has  been  at  liquid-helium 
temperatures  and  in  magnetic  fields  up  to  90  kG  with  scattered  light  from  an 
argon  laser  line.  This  work  has  disclosed  (possibly  for  the  first  time)  what 
appears  to  be  second-,  third-,  and  possibly  fourth-order  stokes  and  antistokes 
lines.  Time  was  also  devoted  to  aligning  and  replacing  components  in  a 
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krypton  laser  system.  In  addition,  the  vacuum  system  associated  with  this 
experimentation  area  was  modified.  Part  of  this  effort  involved  setting  up 
a 69  kG  superconducting  magnet  system  and  designing  metal  and  glass  parts 
to  complement  and  improve  various  experiments  with  a 100  kG  magnet  system. 

7.2  Infrared  Absorption  in  ZnS  and  ZnSe 

Studies  in  infrared  absorption  were  initiated  in  order  to  examine 
the  near  and  far-infrared  absorption  of  ZnSe:Co,  ZnSe:Al,  ZnS:Co,  GaAs,  and 
certain  other  crystals  with  alkali-metal  dopings  as  a function  of  temperature, 
including  cryogenic  temperatures,  with  the  aim  of  sorting  out  and  identifying 
host-lattice  native  defects  and  impurity  defects  having  potential  possibilities 
leading  to  the  production  of  better  infrared  detectors  from  the  large-band- 
gap  compounds. 

Taken  with  a 60-inch  grating  spectrometer,  absorption  measurements  in 
the  near-  and  far-infrared  on  single  crystals  of  pure  and  aluminum  doped  ZnSe, 
have  revealed  the  presence  of  multiphonon  processes  along  with  localized 


vibrational  modes,  the  latter  arising  from  the  A1  doping.  ZnSe  interference 
phenomena  displayed  at  ~30p  from  crystal  temperatures  of  300,  100  and  40°K 


make  possible  a determination  of  the  index  of  refraction. 

Near-infrared  measurements  on  undoped  ZnS  samples  and  4 ZnStCo  samples  have 

revealed  large  absorptions  in  the  regions  1.30  to  1.86  p,  and  2.4  to  3.0  p, 

4 4 7 

approximately,  corresponding  to  T^(F)  and  T9(F)  levels  of  cobalt  (d  ) in 

cubic  ZnS  with  a crystalline  field  parameter,  , of  approximately  375.  At 

liquid-helium  temperatures  sharp  absorption  lines  appear  in  both  the  3 and  1 p 

grating  spectra  of  those  crystals  having  significant  cobalt  impurity.  T>ata 

have  been  obtained  on  various  experimental  samples  in  the  wavelength  regions 

offered  by  the  1,  3,  18,  30  and  45  p gratings,  and,  after  zero-corrections  are 

made,  the  results  are  in  approximate  agreement  with  those  cases  where  published 
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results  appear  in  the  literature.  Water  vapor-lines  were  used  to  check  the 
calibration  of  the  spectrometer. 

7.3  Infrared  Absorption  in  ot-HgS 
7.3.1  Introduction 

Mercury  sulfide  is  a wide  bandgap  semiconductor  which  is  of 

considerable  interest  because  its  unique  properties  offer  excellent  potential 

for  non-linear  optics  and  laser-technology  devices.  It  is  the  most  optically 

active  of  all  known  mineral  compounds,  possessing  the  greatest  birefringence, 

a remarkable  rotary  power,  and  strong  piezoelectric  properties. 
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Sapriel  has  reported  that  a-HgS  has  the  highest  acousto-optical  figure 

of  merit  of  any  known  crystal,  making  it  an  efficient  material  for  acousto- 

optical  applications,  such  as  transducers,  light  modulators,  deflectors,  delay 

lines,  etc.  Ultrasonic  attenuation  and  the  elastic  and  photoelastic  constants 
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in  a-HgS  have  been  measured.  Sapriel  and  Lancon  have  pointed  out  that, 

since  the  electromechanical  coupling  coefficients  of  cinnabar  are  approximately 

two  times  larger  than  those  of  quartz,  piezoelectric  a-HgS  transducers, bonded 

on  a parallelpiped  cinnabar  crystal  would  realize  perfect  accoustical  matching 

and  therefore  yield  wide  bands  for  acousto-optical  deflection.  Thus  a-HgS 

should  serve  well  as  a good  deflectior  of  laser  light.  For  example,  if  one 

assumes  a frequency  bandwidth  of  250  MHz  for  a transducer  made  of  cinnabar, and 

an  acoustic  beam  1/2  mm  square,  it  may  be  shown  that  only  210  mW  of  acoustic 

power  is  necessary  to  deflect  all  of  the  incident  light  of  a He-Ne  laser  beam 

into  62,500  different  positions  with  an  access  time  of  1 ps. 
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Measurements  on  the  red -'colored  a-phase  of  synthetic  and  good  natural 
crystals  indicate  a transparency  of  more  than  50%  between  0.6  and  14  p at  300°K. 
Residual  absorption  in  the  transparent  region  is  known  to  be  less  tha  1 cm 
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for  good  quality  natural  crystals.^"’  Large  indices  of  refraction  have 
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been  measured  over  the  range  0.62-11  u.  The  proximity  of  the  upper  limit 

of  the  transparent  region  to  the  6328$  He-Ne  laser  line,  and  the  large  temper- 

-4 

ature  dependence  of  the  HgS  absorption  edge  (~9xl0  eV/C),  has  already  been 
exploited*^  for  thermal  modulation  of  the  laser  line  (by  the  thermoabsorption 
of  cinnabar). 

Photoelectronic  properties  of  a-HgS  have  been  studied,  chiefly  in  the  more 

available,  good  quality,  natural  crystals  and  these  studies  reveal  a close 
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similarity  to  those  of  other  II-VI  compounds.  Some  results  in  synthetic 
crystals  include  a peak  photoconductivity  at  ~6 000X  (300°K)  for  HgS  (~10^ficm) ; 
the  existence,  within  the  energy  gap,  of  localized  levels  of  electrons  (0.62  eV) 
and  holes  (1.08  eV)  from  the  conduction  band;  and  thresholds  (0.67  eV  and  1.08  eV) 


for  optical  quenching.  In  other  studies  the  conductivity  of  illuminated 

natural  crystals  has  been  observed  to  increase  by  6 orders  of  magnitude. 

Furthermore  the  mobility  in  a direction  perpendicular  to  the  c-axis  has  been 

determined'*'"’'*'  to  originate  from  electrons  at  10  cm^  V ^s  ^ while  others'*'"’^ 

2 -1  -1 

have  noted  a photo-Hall  mobility  of  38  cm  V s in  a crystal  which  may  not 
have  been  oriented. 

The  strong  non-linear  properties  of  cinnabar  are  also  important  for  integrated 


optical  circuits,  parametric  oscillator  upconversion,  difference-frequency 
generation,  tunable  lasers,  and  other  applications.  For  example,  Martin  and 


Thomas  have  suggested  that  the  stimulated  Raman  effect  be  combined  with  the 
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non-linear  electro-optic  effect  in  HgS  in  order  to  conveniently  generate 

pulses  of  radiation  between  9 and  15  h . They  utilize  the  non-linear  electro- 

3+ 

optic  effect  to  generate  a beat  frequency  between  a Nd  /glass  laser  and  the 
laser-activated,  stokes  Raman  emission.  The  difference  frequency  is  equal  to 


the  molecular  vibrational  frequency  of  the  Raman  process. 


* 


The  successful  future  utilization  of  all  of  these  devices,  and  the 

exploration  and  study  of  these  and  additional  properties  of  a-HgS,  rests 

heavily  on  the  growth  of  good  quality  single  crystals,  which,  until  recently, 

has  been  very  difficult.  Below  we  report  optical  measurements  which  help 

to  confirm  the  yield  of  good  cinnabar  crystals.  These  results  form  part 

of  a presentation  at  the  Conference  on  Optical  Properties  of  Highly  Transparent 

Solids  at  Waterville  Valley,  New  Hampshire,  in  1975,  the  proceedings  of  which 
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have  been  published  by  Plenum  Press.  A preliminary  report  of  this  work 
was  reported  earlier. 

7.3.2  Results 

Some  of  the  larger  cinnabar  crystals,  grown  by  the  iodine 
141 

vapor-transport  methods  , were  used  as  samples  for  optical  investigations 
in  the  visable  and  infrared  spectral  regions.  A typical  measured  value  of 
crystal  thickness  was  330  p.  In  one  study,  a 2-meter  spectrograph  was  employed 
to  take  photographic  plates  (4  8/mm)  of  the  edge  absorption  of  a number  of 
crystals  on  which  a tungsten  source  was  focused.  In  this  arrangement  the 
crystals  were  mounted  on  a rod  within  a fixed  temperature  dewar  and  measurements 
were  carried  out  at  room  temperature  and  at  4.2°K.  A densitometer  analysis 
of  the  photographic  plates  enabled  values  of  the  wavelength  at  zero  absorption 
to  be  tabulated  with  the  exposure  times  and  this  led  to  an  extrapolation  of 
the  wavelength  corresponding  to  zero-exposure  time  and  a value  of  6115  X for 
the  approximately  | | absorption  edge  at  300°K  with  a decrease  of  approximately 
2.0  X/K  on  cooling  to  4.2°K.  Additional  transmission  data  taken  on  a Cary 
spectro-photometer  were  used  to  verify  this  value  for  the  edge.  A value  of 


the  reflection  at  0.63  was  calculated  as  26%.  The  observed  temperature  dependence 
is  in  rough  agreement  with  a literature  value^^  of  1.8  X/K,  the  difference  being 
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Ei 


easily  attributed  to  the  pronounced  dichroism  of  the  a-HgS  edge.  In  fact, 
because  of  the  small  size  and  delicate  nature  of  the  crystals,  use  was  made 
of  this  anisotropy  to  check  the  orientation  of  the  crystal  after  mounting. 
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Although  we  did  not  thoroughly  search  for  it  at  very  low  absorption  levels, 

we  did  not  observe  any  discrete  structure  near  the  edge  at  4.2°K.  Typical 

9 

values  of  resistivity  for  our  samples  were  approximately  2 x 10  ftcm  which  may 

12  149 

be  compared  to  10  ftcm  obtained  by  Roberts,  Lind  and  Davis  for  their  undoped 

synthetic  cinnabar.  Their  high  resistivity  crystals,  of  mm  dimensions,  were 

grown  at  temperatures  just  below  the  polymorphic  transition  temperature  by 

a reaction  of  the  HgS  elements  in  sealed  quartz  ampoules  containing  200-400  mm 

Hg  pressure  at  room  temperature. 

Infrared  measurements  were  performed  with  the  crystals  at  temperatures 
of  32,  92,  and  300°K  in  a metal  cryostat  and  mounted  on  a cold  finger,  instrumented 
with  Ge,  Pt,  and  carbon  thermometry.  The  tail  and  barrel  of  the  devar  were 
enclosed  within  a vacuum  spectro-photometer  having  a focal  length  of  60  inches. 

The  spectrometer  has  Czerny-Turner  optics  and  utilizes  interchangeable  7"  x 9", 
kinematically  mounted  reflection  gratings,  blazed  at  1,  3,  9,  18,  30,  45,  60  and 
100  u with  associated  spectral  reststrahlen  filter  plates,  filters,  choppers, 
and  I.R.  sources  together  with  a Golay  cell  detector.  Sapphire,  ZnRe  and 
polyethelene  dewar  windows  were  employed  for  the  low-temperature  measurements 
and  at  30  u a typical  value  of  the  inverse  linear  dispersion  was  experimentally 
computed  to  vary  between  0.028-0.030  u/mm  for  a range  of  grating  angles. 

Measurements  of  the  infrared  transmission  were  carried  out  in  vacuum  over 
a wavelength  range  which  extends  from  the  fundamental  absorption  to  beyond 
the  restrahlen  edge  beginning  at  approximately  27  u.  Some  typical  experimental 
data  representing  some  of  the  more  prominent  absorption  peaks  which  we  have 
observed  are  presented  in  Table  7-1  for  comparison  with  the  work  of  others. ^ ^ 
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Table  7-1  Low  temperature  I-R  data  on  a-HgS  phonon  frequencies! 


E 

A 

E 

l 

1 

T 

L 

T 

T 

L I 

TOO 

Source 

285 

290 

334 

340 

348 

32 

Present  work,  ARL 

283 

290 

340  ' 

344 

350 

300 

152 

Barcello,  et.  al. 

277 

287 

333 

338 

346 

90 

153 

Zallen,  et.  al. 

284 

345 

77 

154 

Riccius,  et.al. 

280 

335 

343 

300 

Poulet,  et.  al.^^ 

From  the  data  comparisons  summarized  above  it  may  be 
concluded  that  the  optical  measurements  verify  the  growth  of  a-HgS  crystals  and 
that  the  iodine  transport  method  and  its  associated  techniques  merit  further 
study  and  development.  Furthermore  it  appears  that  when  the  cinnabar  growth 
methods  have  been  developed  to  the  extent  that  large  single  crystals  are 
available  in  sufficient  supply,  that  a-HgS,  because  of  its  unique  properties, 
will  be  in  great  demand  as  a device  material. 

7 . 4 Characterization  of  GaAs  Samples  by  Infrared  Spectroscopy 
7.4.1  Introduction 

The  need  to  achieve  a more  complete  characterization  and 
identification  of  defects  in  compound  semiconductors  is  an  important  and  urgent 
response  to  the  demands  for  improved  solid-state  microwave  and  optoelectronic, 
III-V  semiconductor  devices.  It  has  now  become  widely  accepted  that  the  goals 
of  long-term  stability,  low  cost,  and  reliable  predictions  of  performance  of 
these  devices  can  be  met  only  with  a thorough  materials  characterization  program 
which  will  not  only  provide  feedback  of  the  usual  characterization  information 
to  the  materials  growth  technology,  but  will  also  investigate  and  utilize  new 
procedures  for  materials  characterization  and  combine  these  with  the  established 
methods.  Infrared  spectroscopy,  utilizing  observations  of  photoconductivity. 


i 

i 


i 


reflection,  luminescence,  and  transmission,  are  valuable  measurement  techniques 
that  can  be  used  in  order  to  provide  new  material  characterization  data  which 
can  be  recognized  and  evaluated  as  having  to  do  with  the  materials'  intrinsic 
properties  and  can  be  separated  from  those  parameter  changes  which  are  introduced 
by  the  effects  of  device  processing.  Infrared  spectroscopy  cannot  only  detect 
the  presence  of  electrically  inactive  impurities  such  as  oxygen  or  carbon, 
but  with  proper  resolution  (accomplished  by  cryogenic  techniques),  it  can 
also  distinguish  between  the  kinds  of  impurites  present. 

7.4.2  Infrared  Photoconductivity  of  GaAs:0 

The  results  of  a series  of  characterization  experiments  on  a 

GaAs:0  sample  demonstrated,  as  erroneous,  a set  of  unusual  conclusions  drawn 
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by  H.  L.  Stocker  in  his  publication  entitled,  "Probable  Observation  of 
Unusual  Local  Vibrational  Modes  in  the  Oscillatory  Photoconductivity  of  Semi- 
insulating  GaAs."  Our  work  was  concerned  chiefly  with  the  characterization 
of  sample  76-3,  a GaAs:0  sample  which  was  reported  to  be  an  oxygen-doped  substrate 
material. 

The  results  of  photoconductivity  observations  at  300°K  under  a 300  V D.C. 
bias  indicated  a gradual  increase  in  photoconductivity  response  beginning  at 
approximately  0.7  eV  and  extending  to  a maximum  of  about  3nA  at  the  silicon- 
filter  cut  off  near  1.0  eV.  Difficulties  experienced  at  nitrogen  temperatures 
with  crystal  breakdown-voltage  oscillations  were  eliminated  when  the  bias  voltage 
was  decreased  to  50  V D.C.  At  300°K  with  the  new  50  V bias  (and  without  the 
silicon  filter),  the  bandgap  energy  was  observed  to  be  represented  by  2 peaks 
separated  by  33  meV,  and  occurring  at  1.42  and  1.39  eV. 

Some  sharp  structure  In  the  photoconductive  response  near  the  bandgap  was 
traced  to  gas  pressure  and  eliminated  by  vacuum  improvements.  Since  silicon 
impurity  has  been  observed  in  the  analysis  of  GaAs  samples  and  since  silicon  is 
known  to  occur  as  an  acceptor  with  an  energy  of  0.03  eV,  it  is  possible  that 
one  of  these  levels  near  the  bandgap  is  related  to  the  presence  of  silicon. 
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at  nitrogen  temperatures  and,  unlike  other  GaAs  samples  measured  previously, 
the  band  remained  even  at  temperatures  below  30°K.  However,  at  these  cryogenic 
temperatures  there  was  some  improvement  in  resolution,  suggesting  a central 
valley  and  hinting  at  the  possibility  that  the  band  was  predominantly  due  to 
two  chief  levels,  with  transition  energies  of  0.80  and  0.89  eV.  It  may  be 
concluded  from  these  and  other  measurements  that:  1)  the  sample  contains 

chromium  in  substantial  amounts;  2)  other  physical  or  chemical  impurities 
must  also  exist;  and  3)  that  the  quality  of  this  crystal  is  relatively  poor. 

In  order  to  become  more  familiar  with  the  characteristics  of  the  oxygen  level 
it  would  be  beneficial  to  investigate  GaAs  material  containing  only  oxygen 
as  a significant  dopant. 

7.4.3  I-R  Photoconductivity  of  GaAs:Cr 

Near-infrared  photoconductivity  spectra  have  been  taken  on  a 
number  of  GaAs:Cr  samples  at  various  temperatures  and  it  has  been  observed  that 
there  are  significant  differences  between  groups  of  crystals  from  different 
manufacturers  as  well  as  important  differences  between  the  results  of  individual 
samples  within  each  group.  The  results  on  one  group  of  samples  are  summarized 
in  Table  7-2  which  lists  the  maximum  photoconductivity  and  the  threshold 
energy  observed  for  these  materials  at  temperatures  ranging  between  300°K  and 
25°K.  In  Table  7-2  the  samples  are  grouped  according  to  their  photoconductive 
response  at  various  experimental  temperatures.  The  four  Cr-doped  substrates 
gave  similar  data  in  that  each  sample  displayed  an  infrared  photoconductive 
spectrum,  the  intensity  of  which  decreased  with  a decrease  of  sample  temperature. 

In  addition,  on  cooling  to  approximately  nitrogen  temperature,  each  spectrum  resolved 
itself  into  a single  peak  due  to  the  chromium  impurity.  At  approximately 
helium  temperature  the  single  peak  had  decreased  slightly  in  intensity  and  had 
shifted  slightly  to  a higher  energy. 
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This  crystal  was  received  in  a broken  and  over-etched  condition. 


Differences  between  the  results  of  the  four  Cr-doped  substrate  samples 


were  chiefly  in  the  intensity  of  the  observed  photoconductive  response. 

Sample  9B  and  broken  Sample  54B  yielded  the  greatest  intensity  at  all  temperatures 
but  particularly  so  at  cryogenic  temperatues  where  the  ratios  of  the  photo- 
conductive  maxima  of  Samples  9B  and  54B  to  those  of  Samples  4B  and  49B, 
respectively,  grew  from  more  than  5 at  temperatures  near  100°K  to  greater 
than  10  at  temperatures  close  to  30°K.  Comparing  9B  with  54B  one  might  conclude 
from  their  cryogenic  spectra  that  they  have  approximately  the  same  amount 
of  chromium  and  that  it  is  the  only  impurity.  Perhaps  9B  has  slightly  more 
Cr  but  the  fact  that  Sample  54B  was  broken  and  overetched  may  invalidate 
this  conclusion. 

Although  Sample  49B  exhibited  the  least  photoconductive  response,  a chromium 
peak  or  band  was  easily  seen  at  cryogenic  temperatures.  If  one  uses  the 
intensity  of  the  chromium  peak  as  an  indication,  it  appears  that  there  is 
significantly  more  chromium  in  Samples  9B  and  54B  than  in  Sample  4B  and  49B. 

Sample  47B  gave  interesting  and  quite  different  results.  It  was  known 
to  be  an  LPE  buffer  on  a semi- insulating  substrate.  On  cooling  to  105°K  the 
intensity  of  the  photoconductive  response  increased  over  three  orders  of 
magnitude  as  the  broad  300°K  band  resolved  itself  into  a narrow  band  due  to 
chromium.  A further  increase  in  intensity  was  noted  when  the  sample  was 
cooled  to  approximately  30°K. 

Furthermore  the  intensity  of  the  maximum  of  the  photocurrent  response 
at  cryogenic  temperatures  was  so  large  that  the  spectrometer  slits  had  to  be 
closed  to  1/2  mm.  Even  then,  at  temperatures  near  30°K  and  with  a 300  volt 
D.C.  bias,  corresponding  to  a maximum  photoconductivity  of  approximately  6 pA, 
current  oscillations  and  breakdown  effects  were  observed.  To  avoid  tbcse 
effects  the  bias  voltage  was  reduced  to  50  volts.  The  results  for  Sample  47B 
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